Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,898 A→G 10.0% M13T (ATG→ACG) ‡ yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8980AG10.0% 401.7 / 28.0 130M13T (ATG→ACG) ‡yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (54/63);  new base G (7/6);  total (61/69)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 132 total aligned reads displayed.

TTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATA  >  NC_000913/383752‑384046
                                                                                                                                                  |                                                                                                                                                    
ttACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                                                                    <  3:1152133/149‑1 (MQ=255)
ttACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                                                                    <  2:366145/149‑1 (MQ=255)
ttACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                                                                    <  1:512413/149‑1 (MQ=255)
  aCCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                                                                  <  1:336755/149‑1 (MQ=255)
   cccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCt                                                                                                                                                 <  1:1018228/149‑1 (MQ=255)
     cGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                                                                               <  1:1744936/149‑1 (MQ=255)
      ggCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                                                                              <  1:1512225/149‑1 (MQ=255)
       gCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                                                                             <  2:287934/149‑1 (MQ=255)
       gCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                                                                             <  2:1755927/149‑1 (MQ=255)
         aaaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                                                                            >  4:1027879/1‑148 (MQ=255)
         aaaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                                                                            <  3:1027879/148‑1 (MQ=255)
          aaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttt                                                                                                                                          <  3:358757/149‑4 (MQ=255)
              tAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                                                                      <  3:812184/149‑1 (MQ=255)
              tAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                                                                      <  4:741108/149‑1 (MQ=255)
              tAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                                                                      >  4:677285/1‑149 (MQ=255)
                aGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    <  3:677285/149‑1 (MQ=255)
                aGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    <  4:1451971/149‑1 (MQ=255)
                  ggTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                                                                  >  2:1074781/1‑149 (MQ=255)
                   gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATg                                                                                                                                              >  1:1748915/1‑136 (MQ=255)
                   gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    <  4:173223/146‑1 (MQ=255)
                   gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    >  3:173223/1‑146 (MQ=255)
                   gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTCGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaa                                                                                                                                 >  2:1010754/1‑149 (MQ=255)
                   gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTCGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaa                                                                                                                                 <  1:1010754/149‑1 (MQ=255)
                    ttACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                                                >  2:110586/1‑149 (MQ=255)
                      aCGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                                                                     <  3:935848/142‑1 (MQ=255)
                          aGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCa                                                                                                                          <  3:1557575/149‑1 (MQ=255)
                           gTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                         <  1:725343/149‑1 (MQ=255)
                           gTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                         <  1:1074781/149‑1 (MQ=255)
                               aTTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAGAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           >  1:553298/1‑143 (MQ=255)
                               aTTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           >  1:781935/1‑143 (MQ=255)
                               aTTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           <  2:781935/143‑1 (MQ=255)
                               aTTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           <  2:553298/143‑1 (MQ=255)
                                  aTCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                                                                  >  4:1167/1‑149 (MQ=255)
                                              cgcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                                                                      <  1:110586/149‑1 (MQ=255)
                                                 cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                   <  2:1323610/149‑1 (MQ=255)
                                                 cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                   >  2:1233185/1‑149 (MQ=255)
                                                   gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                                                >  3:1483488/1‑118 (MQ=255)
                                                   gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                                                <  4:1483488/118‑1 (MQ=255)
                                                     cacGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                                                                               <  3:1167/149‑1 (MQ=255)
                                                           ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTATAAACCATAGCAACATTGCTATGGACGTTCAGATAAATATCCAGGTTCTCATGCAATTTGCTCAAAGAATCAtttt                                                                                         >  3:1600924/1‑149 (MQ=255)
                                                           ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                   <  1:1050724/139‑1 (MQ=255)
                                                           ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                   >  2:1050724/1‑139 (MQ=255)
                                                           ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                         <  1:1233185/149‑1 (MQ=255)
                                                             ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                                                                               >  2:1624194/1‑141 (MQ=255)
                                                             ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                                                                               <  1:1624194/141‑1 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                                                                            >  1:1096912/1‑92 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                                                                   >  1:1301422/1‑149 (MQ=255)
                                                                   ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                                                                           >  3:463314/1‑123 (MQ=255)
                                                                   ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                                                                           <  4:463314/123‑1 (MQ=255)
                                                                       gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                                                                             >  4:806214/1‑149 (MQ=255)
                                                                         aGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAGACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    >  3:922105/1‑92 (MQ=255)
                                                                         aGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAGACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    <  4:922105/92‑1 (MQ=255)
                                                                         aGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    >  1:836109/1‑92 (MQ=255)
                                                                         aGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                    <  2:836109/92‑1 (MQ=255)
                                                                                 tGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTa                                                                                        <  2:1779138/128‑1 (MQ=255)
                                                                                 tGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTa                                                                                        >  1:1779138/1‑128 (MQ=255)
                                                                                 tGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                                                                   >  2:683314/1‑149 (MQ=255)
                                                                                   gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtt                                                                                           >  4:487556/1‑123 (MQ=255)
                                                                                     gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt                                                               >  4:1687882/1‑149 (MQ=255)
                                                                                           aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt                                                         >  1:1208945/1‑149 (MQ=255)
                                                                                                ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGATTCATTTTATGaa                                                                                    >  1:1230309/1‑117 (MQ=255)
                                                                                                ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGATTCATTTTATGaa                                                                                    <  2:1230309/117‑1 (MQ=255)
                                                                                                ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                                                    <  2:1208945/149‑1 (MQ=255)
                                                                                                 tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                                                   >  1:1324288/1‑149 (MQ=255)
                                                                                                  ttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                                                                      >  1:13014/1‑65 (MQ=255)
                                                                                                  ttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                                                                      <  2:13014/65‑1 (MQ=255)
                                                                                                  ttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                                                   >  1:280784/1‑148 (MQ=255)
                                                                                                     ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                                                                               >  4:480379/1‑69 (MQ=255)
                                                                                                     ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                                                                               <  3:480379/69‑1 (MQ=255)
                                                                                                          aCTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCcgc                                          >  4:697758/1‑149 (MQ=255)
                                                                                                           cTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                                                                          <  1:1575306/116‑1 (MQ=255)
                                                                                                           cTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                                                                          >  2:1575306/1‑116 (MQ=255)
                                                                                                               gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt                                                         >  3:241464/1‑129 (MQ=255)
                                                                                                               gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt                                                         <  4:241464/129‑1 (MQ=255)
                                                                                                                tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAtt                                    <  2:1301422/149‑1 (MQ=255)
                                                                                                                tttACGTCCTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCTGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTATAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAtt                                    <  4:1167050/149‑1 (MQ=255)
                                                                                                                      tCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATc                                                                                              >  1:525295/1‑85 (MQ=255)
                                                                                                                        ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTa                                                                                 <  2:1155927/96‑1 (MQ=255)
                                                                                                                         tAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCcgcg                                         <  2:845223/135‑1 (MQ=255)
                                                                                                                           gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt                                                        <  4:678937/118‑1 (MQ=255)
                                                                                                                           gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt                                                        >  3:678937/1‑118 (MQ=255)
                                                                                                                           gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCa                         >  1:989172/1‑149 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAgg                   >  2:1836928/1‑149 (MQ=255)
                                                                                                                                  tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCcgcg                                         >  3:250062/1‑126 (MQ=255)
                                                                                                                                  tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGc                  <  1:835433/149‑1 (MQ=255)
                                                                                                                                      gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGc                  <  1:1042467/143‑1 (MQ=255)
                                                                                                                                      aaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATtgt                                  <  2:509241/129‑1 (MQ=255)
                                                                                                                                       tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgtt                                                <  2:1678914/113‑1 (MQ=255)
                                                                                                                                        ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                                                   >  2:1051854/1‑110 (MQ=255)
                                                                                                                                         cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAc                                                     <  3:739890/107‑1 (MQ=255)
                                                                                                                                         cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATg           >  3:122316/1‑149 (MQ=255)
                                                                                                                                         cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                                                                          <  4:1254389/86‑1 (MQ=255)
                                                                                                                                         cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                                                                          >  3:1254389/1‑86 (MQ=255)
                                                                                                                                         cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATg           >  3:1415080/1‑149 (MQ=255)
                                                                                                                                           tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGa                                                      <  1:377867/104‑1 (MQ=255)
                                                                                                                                           tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGa                                                      >  2:377867/1‑104 (MQ=255)
                                                                                                                                              cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGa                <  3:554835/139‑1 (MQ=255)
                                                                                                                                              cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATc              >  1:888766/1‑141 (MQ=255)
                                                                                                                                              cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATc              <  2:888766/141‑1 (MQ=255)
                                                                                                                                                  aTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCAtata  >  2:1119602/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATA  >  NC_000913/383752‑384046

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