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
A3 F68 I1 R1 24 78.6 3085910 98.1% 3027277 142.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C96.7% 390.4 / 5.7 121intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (62/55);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G96.7% 388.7 / 5.6 121intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (62/55);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913/4296294‑4296519
                                                                                       ||                                                                                                                                              
gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                                 >  8:50863/1‑149 (MQ=0)
 cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                <  3:157092/149‑1 (MQ=2)
 cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                >  1:197614/1‑149 (MQ=2)
  gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               >  4:79794/1‑149 (MQ=9)
  gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               >  3:17626/1‑149 (MQ=9)
  gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               >  2:235365/1‑149 (MQ=9)
  gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               >  5:191812/1‑149 (MQ=9)
  gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               >  7:244488/1‑149 (MQ=9)
     gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                            <  3:91733/149‑1 (MQ=255)
         aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                         >  6:23038/1‑148 (MQ=255)
         aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                         >  3:267588/1‑148 (MQ=255)
         aaCGCGTCTTATCAGGCCTAAGACAGACAGCGCAATAGCCTGATTAAGCGTGATTTTGAAGGTCGGATAAGGCGTTAA‑‑CGCGGCAACCGAAATCAATGCATTATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCTATGGtt                                                                         >  4:88202/1‑148 (MQ=11)
          aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                       >  8:131645/1‑149 (MQ=255)
           cgcgTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                      <  5:257959/149‑1 (MQ=11)
             cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                    <  5:46905/149‑1 (MQ=255)
               tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                  <  2:197614/149‑1 (MQ=255)
                cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                  >  3:16922/1‑148 (MQ=255)
                 ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                 >  7:149681/1‑148 (MQ=255)
                 ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                 >  1:235985/1‑148 (MQ=14)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                  <  6:150156/145‑1 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              >  3:216199/1‑149 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              >  6:247699/1‑149 (MQ=255)
                    tCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                             >  2:238954/1‑149 (MQ=255)
                      aGGGCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCCGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTATCTTAACCGTTGGTTAATTTTTCAATct                                                            <  7:156469/148‑1 (MQ=14)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                           <  7:348655/149‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                          >  8:258903/1‑134 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                          <  7:258903/134‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                           <  7:218511/149‑1 (MQ=255)
                       ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                          <  7:131645/149‑1 (MQ=255)
                        tcatACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCTTTTACGC‑‑CGCATCCGTCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  4:355031/146‑1 (MQ=17)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                          <  1:174173/148‑1 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                <  5:178669/142‑1 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                >  6:178669/1‑142 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  5:23038/149‑1 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     >  1:301385/1‑137 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     <  2:301385/137‑1 (MQ=255)
                          cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         >  4:149105/1‑147 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                              >  6:270090/1‑125 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                              <  5:270090/125‑1 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                       >  3:95121/1‑148 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                       >  1:116898/1‑148 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                       >  3:143282/1‑148 (MQ=255)
                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  6:98525/144‑1 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                   <  3:255600/149‑1 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     >  5:298809/1‑131 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                   >  3:294415/1‑149 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                   >  5:22973/1‑149 (MQ=255)
                                caaaCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 >  5:317909/4‑149 (MQ=255)
                                cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 >  5:222369/1‑149 (MQ=255)
                                cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 <  5:273421/149‑1 (MQ=255)
                                cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 >  7:196598/1‑149 (MQ=255)
                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                <  6:70483/149‑1 (MQ=255)
                                  gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     >  5:38674/1‑127 (MQ=255)
                                  gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     <  6:38674/127‑1 (MQ=255)
                                  gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                >  3:183362/1‑148 (MQ=255)
                                   aCAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              <  1:126874/149‑1 (MQ=255)
                                   aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              <  2:66471/149‑1 (MQ=255)
                                     aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                            <  4:294415/149‑1 (MQ=255)
                                     aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCATAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                             >  5:380034/1‑148 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  7:50863/149‑1 (MQ=255)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              >  2:161938/1‑142 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                    <  4:267588/148‑1 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                    <  4:16922/148‑1 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGATATCAAACCCtt                                    <  2:235985/148‑1 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCATAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                    <  6:380034/148‑1 (MQ=255)
                                               ccTGCTTTAGCGTGTTGTTGTCGGTCGGCTAAGGCGTTTACGC‑‑CGCAGCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  <  6:31730/149‑1 (MQ=21)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  <  6:269449/149‑1 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   >  1:127612/1‑148 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   >  3:277484/1‑148 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   >  6:220930/1‑148 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  <  5:123475/149‑1 (MQ=255)
                                                cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                 <  6:22973/149‑1 (MQ=255)
                                                cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                 >  1:364486/1‑149 (MQ=255)
                                                cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                 >  5:343421/1‑149 (MQ=255)
                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   <  8:34030/145‑1 (MQ=255)
                                                   aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  4:234993/148‑1 (MQ=255)
                                                    tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                             <  3:79794/149‑1 (MQ=255)
                                                       aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                          >  7:23097/1‑149 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       >  3:234012/1‑149 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       >  8:118488/1‑149 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  2:127048/140‑1 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               >  1:127048/1‑140 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACaagatcgg                      >  3:77964/1‑142 (MQ=21)
                                                             aTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                     <  5:103647/148‑1 (MQ=255)
                                                              ttttGTAGGTCGGTTAAGGCGTTTA‑‑CGACGCATCAGAGATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                   <  1:167147/149‑1 (MQ=25)
                                                              ttttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                    >  2:238789/1‑148 (MQ=255)
                                                              ttttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                    >  4:166215/1‑148 (MQ=255)
                                                              ttttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                   <  4:95121/149‑1 (MQ=255)
                                                                ttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                          >  4:196595/1‑92 (MQ=14)
                                                                ttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCACACTGAGCACACCCGTTATCACACCATTCCTTACAGCTCACGTTcct                 >  6:208420/1‑147 (MQ=17)
                                                                   tAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa              <  7:240081/149‑1 (MQ=255)
                                                                     ggTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt            >  1:394375/1‑149 (MQ=255)
                                                                       tCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg           >  4:322271/1‑148 (MQ=255)
                                                                         ggATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg        <  6:191812/149‑1 (MQ=255)
                                                                          gATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc       >  1:134295/1‑149 (MQ=255)
                                                                          gATAAGGCGTTTA‑‑CGACGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc       <  8:149681/149‑1 (MQ=255)
                                                                           aTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa         <  1:395505/146‑1 (MQ=255)
                                                                           aTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa         >  2:395505/1‑146 (MQ=255)
                                                                           aTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg      >  3:170362/1‑149 (MQ=255)
                                                                           aTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg      <  8:367525/149‑1 (MQ=255)
                                                                                       ||                                                                                                                                              
GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913/4296294‑4296519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ 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
ATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCT  >  NC_000913/4296292‑4296530
                                                                                           |                                                                                                                                                       
ATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAC                                                                                            <  NB501016:230:HFVJLAFXY:2:11210:16617:13541/151‑1 (MQ=1)
  GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCG                                                                                          <  NB501016:230:HFVJLAFXY:2:11308:22434:5506/151‑1 (MQ=5)
  GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCG                                                                                          >  NB501016:230:HFVJLAFXY:1:11312:17833:11784/1‑151 (MQ=5)
   CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGT                                                                                         >  NB501016:230:HFVJLAFXY:3:11612:18365:8019/1‑151 (MQ=7)
   CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGT                                                                                         >  NB501016:230:HFVJLAFXY:2:11104:16662:10320/1‑151 (MQ=7)
   CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGT                                                                                         >  NB501016:230:HFVJLAFXY:4:21408:7427:15960/1‑151 (MQ=7)
      CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGG                                                                                      <  NB501016:230:HFVJLAFXY:2:11207:22023:13493/151‑1 (MQ=16)
         TTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAC                                                                                            >  NB501016:230:HFVJLAFXY:4:11502:14233:12377/1‑142 (MQ=1)
          TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTA                                                                                   >  NB501016:230:HFVJLAFXY:2:21206:18966:4843/1‑150 (MQ=22)
            ACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATT                                                                                <  NB501016:230:HFVJLAFXY:3:21412:5811:15260/151‑1 (MQ=27)
              GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTT                                                                              <  NB501016:230:HFVJLAFXY:3:11409:9792:18505/151‑1 (MQ=33)
                 TCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTC                                                                            >  NB501016:230:HFVJLAFXY:2:11104:13210:7711/1‑150 (MQ=37)
                  CTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCA                                                                           >  NB501016:230:HFVJLAFXY:4:11603:14195:6216/1‑150 (MQ=39)
                  CTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCA                                                                           >  NB501016:230:HFVJLAFXY:1:21106:15242:20147/1‑150 (MQ=37)
                    TATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAC                                                                        >  NB501016:230:HFVJLAFXY:2:21107:22285:10487/1‑151 (MQ=45)
                    TATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTC                                                                            >  NB501016:230:HFVJLAFXY:3:11605:13071:1502/1‑147 (MQ=37)
                       gagggCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCCGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTATCTTAACCGTTGGTTAATTTTTCAATCTC                                                                      <  NB501016:230:HFVJLAFXY:4:11604:18449:11448/145‑1 (MQ=37)
                       CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT                                                                     <  NB501016:230:HFVJLAFXY:4:21604:20671:10621/151‑1 (MQ=50)
                       CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT                                                                     <  NB501016:230:HFVJLAFXY:4:21404:9819:4678/151‑1 (MQ=51)
                       CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT                                                                     >  NB501016:230:HFVJLAFXY:2:11106:12485:16229/1‑151 (MQ=51)
                       CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTT                                                                                    <  NB501016:230:HFVJLAFXY:4:21411:23493:4667/136‑1 (MQ=21)
                        aggctattgcgctgtctggcgtagGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTT                                                                    >  NB501016:230:HFVJLAFXY:2:21111:17842:8379/25‑151 (MQ=56)
                        AGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTT                                                                    <  NB501016:230:HFVJLAFXY:4:11511:9333:20097/151‑1 (MQ=53)
                         GGCCTTCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTT                                                                    <  NB501016:230:HFVJLAFXY:1:11308:10180:6902/150‑1 (MQ=50)
                         GGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTC                                                                   <  NB501016:230:HFVJLAFXY:3:11405:19917:9215/151‑1 (MQ=55)
                         GGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAA                                                                          <  NB501016:230:HFVJLAFXY:3:11610:16951:1120/144‑1 (MQ=41)
                         GGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTT                                                                               >  NB501016:230:HFVJLAFXY:1:21207:19431:1378/1‑139 (MQ=31)
                           CCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTC                                                                   <  NB501016:230:HFVJLAFXY:2:11306:15440:16080/149‑1 (MQ=55)
                            CTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAC                                                                 >  NB501016:230:HFVJLAFXY:1:11209:10751:13741/1‑150 (MQ=59)
                            CTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAC                                                                 >  NB501016:230:HFVJLAFXY:2:11305:13092:14623/1‑150 (MQ=59)
                            CTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAC                                                                 >  NB501016:230:HFVJLAFXY:2:11208:19990:7763/1‑150 (MQ=59)
                            CTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTT                                                                                        <  NB501016:230:HFVJLAFXY:3:21503:22907:19369/127‑1 (MQ=10)
                              ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTC                                                                   >  NB501016:230:HFVJLAFXY:3:11507:22226:12330/1‑146 (MQ=56)
                               CGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC                                                             <  NB501016:230:HFVJLAFXY:2:21203:1861:14185/151‑1 (MQ=60)
                               CGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC                                                             >  NB501016:230:HFVJLAFXY:2:21211:26380:13521/1‑151 (MQ=60)
                               CGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC                                                             >  NB501016:230:HFVJLAFXY:3:11405:17985:9022/1‑151 (MQ=60)
                               CGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTT                                                                               >  NB501016:230:HFVJLAFXY:3:21509:19100:8287/1‑133 (MQ=32)
                                 CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGG                                                           >  NB501016:230:HFVJLAFXY:4:11612:23150:3195/1‑151 (MQ=60)
                                 CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGG                                                           <  NB501016:230:HFVJLAFXY:3:21504:8054:12169/151‑1 (MQ=60)
                                 CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGG                                                           >  NB501016:230:HFVJLAFXY:3:21406:1651:4494/1‑151 (MQ=60)
                                 CCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGG                                                           >  NB501016:230:HFVJLAFXY:3:21601:14344:1875/1‑151 (MQ=60)
                                   AGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGT                                                          >  NB501016:230:HFVJLAFXY:2:11312:25144:16897/1‑150 (MQ=60)
                                   AGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTT                                                                               >  NB501016:230:HFVJLAFXY:3:11408:10353:8430/1‑129 (MQ=32)
                                    GACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT                                                        <  NB501016:230:HFVJLAFXY:1:11211:21764:12110/151‑1 (MQ=60)
                                      CAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCATAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATC                                                       >  NB501016:230:HFVJLAFXY:3:21611:14832:18539/1‑150 (MQ=60)
                                         CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC                                                   <  NB501016:230:HFVJLAFXY:4:11409:20815:15506/151‑1 (MQ=60)
                                           CAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC                                                   <  NB501016:230:HFVJLAFXY:3:21411:5392:13213/149‑1 (MQ=60)
                                           CAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT                                                        <  NB501016:230:HFVJLAFXY:1:11306:3178:2401/144‑1 (MQ=60)
                                               AGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTC                                              <  NB501016:230:HFVJLAFXY:1:21302:25006:15845/150‑1 (MQ=60)
                                               AGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTT                                                                              >  NB501016:230:HFVJLAFXY:4:11612:14541:6224/1‑118 (MQ=36)
                                                GCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCC                                             >  NB501016:230:HFVJLAFXY:2:21208:4191:5207/1‑150 (MQ=60)
                                                GCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCC                                             >  NB501016:230:HFVJLAFXY:1:11211:5549:14614/1‑150 (MQ=60)
                                                GCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCA                                            <  NB501016:230:HFVJLAFXY:3:11512:19686:8788/151‑1 (MQ=60)
                                                 CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAT                                           >  NB501016:230:HFVJLAFXY:1:21306:8302:18764/1‑151 (MQ=60)
                                                 CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAT                                           >  NB501016:230:HFVJLAFXY:3:21605:25100:9637/1‑151 (MQ=60)
                                                   TGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCC                                             >  NB501016:230:HFVJLAFXY:4:11406:25352:19732/1‑147 (MQ=60)
                                                     ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAG                                       <  NB501016:230:HFVJLAFXY:2:11205:13207:5620/151‑1 (MQ=60)
                                                        TAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTC                                    >  NB501016:230:HFVJLAFXY:4:11405:4289:3200/1‑151 (MQ=60)
                                                           CGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGG                                 >  NB501016:230:HFVJLAFXY:2:21111:8196:4639/1‑151 (MQ=60)
                                                            GTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAagatcggg                                >  NB501016:230:HFVJLAFXY:2:11204:3301:16651/1‑143 (MQ=60)
                                                            GTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAC                                         >  NB501016:230:HFVJLAFXY:1:11211:13507:12683/1‑142 (MQ=60)
                                                              GATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTC                              >  NB501016:230:HFVJLAFXY:3:21511:7988:15451/1‑151 (MQ=60)
                                                              GATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTT                               <  NB501016:230:HFVJLAFXY:3:11508:4741:11660/150‑1 (MQ=60)
                                                               ATTTTGTAGGTCGGTTAAGGCGTTTA‑‑CGACGCATCAGAGATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA                             <  NB501016:230:HFVJLAFXY:1:11307:8273:1151/151‑1 (MQ=60)
                                                                 TTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTT                                                                                    <  NB501016:230:HFVJLAFXY:2:21103:8438:10628/94‑1 (MQ=24)
                                                                    GTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAG                        <  NB501016:230:HFVJLAFXY:4:21408:8372:1400/151‑1 (MQ=60)
                                                                      AGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTG                      >  NB501016:230:HFVJLAFXY:1:21312:3261:7778/1‑151 (MQ=60)
                                                                           GGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG                 >  NB501016:230:HFVJLAFXY:1:11212:21373:19462/1‑151 (MQ=60)
                                                                            GATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGC                >  NB501016:230:HFVJLAFXY:2:11310:21307:12596/1‑151 (MQ=60)
                                                                            GATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAG                   <  NB501016:230:HFVJLAFXY:1:21312:20889:11483/148‑1 (MQ=60)
                                                                                 GGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTT           >  NB501016:230:HFVJLAFXY:4:11511:25436:15300/1‑151 (MQ=60)
                                                                                 GGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTT           <  NB501016:230:HFVJLAFXY:1:11312:17921:8345/151‑1 (MQ=60)
                                                                                   CGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC         >  NB501016:230:HFVJLAFXY:1:11201:22611:2003/1‑151 (MQ=60)
                                                                                    GTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCT        >  NB501016:230:HFVJLAFXY:3:21503:12501:18040/1‑151 (MQ=60)
                                                                                      ttacgC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTT       <  NB501016:230:HFVJLAFXY:2:21110:2350:8726/145‑1 (MQ=60)
                                                                                       tacgC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTG      <  NB501016:230:HFVJLAFXY:1:21106:7579:17672/146‑1 (MQ=60)
                                                                                        G‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCT  <  NB501016:230:HFVJLAFXY:4:11612:24416:17490/151‑1 (MQ=60)
                                                                                        G‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCT  <  NB501016:230:HFVJLAFXY:3:21604:5221:6966/151‑1 (MQ=60)
                                                                                         cgC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC               >  NB501016:230:HFVJLAFXY:3:21409:13719:2559/3‑139 (MQ=60)
                                                                                           |                                                                                                                                                       
ATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCT  >  NC_000913/4296292‑4296530

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: