Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,490 T→C 18.8% intergenic (+86/‑96) ydcD → / → yncI putative immunity protein for RhsE/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,4900TC18.8% 244.3 / 46.5 111intergenic (+86/‑96)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (55/35);  new base C (9/12);  total (64/47)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTG  >  NC_000913/1530392‑1530588
                                                                                                  |                                                                                                  
tAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtg                                                                                                   <  2:1696110/100‑1 (MQ=255)
 aaaTTGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGCGCAGGGCAAGATCATGAAAATtgt                                                                                                  <  1:1882385/100‑1 (MQ=255)
   aTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGa                                                                                                <  1:123489/100‑1 (MQ=255)
    tGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGAt                                                                                               <  1:1098622/100‑1 (MQ=255)
    tGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGAt                                                                                               <  1:775626/100‑1 (MQ=255)
     ggTCATAATGAACAACACCAGAAAGAGCAAAATATAAAATACATACAGAACGACATGAGAGCAAGGGACGGGGCAGTGCACGATCATGAAAATtgt                                                                                                  >  1:2445283/1‑96 (MQ=255)
      gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCACGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGt                                                                                             <  1:1378224/100‑1 (MQ=255)
       tCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTa                                                                                            <  2:1726606/100‑1 (MQ=255)
       tCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTa                                                                                            <  1:918370/100‑1 (MQ=255)
         aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaaa                                                                                          >  2:1925696/1‑100 (MQ=255)
            aTGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCa                                                                                       >  1:921078/1‑100 (MQ=255)
             tGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAc                                                                                      <  1:1888217/100‑1 (MQ=255)
               aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGCGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                    >  1:2241047/1‑100 (MQ=255)
                acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAt                                                                                   >  1:411206/1‑100 (MQ=255)
                 caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                  >  2:2294400/1‑100 (MQ=255)
                 caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                  >  2:2152700/1‑100 (MQ=255)
                  aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAttt                                                                                 <  2:1043057/100‑1 (MQ=38)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGCCATGAGAGAAAGGAACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAGCACGAttt                                                                                 >  1:511337/1‑99 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                >  2:737427/1‑100 (MQ=255)
                    cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGTCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                               <  1:1376217/100‑1 (MQ=255)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                              >  2:1107958/1‑100 (MQ=255)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                              <  2:2117673/100‑1 (MQ=38)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                              <  2:2544276/100‑1 (MQ=38)
                       cAGAAAGAGAATAATATAAAAAACAGACAGAATGAGAAGAGAGCAAGGGACGCGGAAGTGCAAGATCATGAAAATCGTTAGGTACATCACAATTATCagg                                                                            <  1:357700/100‑3 (MQ=255)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  2:1264063/100‑2 (MQ=37)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  2:886066/100‑2 (MQ=37)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  1:1403369/100‑2 (MQ=38)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACAGGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  2:2040570/100‑2 (MQ=37)
                        aGAAAGAGAAAAATATAGAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  1:2510942/99‑2 (MQ=37)
                        aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                             <  1:1544034/98‑1 (MQ=38)
                        aGAAAGAGAAAAATACAAAATACATACAGAAAGACATGAGAGCAAGGTACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                              <  1:1327594/97‑1 (MQ=37)
                              agaAAAATATGAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAATTCTTGATTTACATCACAATTTTCAtgatcttgc                                                                     <  1:1700109/100‑9 (MQ=21)
                               gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                    <  2:2121754/100‑1 (MQ=255)
                                 aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                  >  1:555110/1‑100 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                <  2:2270204/100‑1 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                <  2:411208/100‑1 (MQ=255)
                                     tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                              <  2:828171/100‑1 (MQ=255)
                                      ataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCt                                                             >  1:663024/1‑100 (MQ=255)
                                      ataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCt                                                             >  2:64122/1‑100 (MQ=255)
                                        aaaaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAc                                                           >  2:45930/1‑100 (MQ=255)
                                          aaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                         >  2:347095/1‑100 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGACGTGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                        <  1:1141147/100‑1 (MQ=255)
                                            tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                       >  1:1130051/1‑100 (MQ=255)
                                            tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                       >  2:177477/1‑100 (MQ=255)
                                            tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                       <  2:769139/100‑1 (MQ=255)
                                                tacaGAAAGACATGAGAGCAAGGGACGGGGCATGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                   >  1:1025934/1‑100 (MQ=255)
                                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                   >  2:339871/1‑100 (MQ=255)
                                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                   >  2:1610668/1‑100 (MQ=255)
                                                  caGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGGTCaa                                                 >  1:1625565/1‑100 (MQ=255)
                                                  caGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                 >  1:2547119/1‑100 (MQ=255)
                                                   aGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                >  2:1199960/1‑100 (MQ=255)
                                                    gAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaaa                                               >  2:680044/1‑100 (MQ=255)
                                                    gAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaaa                                               >  1:120259/1‑100 (MQ=255)
                                                     aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                              <  2:1925482/100‑1 (MQ=255)
                                                     aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                              <  2:1549279/100‑1 (MQ=255)
                                                      aataCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                             <  2:2542145/97‑1 (MQ=255)
                                                      aaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                             <  1:1532639/100‑1 (MQ=255)
                                                        gACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                           <  2:663027/100‑1 (MQ=255)
                                                         aCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTaa                                          <  2:2322558/100‑1 (MQ=255)
                                                         aCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                           <  2:1950970/99‑1 (MQ=255)
                                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                         >  2:882423/1‑100 (MQ=255)
                                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                         >  1:774336/1‑100 (MQ=255)
                                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                       >  1:1036017/1‑100 (MQ=255)
                                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                       >  1:744695/1‑100 (MQ=255)
                                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                       >  2:1061034/1‑100 (MQ=255)
                                                              agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCttt                                     >  2:2308674/1‑100 (MQ=255)
                                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAACTTTg                                    >  1:1984341/1‑100 (MQ=255)
                                                                agCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTgg                                   >  2:898993/1‑100 (MQ=255)
                                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGGCTACAGGTGATCAAAATTTAACCTTTGGa                                  >  2:1871559/1‑100 (MQ=255)
                                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                  <  2:157/100‑1 (MQ=255)
                                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCGTTGGTTTAACGCCTACAGGGGATCAAAATTTAAGCTTTgg                                   >  1:2575313/1‑99 (MQ=255)
                                                                    aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATa                               <  1:27358/100‑1 (MQ=255)
                                                                    aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATa                               <  1:737424/100‑1 (MQ=255)
                                                                       gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                            <  2:483740/100‑1 (MQ=255)
                                                                        aCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAGCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                            <  1:2516252/99‑1 (MQ=255)
                                                                           gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                        >  2:557595/1‑100 (MQ=255)
                                                                           gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                        >  1:1722882/1‑100 (MQ=255)
                                                                            ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                       >  2:2328586/1‑100 (MQ=255)
                                                                            ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                       >  2:1019932/1‑100 (MQ=255)
                                                                             gcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                      <  2:1077245/100‑1 (MQ=255)
                                                                               agggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGAt                    <  2:2547138/100‑1 (MQ=255)
                                                                                gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                   >  2:733275/1‑100 (MQ=255)
                                                                                gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                   >  2:1603313/1‑100 (MQ=255)
                                                                                  gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                 >  2:887803/1‑100 (MQ=255)
                                                                                  gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                 >  2:1821745/1‑100 (MQ=255)
                                                                                  aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                 >  1:2005666/9‑100 (MQ=21)
                                                                                  aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                 >  1:2020236/9‑100 (MQ=21)
                                                                                  aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                 >  1:238525/9‑100 (MQ=25)
                                                                                   caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                >  1:741596/1‑100 (MQ=255)
                                                                                    gcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaa               >  2:1328246/7‑100 (MQ=25)
                                                                                     caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa              >  2:1034809/6‑100 (MQ=25)
                                                                                     aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa              >  2:851253/1‑100 (MQ=255)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg             >  2:120691/5‑100 (MQ=25)
                                                                                       aagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAgg            >  1:740780/4‑100 (MQ=37)
                                                                                         gaTGAAAATCGTGATTTACATCGCAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc          >  1:1501379/2‑100 (MQ=37)
                                                                                           tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTAAAAAGATAAAAAAGGACGc        >  1:1046412/1‑100 (MQ=255)
                                                                                           tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc        >  2:1871875/1‑100 (MQ=255)
                                                                                              aaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACGAAAAAGATAAAAAAGGACGCCAg     <  2:2195246/100‑1 (MQ=255)
                                                                                                 ttGTGATGTAAATCGCGTTTTTCATCTTTGCTTTAACGGCTACAGGTGATCAAAAGTTAACCTTTGGAAGAACTAAAAAGATAAAAAAGGACGCCAGGTg  <  1:421201/100‑1 (MQ=255)
                                                                                                 ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAAAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg  <  1:1680558/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTG  >  NC_000913/1530392‑1530588

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