Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,069,219 G→A 6.0% D89N (GAC→AAC)  aldB‑II → aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,069,2190GA 6.0% 183.1 / 2.0 67D89N (GAC→AAC) aldB‑IIaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/36);  new base A (2/2);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GTCAATGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATCCGCGCCCAGGAGGGTGG  >  NC_002947/3069078‑3069367
                                                                                                                                             |                                                                                                                                                    
gTCAATGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGTTGATAATTGTCTTCCGCATc                                                                                                                                              <  2:1418815/150‑1 (MQ=255)
   aaTGGCCAGCCGATTGCCGAATTCCCGCGCTCCCCGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGaa                                                                                                                                           <  1:383142/150‑1 (MQ=255)
   aaTGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGaa                                                                                                                                           >  2:1509954/1‑150 (MQ=255)
    aTGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAAc                                                                                                                                          <  1:577783/150‑1 (MQ=255)
    aTGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAAc                                                                                                                                          >  2:111164/1‑150 (MQ=255)
     tGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGTCCGCATCGAACa                                                                                                                                         >  1:1030014/1‑150 (MQ=255)
      ggCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAg                                                                                                                                        <  1:986338/150‑1 (MQ=255)
      ggCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGACGAAGCCTGGGGCCAGACCTCGGCGCAAGACCGGGCGCGTGTACTGCTGAAAATTGCCGACCGCATCGAACAg                                                                                                                                        <  1:1571134/150‑1 (MQ=255)
         cAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAAGTGCCGACCGCATCGAACAGAAc                                                                                                                                     >  1:1560866/1‑150 (MQ=255)
          aGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAAcc                                                                                                                                    <  2:533109/150‑1 (MQ=255)
               aTTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGaa                                                                                                                               >  2:175613/1‑150 (MQ=255)
               aTTGCCGAATTCACGCGCTCCACGGCCCAGGACGGCGAGCGCGCCATGGACGCCGCGCATGCCGCCGCCGAAGCCTGGTGCAAGACGTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGaa                                                                                                                               <  2:1055230/150‑1 (MQ=255)
                   ccGAATTCCCGCGCTCCACGGCCCAGGAGGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGAGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGc                                                                                                                           >  2:999872/1‑150 (MQ=255)
                   ccGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGc                                                                                                                           <  2:1579527/150‑1 (MQ=255)
                    cGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCt                                                                                                                          <  2:491837/150‑1 (MQ=255)
                      aaTTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTgg                                                                                                                        <  1:1509954/150‑1 (MQ=255)
                           ccGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGtt                                                                                                                   >  1:606306/1‑150 (MQ=255)
                           ccGCGCTCCACGGCCAAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCATGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGtt                                                                                                                   >  2:204699/1‑150 (MQ=255)
                            cgcgCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTa                                                                                                                  >  2:1230339/1‑150 (MQ=255)
                             gcgcTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTAc                                                                                                                 <  2:6527/150‑1 (MQ=255)
                             gcgcTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTAc                                                                                                                 <  1:934433/150‑1 (MQ=255)
                              cgcTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTAcc                                                                                                                <  2:376890/150‑1 (MQ=255)
                                   cACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAg                                                                                                           <  1:335541/150‑1 (MQ=255)
                                    aCGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGc                                                                                                          >  1:333503/1‑150 (MQ=255)
                                     cGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCt                                                                                                         <  1:717100/150‑1 (MQ=255)
                                      ggCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTAGGTGAACGACAGTTCGCGGGTACTGCTGAAAATTGCCTATCTCAGCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTg                                                                                                        <  2:1161920/150‑1 (MQ=255)
                                      ggCCCACAACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCAGAGGACTGGTCATCAAACAGGAAGAACGACGGGTCGCGGTGATTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTg                                                                                                        <  1:432208/150‑1 (MQ=255)
                                        cccAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCgg                                                                                                                     >  1:1684406/1‑135 (MQ=255)
                                        cccAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCgg                                                                                                                     <  2:1684406/135‑1 (MQ=255)
                                        cccAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTggg                                                                                                      >  2:77164/1‑150 (MQ=255)
                                            ggACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACaa                                                                                                  <  1:1695978/150‑1 (MQ=255)
                                             gACGTCGAGCGCGCCCTGGACGCCGCGCATGCCACCGTCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAAc                                                                                                 <  2:912862/150‑1 (MQ=255)
                                                gTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGc                                                                                              >  1:41932/1‑150 (MQ=255)
                                                 tCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCa                                                                                             >  1:200206/1‑150 (MQ=255)
                                                     gcgcgcCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGc                                                                                         <  1:876408/150‑1 (MQ=255)
                                                            cTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGa                                                                                                     <  1:1477150/131‑1 (MQ=255)
                                                            cTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGa                                                                                                     >  2:1477150/1‑131 (MQ=255)
                                                               gACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGaa                                                                               >  1:1149926/1‑150 (MQ=255)
                                                                    cgcgCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCATGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCtt                                                                          <  1:204699/150‑1 (MQ=255)
                                                                    cgcgCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCtt                                                                          <  1:1419367/150‑1 (MQ=255)
                                                                     gcgcATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACgg                                                                                               >  1:289850/1‑128 (MQ=255)
                                                                     gcgcATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACgg                                                                                               <  2:289850/128‑1 (MQ=255)
                                                                     gcgcATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGGACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGAGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTg                                                                         >  1:1087880/1‑150 (MQ=255)
                                                                            ccgccgccgAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTAAAGCAAGAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCg                                                                  >  2:1417916/1‑150 (MQ=255)
                                                                             cgccgccgAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGa                                                                 <  1:175613/150‑1 (MQ=255)
                                                                             cgccgccgAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGa                                                                 <  1:189138/150‑1 (MQ=255)
                                                                              gccgccgAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGAc                                                                <  1:1394372/150‑1 (MQ=255)
                                                                                 gtcggaGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTg                                                             >  1:1316228/6‑150 (MQ=255)
                                                                                        ccTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTgg                                                      <  2:333503/150‑1 (MQ=255)
                                                                                           ggggCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGccg                                                   <  2:1560866/150‑3 (MQ=255)
                                                                                             ggCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCg                                                 <  2:41932/150‑1 (MQ=255)
                                                                                              gCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCAACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGAc                                                                <  2:1441411/134‑1 (MQ=255)
                                                                                              gCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCAACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGAc                                                                >  1:1441411/1‑134 (MQ=255)
                                                                                              gCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCAACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGAc                                                                >  1:1441431/1‑134 (MQ=255)
                                                                                              gCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCAACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGAc                                                                <  2:1441431/134‑1 (MQ=255)
                                                                                                            cAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTAt                                  <  2:1149926/150‑1 (MQ=255)
                                                                                                                  cGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGcc                            <  2:606306/150‑1 (MQ=255)
                                                                                                                    tGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGc                                                                                         >  2:1549107/1‑87 (MQ=255)
                                                                                                                    tGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGc                                                                                         <  1:1549107/87‑1 (MQ=255)
                                                                                                                           gTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATc                   >  2:1538131/1‑150 (MQ=255)
                                                                                                                            tACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATcc                  <  1:999872/150‑1 (MQ=255)
                                                                                                                               tgctgAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATCcgcg               >  1:760767/1‑150 (MQ=255)
                                                                                                                                  tgAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAAc                                                                             >  1:359410/1‑85 (MQ=255)
                                                                                                                                  tgAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAAc                                                                             <  2:359410/85‑1 (MQ=255)
                                                                                                                                   gAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATCCGCGCCCa           <  2:200206/150‑1 (MQ=255)
                                                                                                                                      aagtGCCGACCTCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAACCTGGGACAACGGCAAGGCCATACGCGAACCCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATCCGCGCCCagga        <  2:1001000/147‑1 (MQ=255)
                                                                                                                                            cGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTCACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGACATCGGACCACTTCCGCTATTTTGCCGGTTGCATCCGCGCCCAGGAGGGTgg  >  2:52815/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GTCAATGGCCAGCCGATTGCCGAATTCCCGCGCTCCACGGCCCAGGACGTCGAGCGCGCCCTGGACGCCGCGCATGCCGCCGCCGAAGCCTGGGGCAAGACCTCGGTGCAAGACCGTGCGCGGGTACTGCTGAAAATTGCCGACCGCATCGAACAGAACCTGGAAGTGCTGGCGGTTACCGAAAGCTGGGACAACGGCAAGGCCATACGCGAAACCTTGAATGCCGACGTGCCGCTGGCAGCGGACCACTTCCGCTATTTTGCCGGTTGCATCCGCGCCCAGGAGGGTGG  >  NC_002947/3069078‑3069367

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: