Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,499 A→C 20.4% intergenic (+95/‑87) 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,4990AC20.4% 242.2 / 61.8 108intergenic (+95/‑87)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base A (54/32);  new base C (9/13);  total (63/45)
Fisher's exact test for biased strand distribution p-value = 8.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 113 total aligned reads displayed.

ATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAA  >  NC_000913/1530401‑1530598
                                                                                                  |                                                                                                   
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)
         aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATTTACATCACAAttt                                                                                           <  2:1555466/100‑1 (MQ=255)
         aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAttt                                                                                           <  2:1043057/100‑1 (MQ=38)
          acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGCCATGAGAGAAAGGAACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAGCACGAttt                                                                                           >  1:511337/1‑99 (MQ=255)
          acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                          >  2:737427/1‑100 (MQ=255)
            acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                                        >  2:1107958/1‑100 (MQ=255)
            acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                        <  2:2544276/100‑1 (MQ=38)
            acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                        <  2:2117673/100‑1 (MQ=38)
              cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  1:1403369/100‑2 (MQ=38)
              cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  2:1264063/100‑2 (MQ=37)
              cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  2:886066/100‑2 (MQ=37)
              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:411208/100‑1 (MQ=255)
                          aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                          <  2:2270204/100‑1 (MQ=255)
                            tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                                        <  2:828171/100‑1 (MQ=255)
                             ataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCt                                                                       >  2:64122/1‑100 (MQ=255)
                             ataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCt                                                                       >  1:663024/1‑100 (MQ=255)
                                 aaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                   >  2:347095/1‑100 (MQ=255)
                                  atacatacAGAAAGACATGAGAGCAAGGGACGTGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                  <  1:1141147/100‑1 (MQ=255)
                                   tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                 >  2:177477/1‑100 (MQ=255)
                                   tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                 >  1:1130051/1‑100 (MQ=255)
                                     catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                               >  2:77831/1‑100 (MQ=255)
                                       tacaGAAAGACATGAGAGCAAGGGACGGGGCATGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                             >  1:1025934/1‑100 (MQ=255)
                                       tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                             >  2:1610668/1‑100 (MQ=255)
                                       tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                             >  2:339871/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)
                                           gAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                          >  1:316098/1‑99 (MQ=255)
                                            aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                                        <  2:1925482/100‑1 (MQ=255)
                                            aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                                        <  1:1526239/100‑1 (MQ=255)
                                             aataCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                       <  2:2542145/97‑1 (MQ=255)
                                             aaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                       <  2:1811917/100‑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                                                   >  1:774336/1‑100 (MQ=255)
                                                 cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                   >  2:882423/1‑100 (MQ=255)
                                                   tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                 >  1:1036017/1‑100 (MQ=255)
                                                   tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                 >  1:744695/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)
                                                        gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                            <  1:915475/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)
                                                                   ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                 >  2:2328586/1‑100 (MQ=255)
                                                                   ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                 >  2:1019932/1‑100 (MQ=255)
                                                                    gcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                                <  2:1077245/100‑1 (MQ=255)
                                                                     cagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                               >  1:719673/1‑100 (MQ=255)
                                                                     cagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                               >  2:292062/1‑100 (MQ=255)
                                                                      agggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGAt                              <  2:2547138/100‑1 (MQ=255)
                                                                       gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                             >  2:2286181/1‑100 (MQ=255)
                                                                       gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                             >  2:733275/1‑100 (MQ=255)
                                                                       gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                             >  2:1603313/1‑100 (MQ=255)
                                                                         gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                           >  2:630182/1‑100 (MQ=255)
                                                                         gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                           >  2:887803/1‑100 (MQ=255)
                                                                         gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                           >  2:1821745/1‑100 (MQ=255)
                                                                         aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                           >  1:2020236/9‑100 (MQ=21)
                                                                         aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                           >  1:2005666/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)
                                                                            naaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  2:2009310/6‑100 (MQ=25)
                                                                            caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  2:1034809/6‑100 (MQ=25)
                                                                            aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  2:851253/1‑100 (MQ=255)
                                                                             gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                       <  1:1211247/100‑1 (MQ=255)
                                                                             aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                       >  2:120691/5‑100 (MQ=25)
                                                                                gaTGAAAATCGTGATTTACATCGCAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                    >  1:1501379/2‑100 (MQ=37)
                                                                                cATGAAAATTGTGATGTAAATCACGATTTTCATATTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCGTTGGAATAACTACAAAGATAAAAAAg                       <  1:1765937/97‑1 (MQ=255)
                                                                                  tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTAAAAAGATAAAAAAGGACGc                  >  1:1046412/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)
                                                                                          gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg          <  2:1492379/100‑1 (MQ=255)
                                                                                               gTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTc     <  1:2286593/100‑1 (MQ=255)
                                                                                                tAAATCACGATTTTCATCTTTGCTTTAACGCCTAAAGGTGATCAAAATTTATCCTTTTGAATAACTAAAGAGATAAAAAAGGACGCCAGGTGAGTATTCa    >  1:386725/1‑100 (MQ=255)
                                                                                                 aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTc     >  2:2061303/1‑98 (MQ=255)
                                                                                                 aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAAGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa   <  1:2256852/100‑1 (MQ=255)
                                                                                                  aaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaaa  >  2:2499482/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAA  >  NC_000913/1530401‑1530598

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