Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 793,336 C→T 6.8% M151I (ATG→ATA) ‡ modF ← molybdate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913793,3360CT 6.8% 320.0 / 7.1 102M151I (ATG→ATA) ‡modFmolybdate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (51/44);  new base T (4/3);  total (55/47)
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.71e-01

Only 100 of 103 total aligned reads displayed.

CTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGT  >  NC_000913/793189‑793485
                                                                                                                                                   |                                                                                                                                                     
atCCGGTATCTCTTCGAATTGATTTAGCACCAGTAGCAGATTAATACAGGAATGATGTATCGAGGCGTGTCGCTCAGCCAGATGCTGACGTGTGGGAACATCCAGTCCATCGAACGGCTCATCAAGAATAATCTAGTCAGGATCCGACAt                                                                                                                                                     <  2:1544061/149‑1 (MQ=255)
      gATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAgcgc                                                                                                                                               >  1:706585/1‑150 (MQ=255)
         ctcATCGAAGCGATTGAGCACCAGTCCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTGAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                            <  2:78982/150‑1 (MQ=255)
         ctcATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                            <  1:334741/150‑1 (MQ=255)
           cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                          <  2:1620659/150‑1 (MQ=255)
           cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                          <  2:1712765/150‑1 (MQ=255)
             tCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGacac                                                                                                                                        <  1:1466984/150‑1 (MQ=255)
             acGAAGCGATTCAGCACCAGTACCAGAGTAATACCGGACTGATGTACCGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGCCTCCGACATCAGCGCCTGacac                                                                                                                                        <  1:1066985/149‑1 (MQ=255)
              cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCATCAAGTCCGGCTCCGACATCAGCGCCTGacacc                                                                                                                                       >  1:23533/1‑149 (MQ=255)
              cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGAcaca                                                                                                                                       >  2:1213858/1‑150 (MQ=255)
              cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGAGCCTGAcaca                                                                                                                                       >  2:751425/1‑150 (MQ=255)
                   cGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                                  <  2:1278611/150‑1 (MQ=255)
                   cGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                                  >  2:540886/1‑150 (MQ=255)
                     tttGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGt                                                                                                                                <  2:1604763/149‑1 (MQ=255)
                     aTTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGt                                                                                                                                <  2:1189104/150‑1 (MQ=255)
                         tgCACCAGTACCAGATTAATACCGGACTGATGTAACGAGGCGAGGCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                                                                            <  2:328031/149‑1 (MQ=255)
                         aGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACATAGCAGGGttttt                                                                                                                            >  2:917000/1‑150 (MQ=255)
                         aGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                                                                            <  1:1548302/150‑1 (MQ=255)
                             ccAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACCTCCAGGCCATCGAACGGCTCATCAAGGATCAACAAGTCAGGCTCCGCCATCAGCGCCTGACACAGCCGGGTTTTTcgcg                                                                                                                        >  2:314981/1‑150 (MQ=255)
                               aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGtc                                                                                                                      >  1:1210168/1‑150 (MQ=255)
                               aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGtc                                                                                                                      >  1:183637/1‑150 (MQ=255)
                               aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCGCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGtc                                                                                                                      <  1:1461243/150‑1 (MQ=255)
                                  aCCAGCGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCg                                                                                                                   <  2:1531095/150‑1 (MQ=255)
                                   ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGc                                                                                                                  >  1:513031/1‑150 (MQ=255)
                                   ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGc                                                                                                                  >  2:735164/1‑150 (MQ=255)
                                        gTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTg                                                                                                             <  1:588042/150‑1 (MQ=255)
                                           aTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCATGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaa                                                                                                          >  1:1124172/1‑150 (MQ=255)
                                              ccGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAgg                                                                                                       >  2:52017/1‑150 (MQ=255)
                                                 gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                                    <  1:1296874/150‑1 (MQ=255)
                                                 gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                                    <  1:1666229/150‑1 (MQ=255)
                                                 gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                                    >  1:1794/1‑150 (MQ=255)
                                                   cTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAttt                                                                                                  <  2:369136/150‑1 (MQ=255)
                                                      aTGTAACGAGTCGAGTAGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCACCAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTaaa                                                                                               <  2:1259827/150‑1 (MQ=255)
                                                        gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGc                                                                                             <  1:657765/150‑1 (MQ=255)
                                                        gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGc                                                                                             <  1:494073/150‑1 (MQ=255)
                                                           aCGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtc                                                                                          <  2:632163/150‑1 (MQ=255)
                                                             gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtcgg                                                                                        >  1:1581523/1‑150 (MQ=255)
                                                             gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGGCTCGCCAGTGGAAAGGTATTTAAAGCgtcgg                                                                                        >  1:1561200/1‑150 (MQ=255)
                                                                gCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtcggtcg                                                                                     <  1:540886/150‑1 (MQ=255)
                                                                 cGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGGATTTAAAGCgtcggtcgg                                                                                    >  1:267771/1‑149 (MQ=255)
                                                                   aGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCTACAAGGCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGGCgagg                                                                                  <  1:426813/150‑1 (MQ=255)
                                                                      cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAATAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAATGGTATTTAAAGCGTCGGTCGaggagg                                                                               >  2:1448953/1‑150 (MQ=255)
                                                                      cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGGAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggcgg                                                                               >  2:1079373/1‑150 (MQ=255)
                                                                      cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAAAAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTGTTTCGCGTCTAGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggagg                                                                               >  2:1442155/1‑150 (MQ=255)
                                                                          cagccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCCGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAATGCGTCGGTCGAGGAGGGCgg                                                                           >  2:406716/1‑150 (MQ=255)
                                                                          cagccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCgg                                                                           <  2:529223/150‑1 (MQ=255)
                                                                            gccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGCCATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAGTTAAAGCGTCGGTCGAGGAGGGCGGTa                                                                         >  1:667548/1‑150 (MQ=255)
                                                                            gccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTa                                                                         <  1:514742/150‑1 (MQ=255)
                                                                             ccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTaa                                                                        <  2:685436/150‑1 (MQ=255)
                                                                                 ctgctgACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATAAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATAcc                                                                    <  1:419585/150‑1 (MQ=255)
                                                                                      gACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTTAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAAct                                                               >  2:213048/1‑150 (MQ=255)
                                                                                      gACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAAct                                                               >  2:1458932/1‑150 (MQ=255)
                                                                                        cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                             >  1:822026/1‑150 (MQ=255)
                                                                                        cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                             >  1:1574081/1‑150 (MQ=255)
                                                                                        cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATAAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                             >  1:827595/1‑150 (MQ=255)
                                                                                         gTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgct                                                            <  2:1272865/150‑1 (MQ=255)
                                                                                         gTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgct                                                            <  2:1558530/150‑1 (MQ=255)
                                                                                               cAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCa                                                      >  1:200942/1‑150 (MQ=255)
                                                                                                 aCATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGt                                                    >  1:1263348/1‑150 (MQ=255)
                                                                                                     ccAGGCCATCGAACGGCTCATCAAGAATCAACCAGTCAGGCTCCGCCATCAGCGCCTGACCCAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGCGGGCGGTAATACCGAACTGCTGCGCCCGTtgca                                                >  2:173395/1‑150 (MQ=255)
                                                                                                      cAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                               <  2:1124172/150‑1 (MQ=255)
                                                                                                      cAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                               <  1:651891/150‑1 (MQ=255)
                                                                                                       aGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatg                                              <  1:1552383/150‑1 (MQ=255)
                                                                                                          ccATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCaa                                           >  1:112386/1‑150 (MQ=255)
                                                                                                           cATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAAc                                          >  1:145981/1‑150 (MQ=255)
                                                                                                           cATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAAc                                          >  2:914785/1‑150 (MQ=255)
                                                                                                                atcGGCTCATCACGACTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTGTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGGCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGt                                     <  1:202454/148‑1 (MQ=255)
                                                                                                                aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGt                                     <  2:548559/150‑1 (MQ=255)
                                                                                                                  cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGc                                   <  1:1213858/150‑1 (MQ=255)
                                                                                                                   ggCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCa                                  >  1:1174894/1‑150 (MQ=255)
                                                                                                                        atctAGAATCAACAAGTCAGGCTCCGATATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCtt                             <  1:961334/150‑1 (MQ=255)
                                                                                                                        atcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCtt                             >  2:1016236/1‑150 (MQ=255)
                                                                                                                         tcaagaatcaagaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                            <  1:1669186/150‑1 (MQ=255)
                                                                                                                         tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCACCGCGTTGCAtgcaac                            <  2:706585/150‑6 (MQ=255)
                                                                                                                         tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                            >  1:741624/1‑150 (MQ=255)
                                                                                                                         tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                            >  1:1449808/1‑150 (MQ=255)
                                                                                                                           aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTAc                          >  2:622300/1‑150 (MQ=255)
                                                                                                                           aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGGATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGAGGTGCATCCTTTag                          >  1:251739/1‑149 (MQ=255)
                                                                                                                            aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                         <  2:1558416/150‑1 (MQ=255)
                                                                                                                            aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                         >  1:1163415/1‑150 (MQ=255)
                                                                                                                                ctgctgtgTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCAt                     >  2:597886/8‑150 (MQ=255)
                                                                                                                                 caacaaGTCCGGCTCCGACAGCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGGCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc                    <  2:183637/150‑1 (MQ=255)
                                                                                                                                 caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTTAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc                    >  2:164613/1‑150 (MQ=255)
                                                                                                                                 caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc                    <  2:1210168/150‑1 (MQ=255)
                                                                                                                                 caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc                    >  1:414981/1‑150 (MQ=255)
                                                                                                                                   acaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCt                  >  2:542976/1‑150 (MQ=255)
                                                                                                                                    tgtgTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTg                 >  1:514519/4‑150 (MQ=255)
                                                                                                                                      aGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGaa               >  1:74895/1‑150 (MQ=255)
                                                                                                                                       gTCAGGCTCCGACATCAGCGCCTGACATAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt              <  1:917000/150‑1 (MQ=255)
                                                                                                                                       gTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt              <  1:1458932/150‑1 (MQ=255)
                                                                                                                                       gTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt              >  1:1215600/1‑150 (MQ=255)
                                                                                                                                       gTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt              >  1:624839/1‑150 (MQ=255)
                                                                                                                                         cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCCTCCGGAATGa            >  2:450666/1‑150 (MQ=255)
                                                                                                                                         cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGa            >  1:1343677/1‑150 (MQ=255)
                                                                                                                                            gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtct         >  1:470830/1‑150 (MQ=255)
                                                                                                                                             cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTAAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtctc        <  1:673023/150‑1 (MQ=255)
                                                                                                                                               ccGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGCACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAg      >  2:1634057/1‑150 (MQ=255)
                                                                                                                                               ccGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAg      <  2:200942/150‑1 (MQ=255)
                                                                                                                                                 gACCTCAGCCCGTGTGACAGCAGGGCTTTTCTCGTCTCTCCAGTGGAAAGGTAATTATTGGGCCGGTCTTGGAGGGCTGTAATACCGTACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGcc    <  2:1070076/150‑1 (MQ=255)
                                                                                                                                                   cATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGt  >  2:1726091/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGT  >  NC_000913/793189‑793485

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