Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 793,340 A→G 5.8% L150P (CTG→CCG)  modF ← molybdate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913793,3400AG 5.8% 305.4 / 7.0 104L150P (CTG→CCG) modFmolybdate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (53/44);  new base G (4/2);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-02

Only 100 of 106 total aligned reads displayed.

GATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTG  >  NC_000913/793195‑793489
                                                                                                                                                 |                                                                                                                                                     
gATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAgcgc                                                                                                                                                   >  1:706585/1‑150 (MQ=255)
   ctcATCGAAGCGATTGAGCACCAGTCCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTGAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                                <  2:78982/150‑1 (MQ=255)
   ctcATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                                <  1:334741/150‑1 (MQ=255)
     cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                              <  2:1712765/150‑1 (MQ=255)
     cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                              <  2:1620659/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:540886/1‑150 (MQ=255)
             cGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                                      <  2:1278611/150‑1 (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)
                       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)
                                     aTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaa                                                                                                              >  1:427268/1‑150 (MQ=255)
                                     aTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCATGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaa                                                                                                              >  1:1124172/1‑150 (MQ=255)
                                        ccGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAgg                                                                                                           >  2:52017/1‑150 (MQ=255)
                                           gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                                        <  1:1666229/150‑1 (MQ=255)
                                           gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                                        <  1:1296874/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)
                                                                cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAATAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAATGGTATTTAAAGCGTCGGTCGaggagg                                                                                   >  2:1448953/1‑150 (MQ=255)
                                                                cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGGAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggcgg                                                                                   >  2:1079373/1‑150 (MQ=255)
                                                                cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAAAAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTGTTTCGCGTCTAGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggagg                                                                                   >  2:1442155/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: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                                                   <  1:651891/150‑1 (MQ=255)
                                                                                                cAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                                   <  2:1124172/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)
                                                                                                                tcatcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATcc                                   >  2:1541777/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)
                                                                                                                                                cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGt   >  2:1315916/1‑150 (MQ=255)
                                                                                                                                                 aGCGCCTGACACAGGAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGTTCTAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg  <  1:1029424/150‑1 (MQ=255)
                                                                                                                                                 aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg  >  1:385035/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTG  >  NC_000913/793195‑793489

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: