Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,477 A→C 6.9% R98R (CGT→CGG yciF ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,4770AC 6.9% 345.8 / 7.4 119R98R (CGT→CGGyciFputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (41/68);  new base C (4/5);  total (46/73)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

CGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCT  >  NC_000913/1315332‑1315625
                                                                                                                                                 |                                                                                                                                                    
cGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCa                                                                                                                                                  >  1:249632/1‑150 (MQ=255)
   tGGCCTTTTCTTCTTCGAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtg                                                                                                                                               <  2:1428573/150‑2 (MQ=255)
     gCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgat                                                                                                                                             <  2:16106/150‑4 (MQ=255)
     gCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCg                                                                                                                                             <  1:227582/150‑1 (MQ=255)
     gCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCg                                                                                                                                             >  1:1236806/1‑150 (MQ=255)
     gCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCAGCTTTACGGTAACCTAATTGGTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCg                                                                                                                                             <  2:1046727/150‑1 (MQ=255)
     gCCTTTTCTTCTTCCAGGGTTTCTTTAAGAAGCTTCGCTGCTTTAAGGTAACCTAATTGTTAAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTTCGGCAATCAGTGCGGCATCACGCACTTCg                                                                                                                                             <  2:1565249/150‑1 (MQ=255)
      ccTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgatt                                                                                                                                            <  2:660340/150‑5 (MQ=255)
      ccTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgatt                                                                                                                                            <  1:910545/150‑5 (MQ=255)
       cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTATCCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAAGCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACtgattg                                                                                                                                           <  1:288962/150‑6 (MQ=255)
       cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACtgattg                                                                                                                                           <  1:244297/150‑6 (MQ=255)
       cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGtt                                                                                                                                           <  2:1426587/150‑1 (MQ=255)
       cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGGAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGtt                                                                                                                                           <  2:1122780/150‑1 (MQ=255)
        ttttCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgattgc                                                                                                                                          <  2:1178724/150‑7 (MQ=255)
        ttttCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGttt                                                                                                                                          <  2:1581587/150‑1 (MQ=255)
            cttcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTtct                                                                                                                                      <  1:606789/150‑1 (MQ=255)
             ttcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTtctc                                                                                                                                     <  1:1387332/150‑1 (MQ=255)
             ttcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTtctc                                                                                                                                     <  2:317240/150‑1 (MQ=255)
              tcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCg                                                                                                                                    >  2:774339/1‑150 (MQ=255)
               cttcCCGGGGTGCTTTCAGAAGCTTCGCTGCTTTCCGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCgg                                                                                                                                   <  1:334629/150‑1 (MQ=255)
                 tcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTa                                                                                                                                 >  2:121754/1‑150 (MQ=255)
                   cAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACt                                                                                                                               >  1:70123/1‑150 (MQ=255)
                      ggTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGCCTTTCTGTGCTGCGGCTCTCAGTGCGGCAGCACGCAGTTCGTTTTTCTCGGTACtttg                                                                                                                            >  2:1164574/1‑149 (MQ=255)
                            tttcAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGAc                                                                                                                      >  2:410117/1‑150 (MQ=255)
                            tttcAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTATCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGAc                                                                                                                      >  1:755703/1‑150 (MQ=255)
                                aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCCCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCa                                                                                                                  >  1:1580829/1‑150 (MQ=255)
                                aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCa                                                                                                                  >  2:1195538/1‑150 (MQ=255)
                                aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCa                                                                                                                  >  2:708605/1‑150 (MQ=255)
                                aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCa                                                                                                                  <  1:1155035/150‑1 (MQ=255)
                                aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTCCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCa                                                                                                                  >  2:1097066/1‑150 (MQ=255)
                                  aaGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCAtt                                                                                                                <  2:771295/150‑1 (MQ=255)
                                   aGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTa                                                                                                               <  1:1326975/150‑1 (MQ=255)
                                    gCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAg                                                                                                              <  2:619602/150‑1 (MQ=255)
                                      ttCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGct                                                                                                            <  1:951844/150‑1 (MQ=255)
                                        cGCTGCTTTACGTTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTTAATGGCAATCTCAAAATGCTCGACTTTCTGTGCTGCTGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGTTACTTTCGATGACCTCATTAGcttc                                                                                                          <  1:221864/150‑1 (MQ=255)
                                         gctgctTTTCGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTGTCTCGGTACTTTCGATGACCTCATTAGcttct                                                                                                         <  1:47539/150‑1 (MQ=255)
                                          ctgctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGcttctt                                                                                                        <  1:1150195/150‑1 (MQ=255)
                                           tgctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCttcttc                                                                                                       <  2:140598/150‑1 (MQ=255)
                                           tgctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCttcttc                                                                                                       <  2:249632/150‑1 (MQ=255)
                                           tgctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCttcttc                                                                                                       <  1:201719/150‑1 (MQ=255)
                                               ttACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATa                                                                                                   <  1:1597306/150‑1 (MQ=255)
                                                tACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATaa                                                                                                  <  1:1546186/150‑1 (MQ=255)
                                                   ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTCGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGAc                                                                                               <  2:1280912/150‑1 (MQ=255)
                                                     tAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGAGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCt                                                                                             >  2:968504/1‑150 (MQ=255)
                                                       aCCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTc                                                                                           <  2:503783/150‑1 (MQ=255)
                                                        ccTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCTGCATGCCGCACGTCGTGGTTGGCGGTAGTGTCGATGAGCGCATGAGCTTCTGCTGTAAGACCGTcc                                                                                          >  2:275895/1‑150 (MQ=255)
                                                           aaTTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCAtt                                                                                       <  1:315711/150‑1 (MQ=255)
                                                           aaTTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCAtt                                                                                       <  1:1464592/150‑1 (MQ=255)
                                                            aTTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACATTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTg                                                                                      <  2:975528/150‑1 (MQ=255)
                                                             ttGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGc                                                                                     <  2:225118/150‑1 (MQ=255)
                                                                ttCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcac                                                                                  <  1:483345/150‑1 (MQ=255)
                                                                ttCAGCCAGCGTCGCAAATGTCCCGTAACTGTCAATCACATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGATTTTCTCGGTACTATCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCACTGCcac                                                                                  <  1:107327/150‑1 (MQ=255)
                                                                 tCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                 <  1:780229/150‑1 (MQ=255)
                                                                  cagccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcacac                                                                                >  1:283563/1‑150 (MQ=255)
                                                                  cagccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcacac                                                                                >  1:1335202/1‑150 (MQ=255)
                                                                    gccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGCCTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCAATAGCTTCTTCTATAAGACCTTCCATTGCCACACAt                                                                              >  2:1164871/1‑150 (MQ=255)
                                                                    gccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACAt                                                                              >  1:1675525/1‑150 (MQ=255)
                                                                     ccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACAtt                                                                             <  2:490552/150‑1 (MQ=255)
                                                                        gcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCa                                                                          <  1:1219964/150‑1 (MQ=255)
                                                                        gcGTCGCTAATGTCCCATAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCa                                                                          <  1:147686/150‑1 (MQ=255)
                                                                            cGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACCCATTTCATgcg                                                                      >  2:761343/1‑150 (MQ=255)
                                                                             gCTAATGTCCCGTAACTGGCAATCTCATAATGCACGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCCTTTTTCTCGGTACTTTCGATTACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATgcgc                                                                     >  1:1057074/1‑150 (MQ=255)
                                                                              cTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGCCCTTCCATTGCCACACATTTCATgcgcg                                                                    >  2:285504/1‑149 (MQ=255)
                                                                              cTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCt                                                                    >  2:1708654/1‑150 (MQ=255)
                                                                               tAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCtt                                                                   <  1:437198/150‑1 (MQ=255)
                                                                                aaTGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTGTTCTCGGTACTTTCGATGACCTCATTTGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTa                                                                  >  2:481261/1‑150 (MQ=255)
                                                                                     cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCCCACATTTCATGCGCTTAAtttt                                                             >  1:1537674/1‑150 (MQ=255)
                                                                                      ccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATTACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAGTTTc                                                            >  1:182842/1‑150 (MQ=255)
                                                                                       cGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCa                                                           <  2:1529981/150‑1 (MQ=255)
                                                                                       cGTAACTGGCAATATCATAATGCTAGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCa                                                           <  2:1094373/150‑1 (MQ=255)
                                                                                            cTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTc                                                      >  1:1146030/1‑150 (MQ=255)
                                                                                             tGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCg                                                     <  2:70123/150‑1 (MQ=255)
                                                                                                cAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtt                                                  >  1:1312545/1‑150 (MQ=255)
                                                                                                 aaTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTc                                                 <  1:550986/150‑1 (MQ=255)
                                                                                                  aTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                <  1:634844/150‑1 (MQ=255)
                                                                                                  aTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                <  2:17285/150‑1 (MQ=255)
                                                                                                  aTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                <  2:1236806/150‑1 (MQ=255)
                                                                                                   tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTATCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCgg                                               <  1:373243/150‑1 (MQ=255)
                                                                                                   tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCgg                                               <  1:1554979/150‑1 (MQ=255)
                                                                                                        tAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCAGTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTcc                                          >  2:324011/1‑150 (MQ=255)
                                                                                                         aaTGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCca                                         <  2:1730892/150‑1 (MQ=255)
                                                                                                             cTCGACTTTCTTTGCTGCGGCAATCAGTGCGGCATCACGCAATTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                     <  1:759147/150‑1 (MQ=255)
                                                                                                             cTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                     >  1:937617/1‑150 (MQ=255)
                                                                                                             cTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                     >  1:199160/1‑150 (MQ=255)
                                                                                                              tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACt                                    <  1:708605/150‑1 (MQ=255)
                                                                                                              tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACt                                    <  1:1097066/150‑1 (MQ=255)
                                                                                                               cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTGTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGCTTCGGCCTCCACCACtt                                   >  1:287721/1‑150 (MQ=255)
                                                                                                                 aCTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGa                                 >  1:1311274/1‑150 (MQ=255)
                                                                                                                   tttCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                               <  2:932852/150‑1 (MQ=255)
                                                                                                                   tttCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                               <  1:338693/150‑1 (MQ=255)
                                                                                                                   tttCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                               <  2:34390/150‑1 (MQ=255)
                                                                                                                       tgtgCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATa                           >  2:422002/1‑150 (MQ=255)
                                                                                                                                gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                  <  1:1261886/150‑1 (MQ=255)
                                                                                                                                gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                  <  1:909575/150‑1 (MQ=255)
                                                                                                                                gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                  <  1:407715/150‑1 (MQ=255)
                                                                                                                                gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                  >  2:1695178/1‑150 (MQ=255)
                                                                                                                                        tGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATg          <  2:1388515/150‑1 (MQ=255)
                                                                                                                                          cGGCATCACGCACTTGGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAg        <  2:1226141/150‑1 (MQ=255)
                                                                                                                                            gCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATCCGCTAAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGtt      <  2:1335202/150‑1 (MQ=255)
                                                                                                                                                cACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATAGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCt  <  2:1588326/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCT  >  NC_000913/1315332‑1315625

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: