Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,422,556 A→T 6.0% I41F (ATC→TTC) ‡ yhdZ → putative amino acid ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,422,5560AT 6.0% 231.3 / 4.7 67I41F (ATC→TTC) ‡yhdZputative amino acid ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (2/2);  total (33/34)
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.53e-01

CCGGGCGTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGT  >  NC_000913/3422409‑3422704
                                                                                                                                                   |                                                                                                                                                    
ccGGGCGTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGAt                                                                                                                                                     <  4:343815/149‑1 (MQ=255)
      gTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCt                                                                                                                                               <  1:76602/149‑1 (MQ=255)
        acacCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCtgt                                                                                                                                             <  2:115653/149‑1 (MQ=255)
           ccGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCtgtgtg                                                                                                                                          >  2:462049/1‑149 (MQ=255)
                   aGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCa                                                                                                                                  <  1:126267/149‑1 (MQ=255)
                     gACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAgg                                                                                                                                >  1:273782/1‑149 (MQ=255)
                       cacTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGtt                                                                                                                              >  3:435723/1‑149 (MQ=255)
                        acTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTc                                                                                                                             >  2:103712/1‑149 (MQ=255)
                          tATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCg                                                                                                                           >  4:308479/1‑149 (MQ=255)
                              aGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTaa                                                                                                                       >  2:334288/1‑149 (MQ=255)
                                   aaTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGa                                                                                                                  >  3:171070/1‑149 (MQ=255)
                                    aTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGAc                                                                                                                 >  3:137003/1‑149 (MQ=255)
                                    aTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGAc                                                                                                                 <  2:96798/149‑1 (MQ=255)
                                      tttACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACaa                                                                                                               <  1:568508/149‑1 (MQ=255)
                                               aCCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGtt                                                                                                      >  1:134923/1‑149 (MQ=255)
                                               aCCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGtt                                                                                                      <  4:171070/149‑1 (MQ=255)
                                                 cTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGt                                                                                                    <  2:318168/149‑1 (MQ=255)
                                                  tGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTa                                                                                                   <  3:308479/149‑1 (MQ=255)
                                                         gcgATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAatca                                                                                            >  2:355986/1‑149 (MQ=255)
                                                           gatgatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAatcatc                                                                                          >  3:221339/1‑149 (MQ=255)
                                                           gatgatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAatcatc                                                                                          <  1:54426/149‑1 (MQ=255)
                                                             tgatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTg                                                                                        <  2:215906/149‑1 (MQ=255)
                                                             tgatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTg                                                                                        <  1:334288/149‑1 (MQ=255)
                                                              gatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTgg                                                                                       >  1:172277/1‑149 (MQ=255)
                                                               atTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGga                                                                                      >  1:451192/1‑149 (MQ=255)
                                                               atTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGga                                                                                      <  1:194498/149‑1 (MQ=255)
                                                                ttACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGgaa                                                                                     <  4:471356/149‑1 (MQ=255)
                                                                ttACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGgaa                                                                                     <  1:462049/149‑1 (MQ=255)
                                                                     cTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACa                                                                                >  2:148680/1‑149 (MQ=255)
                                                                         aaaaCGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCaa                                                                            >  3:20587/1‑149 (MQ=255)
                                                                              gTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAggg                                                                       >  4:12086/1‑149 (MQ=255)
                                                                               tCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGa                                                                      <  4:468624/149‑1 (MQ=255)
                                                                               tCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGAc                                                                     >  3:177011/1‑149 (MQ=255)
                                                                                cAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGAc                                                                     <  3:507403/149‑1 (MQ=255)
                                                                                cAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGAc                                                                     <  2:308138/149‑1 (MQ=255)
                                                                                 aaTAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACg                                                                    >  3:125995/1‑149 (MQ=255)
                                                                                  aTAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACgg                                                                   >  2:427463/1‑149 (MQ=255)
                                                                                  aTAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACgg                                                                   >  4:328781/1‑149 (MQ=255)
                                                                                                    tCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATc                                                 <  3:427367/149‑1 (MQ=255)
                                                                                                       aTGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTTTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGaa                                              <  4:314286/149‑1 (MQ=255)
                                                                                                       aTGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATAGCATCGaa                                              <  4:377389/149‑1 (MQ=255)
                                                                                                        tGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAAc                                             >  3:120944/1‑149 (MQ=255)
                                                                                                        tGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAAACATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAAc                                             >  4:142908/1‑149 (MQ=255)
                                                                                                            ttGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTaa                                         <  2:172277/149‑1 (MQ=255)
                                                                                                                aaaaTATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTTAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGaa                                     <  1:103712/149‑1 (MQ=255)
                                                                                                                aaaaTATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGaa                                     <  3:328781/149‑1 (MQ=255)
                                                                                                                aaaaTATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGaa                                     >  3:474270/1‑149 (MQ=255)
                                                                                                                  aaTATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGa                                   <  4:20587/149‑1 (MQ=255)
                                                                                                                            tAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAat                         >  2:424683/1‑149 (MQ=255)
                                                                                                                                 gTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAAc                                                                                 <  3:290480/88‑1 (MQ=255)
                                                                                                                                 gTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAAc                                                                                 >  4:290480/1‑88 (MQ=255)
                                                                                                                                 gTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGa                    <  4:137003/149‑1 (MQ=255)
                                                                                                                                     aaCCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAgcgc                >  2:333627/1‑149 (MQ=255)
                                                                                                                                     aaCCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAgcgc                <  2:127754/149‑1 (MQ=255)
                                                                                                                                      ccacacAGAACGATCCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAatcatc                                                                                          >  2:84773‑M1/7‑74 (MQ=255)
                                                                                                                                      ccacacAGAACGATCCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAatcatc                                                                                          <  1:84773‑M1/68‑1 (MQ=255)
                                                                                                                                        acacAGAACGATCCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTaa                                         <  3:457147‑M1/117‑1 (MQ=255)
                                                                                                                                        acacAGAACGATCCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTaa                                         >  4:457147‑M1/5‑121 (MQ=255)
                                                                                                                                          ggAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcag           <  2:451192/149‑1 (MQ=255)
                                                                                                                                             GAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggca        <  1:538389‑M1/149‑1 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAAt                                        <  1:151970‑M1/115‑1 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAAt                                        >  2:151970‑M1/1‑115 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGc                            >  3:206591‑M1/1‑127 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGc                            <  4:206591‑M1/127‑1 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg   >  3:339928/1‑149 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg   >  4:521804/1‑149 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg   >  4:53528/1‑149 (MQ=255)
                                                                                                                                                   aTCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGt  >  1:534603/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGGCGTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGT  >  NC_000913/3422409‑3422704

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