Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,797,377 (T)5→4 coding (517/960 nt) nrdF → class 1b ribonucleoside‑diphosphate reductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,797,3730T.94.7% 212.8 / 8.7 57coding (513/960 nt)nrdFclass 1b ribonucleoside‑diphosphate reductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base . (13/41);  total (14/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

ATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAG  >  NZ_CP009273/2797285‑2797458
                                                                                        |                                                                                     
aTCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAtt                                                                                      >  2:942029/1‑90 (MQ=255)
aTCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAtt                                                                                      >  1:550299/1‑90 (MQ=255)
     ggTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttcca                                                                                 >  2:1011744/1‑87 (MQ=255)
      gTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttccag                                                                                >  2:716282/1‑86 (MQ=255)
       tgatgaTCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGc                                                                              <  2:550299/90‑1 (MQ=255)
        gatgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGcc                                                                             >  2:1561013/1‑90 (MQ=255)
            atCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCgg                                                                         <  1:716282/90‑1 (MQ=255)
              ccGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGaa                                                                       <  2:1563230/90‑1 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  1:101636/1‑90 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  2:1319808/1‑90 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  1:863049/1‑90 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  2:916483/1‑90 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTCCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  2:470744/1‑90 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:741124/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:1481524/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:327163/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:53096/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  1:1546819/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:80861/90‑1 (MQ=255)
                      aaaGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGAcc                                                               <  2:1445175/90‑1 (MQ=255)
                         gAAAATCGCCAGTCTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAAt                                                            <  1:1372111/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGGCCAATACcgc                                                       <  1:982125/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  1:930294/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  1:63631/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:1428604/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:16446/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:908344/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:275569/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  1:798158/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:896625/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:660417/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  1:916483/90‑1 (MQ=255)
                                   aGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGAcc                                                  >  1:944033/1‑90 (MQ=255)
                                         tttCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATcc                                            >  2:743582/1‑90 (MQ=255)
                                              tGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTg                                       >  1:1259850/1‑90 (MQ=255)
                                                    ttttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATc                                 >  1:352311/1‑90 (MQ=255)
                                                     tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                                <  2:1277338/90‑1 (MQ=255)
                                                     tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                                <  2:83861/90‑1 (MQ=255)
                                                     tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                                <  1:1561013/90‑1 (MQ=255)
                                                            tCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGa                            <  2:608261/88‑1 (MQ=255)
                                                            tCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGa                            >  1:608261/1‑88 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                                <  1:1371663/80‑1 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                                >  2:1371663/1‑80 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                      >  1:1279539/1‑90 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                      <  1:1319808/90‑1 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                      <  2:418572/90‑1 (MQ=255)
                                                                     gTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCa                <  2:1050539/90‑1 (MQ=255)
                                                                     gTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCa                <  1:949892/90‑1 (MQ=255)
                                                                      tttCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCa                <  2:707107/90‑1 (MQ=255)
                                                                      tttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc               <  2:759893/90‑1 (MQ=255)
                                                                       ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg              <  1:579747/90‑1 (MQ=255)
                                                                       ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg              <  2:124172/90‑1 (MQ=255)
                                                                       ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg              <  1:1494487/90‑1 (MQ=255)
                                                                       ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg              <  2:961138/90‑1 (MQ=255)
                                                                         cTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGt            >  2:687835/1‑90 (MQ=255)
                                                                          tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt           <  1:955879/90‑1 (MQ=255)
                                                                          tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt           <  2:944033/90‑1 (MQ=255)
                                                                          tGGCCGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt           <  2:372954/90‑1 (MQ=255)
                                                                           ggCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTta          <  2:40449/90‑1 (MQ=255)
                                                                                   aTGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAg  <  2:558226/90‑1 (MQ=255)
                                                                                   aTGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAg  <  2:1259850/90‑1 (MQ=255)
                                                                                        |                                                                                     
ATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAG  >  NZ_CP009273/2797285‑2797458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG  >  NZ_CP009273/2797280‑2797459
                                                                                                 |                                                                                  
ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCC                                                                                  >  SRR3722078.1186056/1‑100 (MQ=60)
  ATTATCGCGGTGATGATAAGATGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGC                                                                               <  SRR3722078.83366/100‑1 (MQ=60)
           GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAG                                                                      >  SRR3722078.102841/1‑100 (MQ=60)
           GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAG                                                                      >  SRR3722078.874602/1‑100 (MQ=60)
                 ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACC                                                                <  SRR3722078.725104/100‑1 (MQ=60)
                      CTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATAC                                                           <  SRR3722078.1572021/100‑1 (MQ=60)
                              GAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACC                                                   >  SRR3722078.957207/1‑100 (MQ=60)
                              GAAAATCGCCAGTCTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACC                                                   <  SRR3722078.1394217/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGGCCAATACCGCGGACCTGATC                                              <  SRR3722078.996303/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                              <  SRR3722078.64380/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                              <  SRR3722078.808659/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                              <  SRR3722078.928931/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                              <  SRR3722078.943109/100‑1 (MQ=60)
                                         GTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTG                                        >  SRR3722078.1279822/1‑100 (MQ=60)
                                               TTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATC                                  >  SRR3722078.356615/1‑100 (MQ=60)
                                                       TCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATG                           >  SRR3722078.615289/1‑100 (MQ=60)
                                                          TTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGC                       <  SRR3722078.1586418/100‑1 (MQ=60)
                                                          TTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGC                       <  SRR3722078.1393753/100‑1 (MQ=60)
                                                          TTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGC                       >  SRR3722078.1299893/1‑100 (MQ=60)
                                                                    ATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGT             <  SRR3722078.1340901/100‑1 (MQ=60)
                                                                          GTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTAC       <  SRR3722078.963218/100‑1 (MQ=60)
                                                                            TTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAT     <  SRR3722078.586509/100‑1 (MQ=60)
                                                                            TTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAT     <  SRR3722078.1518852/100‑1 (MQ=60)
                                                                               TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTATATAGG  <  SRR3722078.969381/100‑1 (MQ=60)
                                                                                                 |                                                                                  
ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG  >  NZ_CP009273/2797280‑2797459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: