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 F72 I1 R1 14 98.9 3382490 97.5% 3297927 140.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA97.4% 256.7 / ‑3.6 79P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (32/43);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184400‑4184681
                                                                                                                                               |                                                                                                                                          
atatCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATAAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCATCTgg                                                                                                                                       >  1:412258/1‑149 (MQ=255)
 tatCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       <  1:219207/148‑1 (MQ=255)
     tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                                  >  1:360943/1‑149 (MQ=255)
      ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                 >  3:66907/1‑149 (MQ=255)
         ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                              <  7:309874/149‑1 (MQ=255)
          gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                             <  7:301704/149‑1 (MQ=255)
           ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                            <  4:66907/149‑1 (MQ=255)
           ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                            <  2:360943/149‑1 (MQ=255)
             acaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGCGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                           <  5:51728/146‑1 (MQ=255)
             aaaCGCCGTATCCAGACTGTTGAAAAGATGGCAGGGCGTCACGGTAACTAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                           <  4:213556/148‑1 (MQ=255)
              aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                         <  4:360815/149‑1 (MQ=255)
              aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCAGCGCGTACCGTCTCg                                                                                                                          >  8:216851/1‑148 (MQ=255)
                cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                       <  5:232866/149‑1 (MQ=255)
                       tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATc                                                                                                                >  6:402613/1‑149 (MQ=255)
                          aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                         >  2:281643/1‑137 (MQ=255)
                          aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                         <  1:281643/137‑1 (MQ=255)
                          aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                             <  6:361358/149‑1 (MQ=255)
                          aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAAATCGGt                                                                                                             <  7:23368/149‑1 (MQ=255)
                             cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                          <  5:140699/149‑1 (MQ=255)
                               ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                             >  1:344193/1‑144 (MQ=255)
                                gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                         <  4:106643/147‑1 (MQ=255)
                                       gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       >  7:45054/1‑110 (MQ=255)
                                       gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       <  8:45054/110‑1 (MQ=255)
                                       gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  7:286123/149‑1 (MQ=255)
                                       gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                >  1:240800/1‑149 (MQ=255)
                                       gATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  4:310822/149‑1 (MQ=255)
                                        aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                               <  2:344193/149‑1 (MQ=255)
                                           gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                            >  4:124596/1‑149 (MQ=255)
                                            cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                           <  1:177717/149‑1 (MQ=255)
                                             aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                          >  4:20680/1‑149 (MQ=255)
                                                 cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                      <  5:331828/149‑1 (MQ=255)
                                                     aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                                   <  3:338316/148‑1 (MQ=255)
                                                         tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                              <  3:170331/149‑1 (MQ=255)
                                                          aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                             >  1:46562/1‑149 (MQ=255)
                                                              cggctgtAATTTCTAAGATCAACCCGATCGAAGATACGCCTTCCGATGAAAACGGTACTCCGGTACCCATCGTACTTCACCCGCTGGGCGACCCGTCTCGTATGAACATCGGTCAGATCCTCGACCCCCACCTGGGTATGGCTGCGata                                                                         <  6:315870/145‑3 (MQ=255)
                                                              aGGGTGTAATTTCTAAGATCAACCCGATCGTAGATATGCCTTACGATGAAAACTGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACTTCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  7:216851/149‑1 (MQ=255)
                                                              aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  2:46562/149‑1 (MQ=255)
                                                               gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                        <  6:83318/149‑1 (MQ=255)
                                                                    tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                                   <  3:20680/149‑1 (MQ=255)
                                                                     aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                  >  3:282807/1‑149 (MQ=255)
                                                                      aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  5:402613/149‑1 (MQ=255)
                                                                          cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACACACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              >  5:212853/1‑148 (MQ=255)
                                                                              gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                         <  4:140057/149‑1 (MQ=255)
                                                                               aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                             >  2:335242/1‑144 (MQ=255)
                                                                                tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                        >  4:219184/1‑148 (MQ=255)
                                                                                   aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                    >  6:83901/1‑149 (MQ=255)
                                                                                   aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                    >  5:96886/1‑149 (MQ=255)
                                                                                     ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                  >  7:196736/1‑149 (MQ=255)
                                                                                            aaGATATGCCTTACGTTGAAAACTGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAACCCCCTGGGTATGGCTGCGAAAGGTATCGGAGACAAGATCAAAGCCATGCt                                            >  2:336521/1‑148 (MQ=255)
                                                                                            aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                            >  6:352079/1‑148 (MQ=255)
                                                                                            aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                            >  3:370784/1‑148 (MQ=255)
                                                                                             aGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                           >  8:423412/1‑132 (MQ=255)
                                                                                             aGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                           <  7:423412/132‑1 (MQ=255)
                                                                                                tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                       <  8:196736/149‑1 (MQ=255)
                                                                                                 atGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                        >  7:404171/1‑147 (MQ=255)
                                                                                                  tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                     <  3:219184/149‑1 (MQ=255)
                                                                                                    ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                   <  6:376232/149‑1 (MQ=255)
                                                                                                     cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                  >  6:50510/1‑149 (MQ=255)
                                                                                                         cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                               <  5:83901/148‑1 (MQ=255)
                                                                                                          gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                             <  1:97628/149‑1 (MQ=255)
                                                                                                          gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                              >  1:157813/1‑148 (MQ=255)
                                                                                                          gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGCTCAGATCCTCGAAACCCACCTCGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                              >  2:125321/1‑148 (MQ=255)
                                                                                                           ttGATAACGGTACGGCGGTAGACATCGTACTGAACCTGCTGGGCGTACCGTCTAGTATGAACATCGGTCAGATCCTCGAAACCCTCCTGGGTATGGCTGCGAATGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGgaagaa                            <  8:404171/148‑1 (MQ=255)
                                                                                                            tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                        <  3:140087/56‑1 (MQ=255)
                                                                                                            tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                        >  4:140087/1‑56 (MQ=255)
                                                                                                               aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                        >  7:290295/1‑149 (MQ=255)
                                                                                                                   ggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                     >  1:84749/1‑148 (MQ=255)
                                                                                                                    gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                   <  5:352079/149‑1 (MQ=255)
                                                                                                                    gTACGCCGGTAGACATAGTACTGAACCAGCTTGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                   <  1:336521/149‑1 (MQ=255)
                                                                                                                         cctgTAGACATCGTACTGAACCTGCTGGGCGTACCTTATCGTATGAACATCGGTCAGATCCTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAATCAGCAGCAAGAAGTCGCGAAACTgcg              <  3:325960/146‑1 (MQ=255)
                                                                                                                         ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg              >  8:47702/1‑149 (MQ=255)
                                                                                                                          cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc             <  5:50510/149‑1 (MQ=255)
                                                                                                                          cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc             <  4:282807/149‑1 (MQ=255)
                                                                                                                                  aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt     <  1:335242/149‑1 (MQ=255)
                                                                                                                                   tCGTACTGACCCAGCGGGGCGGACCGTCTCGTCTGACCATCGGTCAGATCCGCGAATCCCAGCTGGTTATGGCGGCGAAAGGTCGCGGCCACTAGCGGAACGCCATGCGTAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc    <  8:310996/149‑1 (MQ=255)
                                                                                                                                   tCGTACTGAACCCGCTGGGCGAACCGTCTCGTATGAACCTCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGGCGCGAAACTGCGCGAATTCATc    <  2:412258/149‑1 (MQ=255)
                                                                                                                                     gTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCa  <  8:290295/149‑1 (MQ=255)
                                                                                                                                         tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                          >  5:406204/1‑89 (MQ=255)
                                                                                                                                         tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                          <  6:406204/89‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
ATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184400‑4184681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTAC  >  NC_000913/4184394‑4184691
                                                                                                                                                     |                                                                                                                                                    
gtggctcttccgatctGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAG                                                                                                                                                     <  NB501016:230:HFVJLAFXY:4:11609:21836:10989/135‑1 (MQ=60)
 AAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGC                                                                                                                                                    >  NB501016:230:HFVJLAFXY:2:21306:17101:3159/1‑151 (MQ=60)
     TATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATAAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCATCTGGG                                                                                                                                                >  NB501016:230:HFVJLAFXY:1:21310:13203:3610/1‑151 (MQ=60)
      ATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGG                                                                                                                                                <  NB501016:230:HFVJLAFXY:1:21101:7329:10643/150‑1 (MQ=60)
          CTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAC                                                                                                                                           >  NB501016:230:HFVJLAFXY:1:21301:7514:7992/1‑151 (MQ=60)
           TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACC                                                                                                                                          >  NB501016:230:HFVJLAFXY:2:11201:16165:5786/1‑151 (MQ=60)
              ggcTTAAACGCCTTATCAAGCCTGGTGACAAGATGGCAGGTCGTCACGGTATCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAACACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTC                                                                                                                                       <  NB501016:230:HFVJLAFXY:4:21502:16222:13338/148‑1 (MQ=60)
              CGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTC                                                                                                                                       <  NB501016:230:HFVJLAFXY:4:21502:16213:13332/151‑1 (MQ=60)
               GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCT                                                                                                                                      <  NB501016:230:HFVJLAFXY:4:21412:26478:14851/151‑1 (MQ=60)
                  TACACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGCGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCG                                                                                                                                    <  NB501016:230:HFVJLAFXY:3:11409:6893:19935/150‑1 (MQ=60)
                     ACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATG                                                                                                                                <  NB501016:230:HFVJLAFXY:3:21404:21125:15883/151‑1 (MQ=60)
                               CAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAAATCGGTC                                                                                                                      <  NB501016:230:HFVJLAFXY:4:11404:14459:13370/151‑1 (MQ=60)
                               CAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTA                                                                                                                                  <  NB501016:230:HFVJLAFXY:1:21111:24510:10843/139‑1 (MQ=60)
                                  CCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGA                                                                                                                   <  NB501016:230:HFVJLAFXY:3:11601:9497:8436/151‑1 (MQ=60)
                                    TGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTC                                                                                                                      >  NB501016:230:HFVJLAFXY:1:21210:21039:8460/1‑146 (MQ=60)
                                            AGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAC                                                                                                         >  NB501016:230:HFVJLAFXY:1:21104:17608:16950/1‑151 (MQ=60)
                                            AGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAC                                                                                                         <  NB501016:230:HFVJLAFXY:4:21410:19352:9109/151‑1 (MQ=60)
                                            AGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGG                                                                                                                                                >  NB501016:230:HFVJLAFXY:4:11407:12022:18164/1‑112 (MQ=60)
                                                 GCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACC                                                                                                    <  NB501016:230:HFVJLAFXY:1:11306:11808:13051/151‑1 (MQ=60)
                                                      GCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGT                                                                                               <  NB501016:230:HFVJLAFXY:3:21510:22096:15980/151‑1 (MQ=60)
                                                          CACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATG                                                                                            <  NB501016:230:HFVJLAFXY:2:21302:5716:5153/150‑1 (MQ=60)
                                                              GTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGC                                                                                       <  NB501016:230:HFVJLAFXY:2:11307:15163:15279/151‑1 (MQ=60)
                                                               TAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCG                                                                                      >  NB501016:230:HFVJLAFXY:1:11108:18652:5344/1‑151 (MQ=60)
                                                                   AAGGGTGTAATTTCTAAGATCAACCCGATCGTAGATATGCCTTACGATGAAAACTGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACTTCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAG                                                                                  <  NB501016:230:HFVJLAFXY:4:11611:3721:11518/151‑1 (MQ=60)
                                                                         GTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCG                                                                            <  NB501016:230:HFVJLAFXY:2:11104:17395:14836/151‑1 (MQ=60)
                                                                          TAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGG                                                                           >  NB501016:230:HFVJLAFXY:2:21204:17452:1968/1‑151 (MQ=60)
                                                                           AATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGC                                                                          <  NB501016:230:HFVJLAFXY:3:21610:19859:11821/151‑1 (MQ=60)
                                                                               TCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAACACACCTGGGTATGGCTGCGAAAGGTATCGGCGAC                                                                       >  NB501016:230:HFVJLAFXY:3:21401:13112:5523/1‑150 (MQ=60)
                                                                                        AACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACG                                                             >  NB501016:230:HFVJLAFXY:3:11505:16494:14327/1‑151 (MQ=60)
                                                                                          CCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCC                                                           >  NB501016:230:HFVJLAFXY:4:11608:11560:4211/1‑151 (MQ=60)
                                                                                                 GAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTG                                                     >  NB501016:230:HFVJLAFXY:2:21307:16591:16252/1‑150 (MQ=60)
                                                                                                  AAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAG                                                                    <  NB501016:230:HFVJLAFXY:4:21609:2211:15378/134‑1 (MQ=60)
                                                                                                      TATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAC                                                 >  NB501016:230:HFVJLAFXY:4:21606:19315:9800/1‑149 (MQ=60)
                                                                                                       ATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGC                                              <  NB501016:230:HFVJLAFXY:2:21104:16413:4039/151‑1 (MQ=60)
                                                                                                              ACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAG                                        <  NB501016:230:HFVJLAFXY:3:11503:7301:9420/150‑1 (MQ=60)
                                                                                                               CGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAA                                      <  NB501016:230:HFVJLAFXY:1:11204:9557:15266/151‑1 (MQ=60)
                                                                                                               CGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGA                                       >  NB501016:230:HFVJLAFXY:3:11602:13023:7803/1‑150 (MQ=60)
                                                                                                               CGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGA                                       >  NB501016:230:HFVJLAFXY:1:11303:2166:5890/1‑150 (MQ=60)
                                                                                                                 ATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAT                                                                                                                                 <  NB501016:230:HFVJLAFXY:2:11302:14081:12031/58‑1 (MQ=60)
                                                                                                                    AAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGC                                 >  NB501016:230:HFVJLAFXY:4:21411:16639:1667/1‑151 (MQ=60)
                                                                                                                       ACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGT                                                                                                                                        >  NB501016:230:HFVJLAFXY:3:11602:20272:10748/1‑45 (MQ=60)
                                                                                                                        CGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAA                              >  NB501016:230:HFVJLAFXY:1:11202:3358:10860/1‑150 (MQ=60)
                                                                                                                         GGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAC                            <  NB501016:230:HFVJLAFXY:3:21602:25235:5475/151‑1 (MQ=60)
                                                                                                                         GGTACGCCGGTAGACATAGTACTGAACCAGCTTGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAC                            <  NB501016:230:HFVJLAFXY:1:21208:25420:20359/151‑1 (MQ=60)
                                                                                                                              GCCTGTAGACATCGTACTGAACCTGCTGGGCGTACCTTATCGTATGAACATCGGTCAGATCCTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAATCAGCAGCAAGAAGTCGCGAAACTGCGC                       <  NB501016:230:HFVJLAFXY:2:21212:16919:1349/151‑1 (MQ=60)
                                                                                                                               CCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCG                      <  NB501016:230:HFVJLAFXY:3:11409:19108:15874/151‑1 (MQ=60)
                                                                                                                                       CATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC              <  NB501016:230:HFVJLAFXY:1:21208:24687:15851/151‑1 (MQ=60)
                                                                                                                                              CTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGA                                                                   >  NB501016:230:HFVJLAFXY:3:21611:26021:4017/1‑91 (MQ=60)
                                                                                                                                                  ACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGCTCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCA                <  NB501016:230:HFVJLAFXY:3:21604:16420:3497/138‑1 (MQ=60)
                                                                                                                                                  ACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGT    >  NB501016:230:HFVJLAFXY:1:11210:24495:8869/1‑150 (MQ=60)
                                                                                                                                                   CCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTAC  >  NB501016:230:HFVJLAFXY:1:11205:15090:16309/1‑151 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
TAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTAC  >  NC_000913/4184394‑4184691

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: