Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,505 G→A 19.0% intergenic (+101/‑81) 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,5050GA19.0% 263.8 / 47.2 115intergenic (+101/‑81)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (57/36);  new base A (9/13);  total (66/49)
Fisher's exact test for biased strand distribution p-value = 9.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 115 total aligned reads displayed.

AACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTG  >  NC_000913/1530407‑1530603
                                                                                                  |                                                                                                  
aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGCGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                                   >  1:2241047/1‑100 (MQ=255)
 acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAt                                                                                                  >  1:411206/1‑100 (MQ=255)
  caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                                 >  2:2152700/1‑100 (MQ=255)
   aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATTTACATCACAAttt                                                                                                <  2:1555466/100‑1 (MQ=255)
    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)
        cAGAAAGAGAATAATATAAAAAACAGACAGAATGAGAAGAGAGCAAGGGACGCGGAAGTGCAAGATCATGAAAATCGTTAGGTACATCACAATTATCagg                                                                                           <  1:357700/100‑3 (MQ=255)
        cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                           <  2:886066/100‑2 (MQ=37)
        cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                           <  2:1264063/100‑2 (MQ=37)
        cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                           <  1:1403369/100‑2 (MQ=38)
         aGAAAGAGAAAAATATAGAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                           <  1:2510942/99‑2 (MQ=37)
         aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                            <  1:1544034/98‑1 (MQ=38)
               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)
                         aaaaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAc                                                                          >  2:45930/1‑100 (MQ=255)
                           aaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                        >  2:347095/1‑100 (MQ=255)
                            atacatacAGAAAGACATGAGAGCAAGGGACGTGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                       <  1:1141147/100‑1 (MQ=255)
                             tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                      <  2:769139/100‑1 (MQ=255)
                             tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                      >  1:1130051/1‑100 (MQ=255)
                             tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                      >  2:177477/1‑100 (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                                                              >  1:120259/1‑100 (MQ=255)
                                     gAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaaa                                                              >  2:680044/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)
                                      aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                                             <  2:1549279/100‑1 (MQ=255)
                                       aataCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                            <  2:2542145/97‑1 (MQ=255)
                                       aaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                            >  1:230277/1‑100 (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                                                      >  2:1061034/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)
                                                        gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                                           <  2:483740/100‑1 (MQ=255)
                                                         aCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAGCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                                           <  1:2516252/99‑1 (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)
                                                               cagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                                    >  2:292062/1‑100 (MQ=255)
                                                               cagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                                    >  1:719673/1‑100 (MQ=255)
                                                                agggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGAt                                   <  2:2547138/100‑1 (MQ=255)
                                                                 gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                                  >  2:1603313/1‑100 (MQ=255)
                                                                 gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                                  >  2:733275/1‑100 (MQ=255)
                                                                 gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                                  >  2:2286181/1‑100 (MQ=255)
                                                                   gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                                >  2:887803/1‑100 (MQ=255)
                                                                   gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                                >  2:1821745/1‑100 (MQ=255)
                                                                   gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                                >  2:630182/1‑100 (MQ=255)
                                                                   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)
                                                                        aagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAgg                           >  1:740780/4‑100 (MQ=37)
                                                                          gaTGAAAATCGTGATTTACATCGCAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                         >  1:1501379/2‑100 (MQ=37)
                                                                          cATGAAAATTGTGATGTAAATCACGATTTTCATATTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCGTTGGAATAACTACAAAGATAAAAAAg                            <  1:1765937/97‑1 (MQ=255)
                                                                            tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTAAAAAGATAAAAAAGGACGc                       >  1:1046412/1‑100 (MQ=255)
                                                                            tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                       >  2:1871875/1‑100 (MQ=255)
                                                                               aaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACGAAAAAGATAAAAAAGGACGCCAg                    <  2:2195246/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)
                                                                                           aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa        <  2:2476283/100‑1 (MQ=255)
                                                                                           aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAAGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa        <  1:2256852/100‑1 (MQ=255)
                                                                                            aaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaaa       >  2:2499482/1‑100 (MQ=255)
                                                                                             aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAg      <  1:710130/100‑1 (MQ=255)
                                                                                                 cGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTg  <  2:1941053/100‑1 (MQ=255)
                                                                                                 cGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTg  >  1:1869870/1‑100 (MQ=255)
                                                                                                 cGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTg  >  1:592459/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTG  >  NC_000913/1530407‑1530603

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: