Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,422,555 G→A 6.2% R40R (CGG→CGA 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,5550GA 6.2% 204.0 / 3.3 64R40R (CGG→CGAyhdZputative amino acid ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (2/2);  total (30/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 = 7.03e-01

CCGGGCGTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAG  >  NC_000913/3422409‑3422703
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                <  2:96798/149‑1 (MQ=255)
                                    aTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGAc                                                                                                                >  3:137003/1‑149 (MQ=255)
                                      tttACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACaa                                                                                                              <  1:568508/149‑1 (MQ=255)
                                               aCCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGtt                                                                                                     <  4:171070/149‑1 (MQ=255)
                                               aCCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGtt                                                                                                     >  1:134923/1‑149 (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:194498/149‑1 (MQ=255)
                                                               atTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGga                                                                                     >  1:451192/1‑149 (MQ=255)
                                                                ttACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGgaa                                                                                    <  1:462049/149‑1 (MQ=255)
                                                                ttACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGgaa                                                                                    <  4:471356/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                                                                  >  4:328781/1‑149 (MQ=255)
                                                                                  aTAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACgg                                                                  >  2:427463/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:127754/149‑1 (MQ=255)
                                                                                                                                     aaCCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAgcgc               >  2:333627/1‑149 (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                                       >  2:151970‑M1/1‑115 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAAt                                       <  1:151970‑M1/115‑1 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGc                           <  4:206591‑M1/127‑1 (MQ=255)
                                                                                                                                               ACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGc                           >  3:206591‑M1/1‑127 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg  >  3:339928/1‑149 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg  >  4:521804/1‑149 (MQ=255)
                                                                                                                                                  gATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAg  >  4:53528/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGGGCGTACACCGCATTGAGGACACTATGAGCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAG  >  NC_000913/3422409‑3422703

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: