Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,819,810 A→T 8.0% L454* (TTA→TAA)  ligB ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,819,8100AT 8.0% 222.1 / 5.6 83L454* (TTA→TAA) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base A (30/46);  new base T (2/5);  total (32/51)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCATTAAGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCA  >  NC_000913/3819661‑3819959
                                                                                                                                                     |                                                                                                                                                     
tGGCATTAAGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTta                                                                                                                                                       <  1:1005035/150‑1 (MQ=255)
   cATTAAGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTtaata                                                                                                                                                    >  2:1018505/1‑150 (MQ=255)
    agtattcgccgcCCGGGTTAGCGTTATTCCCATTTCCATTACGCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGAATCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTATGTAATTGCTCTGGCGTtaataa                                                                                                                                                   <  2:315888/144‑1 (MQ=255)
    aTTAAGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAACGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTtaataa                                                                                                                                                   >  1:1692110/1‑150 (MQ=255)
      tcaGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTACAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATaaaa                                                                                                                                                 <  1:413873/148‑1 (MQ=255)
        aGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAAGTGCTCTGGCGTTAATAAAAGc                                                                                                                                               >  2:1271571/1‑150 (MQ=255)
          cgccgcCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGGTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCa                                                                                                                                             <  1:1608495/150‑1 (MQ=255)
             cgccCGGGTTAGCGGTATTCCCATTGCCATCACCGAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgag                                                                                                                                          <  1:443624/150‑3 (MQ=255)
             cgccCGGGTTAGCGGTAATCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgag                                                                                                                                          <  2:598971/150‑3 (MQ=255)
              gccCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagc                                                                                                                                         <  2:1632354/150‑4 (MQ=255)
               cccGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagca                                                                                                                                        <  2:709801/150‑5 (MQ=255)
                 cGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGTCTGCTTACGAGCCAGATTAAACTGATGCCATATCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAg                                                                                                                                      <  2:1833401/150‑1 (MQ=255)
                 cGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCTGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAg                                                                                                                                      <  1:291786/150‑1 (MQ=255)
                 cGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagcaat                                                                                                                                      <  1:1157835/150‑7 (MQ=255)
                    gTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAGAAGGCTGCTTACGAGCCAGATTAAACTGATGCCAGAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGGGGTCTGTAATTGCTCTGGCGTTATTAAAAGGCAGGAATAGcta                                                                                                                                   >  1:132252/1‑150 (MQ=255)
                      tAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATg                                                                                                                                 <  1:1854931/150‑1 (MQ=255)
                         cGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCACGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTc                                                                                                                              <  1:288283/150‑1 (MQ=255)
                             aTTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAg                                                                                                                          >  2:1699797/1‑150 (MQ=255)
                              ttCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGc                                                                                                                         <  1:1495356/150‑1 (MQ=255)
                              ttCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGGTAATAAAAGCCAGGAAAAGATATGTTCAAAgt                                                                                                                         >  1:1914950/1‑149 (MQ=255)
                                cccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGa                                                                                                                       >  2:1653895/1‑150 (MQ=255)
                                cccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGa                                                                                                                       >  2:489962/1‑150 (MQ=255)
                                  cATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCGAAGCGATg                                                                                                                     <  1:1759148/150‑1 (MQ=255)
                                  cATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATg                                                                                                                     <  1:931211/150‑1 (MQ=255)
                                     tGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGt                                                                                                                  <  2:1366214/150‑1 (MQ=255)
                                      gCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTc                                                                                                                 <  1:1303930/150‑1 (MQ=255)
                                      gCCATCACCCAGCGAGAAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTc                                                                                                                 <  2:1930081/150‑1 (MQ=255)
                                       ccATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCt                                                                                                                >  1:355631/1‑150 (MQ=255)
                                        cATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCCGTTATTGCTCAGGCGTTAAATAAAGCCAGGTAAATATATGTTAAAAGAGATGAGTcgg                                                                                                               >  1:1913980/1‑148 (MQ=255)
                                         aTCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTGTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGa                                                                                                              >  2:435186/1‑150 (MQ=255)
                                           cACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATg                                                                                                            >  2:1809828/1‑150 (MQ=255)
                                            aCCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCtgctgc                                                                                                           >  2:1023449/1‑150 (MQ=255)
                                            aCCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGc                                                                                                           <  1:505060/150‑1 (MQ=255)
                                             cccAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCa                                                                                                          >  2:1783501/1‑150 (MQ=255)
                                             cccAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCa                                                                                                          >  1:1301168/1‑150 (MQ=255)
                                                aGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCACCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAgcg                                                                                                       >  2:1025357/1‑150 (MQ=255)
                                                    aGTAACAGGCTGCTTACGAGCCAGATTAAACTGATGACATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAgcgcgcg                                                                                                   <  1:1281572/150‑1 (MQ=255)
                                                    aGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAgcgcgcg                                                                                                   <  2:1083902/150‑1 (MQ=255)
                                                      tAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTACTAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGcc                                                                                                 <  1:978886/150‑1 (MQ=255)
                                                              tgctTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGcc                                                                                         <  1:1271229/150‑1 (MQ=255)
                                                                    cGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTAATTTTCGCGATCCCCGGCGTGTTCTGTAAATGCTCTGGCGTTAATAAAAGCCAGGAAAACATATGTTAAAAGCGATGAGTCTGATGCTCCTCGCGCCAAATGGCCTCACCaa                                                                                  >  2:1602225/1‑150 (MQ=255)
                                                                     gAGCCAGATTAAACTGATGCCATAGCTGCGCATTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTAAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCaa                                                                                  <  1:1259987/150‑1 (MQ=255)
                                                                        ccAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGc                                                                               <  2:712918/150‑1 (MQ=255)
                                                                        ccAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGc                                                                               <  2:402981/150‑1 (MQ=255)
                                                                        ccAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGc                                                                               <  2:406538/150‑1 (MQ=255)
                                                                         cAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGcc                                                                              <  2:1258756/150‑1 (MQ=255)
                                                                             ttAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGTTCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATcc                                                                          <  2:1331988/150‑1 (MQ=255)
                                                                             ttAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATcc                                                                          <  1:1912515/150‑1 (MQ=255)
                                                                             ttAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATcc                                                                          <  1:1699798/150‑1 (MQ=255)
                                                                              tAAACTGATGCCATAGCTGCGCACTTTTACTTATCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCa                                                                         <  2:1582039/150‑1 (MQ=255)
                                                                               aaaCTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGGCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAg                                                                        <  2:1440073/150‑1 (MQ=255)
                                                                                 aCTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGcc                                                                      <  2:1695622/150‑1 (MQ=255)
                                                                                       cccATATCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGGTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGaa                                                                <  2:1376573/149‑1 (MQ=255)
                                                                                         cATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGTCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAAcc                                                              <  1:1472164/150‑1 (MQ=255)
                                                                                         cATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTCGTATGTAATTGCTCTGGCGTTTATATAAGCCAGGAAAAGCTATGTTCAAAGGGATGAGTCTGATGCACCGCGCGCCAACCGGCCTCACCACTGCCATCCCGCCCAAgctcg                                                              >  1:1562353/1‑146 (MQ=255)
                                                                                         aaTAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAAcc                                                              <  2:289890/149‑1 (MQ=255)
                                                                                             gCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGtt                                                          <  2:851623/150‑1 (MQ=255)
                                                                                                 cgcACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGc                                                      <  2:1512949/150‑1 (MQ=255)
                                                                                                  gcACTTTTACTTTTCGCGATACCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGcc                                                     <  2:434260/150‑1 (MQ=255)
                                                                                                     cttttacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGccccc                                                  >  2:1627471/1‑150 (MQ=255)
                                                                                                      ttttacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCcca                                                 <  2:955390/150‑1 (MQ=255)
                                                                                                       tttacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccag                                                >  1:1659700/1‑150 (MQ=255)
                                                                                                         tacttttCGCGATCCCCGGCGTGTTGTGTAATTGCTCTGGCGTTTATAAAAGCCAGGAAAAGATATGTTCCAAGCGATGAGTCTGATGCATCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCGCccagcc                                              >  2:600037/1‑150 (MQ=255)
                                                                                                         tacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCCTCCAGCCCAAGAACCGGTTTTGCCCccagcc                                              >  1:1883828/1‑150 (MQ=255)
                                                                                                            ttttCGCGATCCCCGGCGTGGTCTGTAATTGCTCTGGCGTTAATAAACGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGa                                           <  1:1783502/150‑1 (MQ=255)
                                                                                                             tttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGAc                                          <  1:435186/150‑1 (MQ=255)
                                                                                                             tttCGCGATCCCCGGCGGGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGAc                                          <  1:890619/150‑1 (MQ=255)
                                                                                                              ttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACt                                         >  1:2040271/1‑150 (MQ=255)
                                                                                                               tCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAACGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTa                                        <  2:810320/150‑1 (MQ=255)
                                                                                                               tCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTa                                        <  2:802992/150‑1 (MQ=255)
                                                                                                                     tCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTg                                  <  1:1973153/150‑1 (MQ=255)
                                                                                                                       cccGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCTCGCCAACCGGCCTCACCAATGCCATCCCGCCCAAGAACCTGTTTTTCCCCCAGCCAGACTAAGCGTGaa                                >  2:192879/1‑150 (MQ=255)
                                                                                                                       cccGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTgca                                >  1:414214/1‑148 (MQ=255)
                                                                                                                       cccGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGaa                                >  1:722298/1‑150 (MQ=255)
                                                                                                                         cGGCGTGTGCTGTAATTGCTCTGGCGTAAATAAAAGCCAGGAAAAGATATCATCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAAt                              <  2:1107211/150‑1 (MQ=255)
                                                                                                                             gtgtTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGGCAGACTAAGCGTGAAATGAAc                          <  2:622649/150‑1 (MQ=255)
                                                                                                                               gtTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCAACAATGCCATCCAGCCCAAGAACCTGATTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTg                        >  1:756353/1‑150 (MQ=255)
                                                                                                                                   tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTcc                    >  1:597312/1‑150 (MQ=255)
                                                                                                                                   tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGAGAAATGAACTGTTcc                    >  1:1745517/1‑150 (MQ=255)
                                                                                                                                     tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTg                  >  2:1562277/1‑150 (MQ=255)
                                                                                                                                     tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTg                  <  1:1186220/150‑1 (MQ=255)
                                                                                                                                     tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCACAAGAACCTGTTTTGCCCCCAGCCAGACTTAGCGTGAAATGAACTGTTCCTg                  >  1:1340286/1‑150 (MQ=255)
                                                                                                                                       aTTGCTCTGGCGGTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTTCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGAc                <  2:275191/150‑1 (MQ=255)
                                                                                                                                          gCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACGAATGCCATCCAGCCCAAGCACCTGGGTTGCCCCCCGCCAGACGAAGCGGGACATGACCTGTTCCTGAcgac             >  2:317443/1‑147 (MQ=255)
                                                                                                                                           ctctGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCACAATAAACTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAAc            <  1:442257/150‑1 (MQ=255)
                                                                                                                                              tGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATc         >  1:816461/1‑150 (MQ=255)
                                                                                                                                                 cGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAACCATCAGa      >  2:1840022/1‑150 (MQ=255)
                                                                                                                                                     aataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCa  >  1:1102044/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCATTAAGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCA  >  NC_000913/3819661‑3819959

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