Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 994,627 G→T 6.4% A69E (GCG→GAG)  ssuC ← aliphatic sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913994,6270GT 6.4% 211.1 / 4.7 79A69E (GCG→GAG) ssuCaliphatic sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (48/26);  new base T (2/3);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01

AAAATCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATG  >  NC_000913/994480‑994773
                                                                                                                                                   |                                                                                                                                                  
aaaaTCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGATTTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGcgcc                                                                                                                                                  <  2:19522/150‑1 (MQ=255)
 aaaTCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCGAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGcgccc                                                                                                                                                 >  1:118877/1‑150 (MQ=255)
 aaaTCACCAGCGGGATCAGCCCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTTTCCAGCAGCCGTTCTCCCCACCGAGATAGCCCGCTAATCAGTCCCAAAATCAACACCAGCGATCCGCCAATTGAAAAGCCAATCAGcgccc                                                                                                                                                 <  2:690638/150‑1 (MQ=255)
  aaTCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGcgcccg                                                                                                                                                <  1:392209/150‑1 (MQ=255)
      aCCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAACAGCCAATCAGCGCCCGCCAg                                                                                                                                            >  2:319883/1‑150 (MQ=255)
        cAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGa                                                                                                                                          >  2:269266/1‑150 (MQ=255)
            gggATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTg                                                                                                                                      <  1:480368/150‑1 (MQ=255)
               aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAtt                                                                                                                                   <  1:1188566/150‑2 (MQ=255)
               aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATc                                                                                                                                   >  2:979564/1‑150 (MQ=255)
                tCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCg                                                                                                                                  >  2:648551/1‑150 (MQ=255)
                 cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAttgg                                                                                                                                 <  1:1042183/150‑4 (MQ=255)
                 cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGc                                                                                                                                 >  2:785249/1‑150 (MQ=255)
                   gcgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCa                                                                                                                               >  1:495400/1‑150 (MQ=255)
                    cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCac                                                                                                                              >  1:453869/1‑149 (MQ=255)
                    cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAg                                                                                                                              >  2:303823/1‑150 (MQ=255)
                     gcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCCACATCGACCCCAGCGATCCGCCAATTGAAAAGCCCATCAGCGACCGCCAGGTGCTGATCGCCAGa                                                                                                                             >  2:191565/1‑150 (MQ=255)
                            gCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGcc                                                                                                                      <  2:1225683/150‑1 (MQ=255)
                               gCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCcaca                                                                                                                   >  2:534701/1‑150 (MQ=255)
                                  cGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCCACCCCAGCGATCCGCCCATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGtt                                                                                                                >  2:74797/1‑150 (MQ=255)
                                     tGCGCAACATCTGAATTGATGTATACAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCACTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                             <  1:138741/150‑1 (MQ=255)
                                     tGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                             <  1:143640/150‑1 (MQ=255)
                                       cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccc                                                                                                           >  1:258969/1‑149 (MQ=255)
                                       cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccg                                                                                                           >  2:32158/1‑150 (MQ=255)
                                           aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                       <  1:1237506/150‑1 (MQ=255)
                                           aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                       >  1:496279/1‑150 (MQ=255)
                                                tGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGgag                                                                                                  >  2:516998/1‑150 (MQ=255)
                                                     tGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGt                                                                                             >  2:735265/1‑150 (MQ=255)
                                                       aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                           <  1:333679/150‑1 (MQ=255)
                                                       aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                           <  2:399480/150‑1 (MQ=255)
                                                       aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                           >  1:1058682/1‑150 (MQ=255)
                                                        ggTATCGAGCAGCAGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCTATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                          <  2:484813/150‑1 (MQ=255)
                                                        ggTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                          <  1:1240465/150‑1 (MQ=255)
                                                             ccAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACg                                                                                     >  2:1109676/1‑150 (MQ=255)
                                                               agcagcCGTTCTCCCCAGCGAGATAGCCCGCTACTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGCAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGcc                                                                                   >  2:623981/1‑150 (MQ=255)
                                                                       gtcaccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGGCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGAAGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTcaccac                                                                           <  1:1002276/146‑1 (MQ=255)
                                                                         ctcCCAAGCGAGATAGCCAGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAAATGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTGGCCGCTGGCGGAGAGGGTCCAGCCCGCCGTCACCACTc                                                                         >  2:150666/1‑150 (MQ=255)
                                                                           ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCt                                                                       >  1:647256/1‑150 (MQ=255)
                                                                           ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                       <  1:419529/150‑1 (MQ=255)
                                                                              cAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTcc                                                                    >  2:680587/1‑150 (MQ=255)
                                                                                gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCCCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCgg                                                                  >  2:1145642/1‑150 (MQ=255)
                                                                                gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                  <  2:1283830/150‑1 (MQ=255)
                                                                                gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                  <  1:1113118/150‑1 (MQ=255)
                                                                                  gcgaTAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCATCCGGTg                                                                <  2:1119119/148‑1 (MQ=255)
                                                                                  gagaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTg                                                                <  1:1203682/150‑1 (MQ=255)
                                                                                   agaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAACCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGa                                                               <  1:1119041/150‑1 (MQ=255)
                                                                                    gaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                              <  1:724780/150‑1 (MQ=255)
                                                                                     aTAGCCCGCTTATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATGAGCGCCCGCCAGTAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACGACTCCTTCCGGTGAAg                                                             <  2:508868/150‑1 (MQ=255)
                                                                                     aTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAg                                                             <  2:519425/150‑1 (MQ=255)
                                                                                        gCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGGGAAGGcc                                                          >  1:620822/1‑149 (MQ=255)
                                                                                        gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                          >  1:1155344/1‑150 (MQ=255)
                                                                                        gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                          >  1:34213/1‑150 (MQ=255)
                                                                                        gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                          >  2:480548/1‑150 (MQ=255)
                                                                                         cccGCTAATCAGGCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCGGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCaa                                                         <  1:1316589/150‑1 (MQ=255)
                                                                                             cTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATa                                                     <  1:269266/150‑1 (MQ=255)
                                                                                                  cAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGt                                                >  2:477799/1‑150 (MQ=255)
                                                                                                      cccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCTGGTGAAGGCAAAATACGCGTCGAc                                            >  1:447409/1‑150 (MQ=255)
                                                                                                       ccAAAATCAACCCCAGCGCTCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGGGAAGGCAAAATACGCGTCGAcc                                           >  2:304202/1‑149 (MQ=255)
                                                                                                       ccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcc                                           >  2:277297/1‑149 (MQ=255)
                                                                                                        cAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcag                                          <  2:647256/150‑1 (MQ=255)
                                                                                                         aaaaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGCcagc                                         >  2:961216/1‑150 (MQ=255)
                                                                                                           aaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagcca                                       <  2:184904/150‑1 (MQ=255)
                                                                                                            aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGGGCTGATCGCCAGATGCTGCAACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCGTGTGCAGGCAAAATACGCGTCGAcagccag                                      >  2:1167636/1‑150 (MQ=255)
                                                                                                            aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccag                                      >  1:302837/1‑150 (MQ=255)
                                                                                                            aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccac                                      >  2:479478/1‑149 (MQ=255)
                                                                                                                  cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCAAcc                                >  2:30258/1‑150 (MQ=255)
                                                                                                                  cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAAcc                                >  2:77077/1‑150 (MQ=255)
                                                                                                                  cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTGCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCCCTCCTTCCGGTGAAGGCATACTACGCGACGCCAGCCAGCCAAcc                                >  2:215651/1‑150 (MQ=255)
                                                                                                                    cAGCGATCCGCCAAGTGAAACGCCAATCAGCGCCCGCCAGGAGCGGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGCGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCga                              <  1:1175630/150‑1 (MQ=255)
                                                                                                                           ccGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCCACCGAGGAGGcc                       >  2:730798/1‑150 (MQ=255)
                                                                                                                              ccAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGt                    >  1:1214279/1‑150 (MQ=255)
                                                                                                                               cAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGtt                   >  2:1126936/1‑150 (MQ=255)
                                                                                                                                 aTTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCCGCCAACCGAGGAGGCCAGTTGc                 >  1:972322/1‑150 (MQ=255)
                                                                                                                                  ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                >  1:497874/1‑150 (MQ=255)
                                                                                                                                  ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                <  2:729963/150‑1 (MQ=255)
                                                                                                                                   tGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCCACCGAGGAGGCCAGTTGccc               >  1:1435/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCcaca             >  1:387874/1‑150 (MQ=255)
                                                                                                                                      aaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCacac            <  1:648551/150‑1 (MQ=255)
                                                                                                                                       aaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACAcc           <  2:495400/150‑1 (MQ=255)
                                                                                                                                        aGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCg          >  1:1047991/1‑150 (MQ=255)
                                                                                                                                         ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCACCCGAGGAGGCCAGTTGCCACACCGc         >  2:647886/5‑150 (MQ=255)
                                                                                                                                         gCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACCCCGc         >  2:36688/1‑150 (MQ=255)
                                                                                                                                           cgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCACCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGccc       >  2:1135920/3‑149 (MQ=255)
                                                                                                                                                aGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGCCAAAAGACGCGTCGACAGCCCGCCAACCGAGGAGGCCCGGTGCCACCCCGCCACGATg  >  2:784463/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAATCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATG  >  NC_000913/994480‑994773

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