Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 994,630 C→T 6.2% R68Q (CGG→CAG)  ssuC ← aliphatic sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913994,6300CT 6.2% 229.8 / 2.9 79R68Q (CGG→CAG) ssuCaliphatic sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (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 = 6.03e-01

AAATCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACC  >  NC_000913/994481‑994779
                                                                                                                                                     |                                                                                                                                                     
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:1058682/1‑150 (MQ=255)
                                                      aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                 <  2:399480/150‑1 (MQ=255)
                                                      aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                 <  1:333679/150‑1 (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                                                                        <  1:1113118/150‑1 (MQ=255)
                                                                               gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                        <  2:1283830/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)
                                                                                                                                                     cgccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCAcc  <  2:557747/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATCACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACC  >  NC_000913/994481‑994779

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: