Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,566,646 T→A 10.5% K201* (AAA→TAA) ‡ yjiN ← DUF445 domain‑containing protein YjiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,566,6460TA10.5% 247.1 / 16.0 86K201* (AAA→TAA) ‡yjiNDUF445 domain‑containing protein YjiN
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (2/7);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

AATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCA  >  NC_000913/4566500‑4566795
                                                                                                                                                  |                                                                                                                                                     
aaTATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATc                                                                                                                                                    >  2:1477704/1‑150 (MQ=255)
aaTATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATc                                                                                                                                                    >  1:1192385/1‑150 (MQ=255)
 atatCATCAAGCAAAGAATCCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcg                                                                                                                                                   <  2:625713/150‑1 (MQ=255)
   atcatcAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcgcg                                                                                                                                                 >  2:843514/1‑150 (MQ=255)
     catcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgca                                                                                                                                               <  2:1419933/150‑2 (MQ=255)
     catcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgca                                                                                                                                               <  2:916589/150‑2 (MQ=255)
      atcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaa                                                                                                                                              <  2:851212/150‑3 (MQ=255)
       tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTgg                                                                                                                                             <  1:1385556/150‑1 (MQ=255)
       tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaag                                                                                                                                             <  1:1292854/150‑4 (MQ=255)
        caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGga                                                                                                                                            <  1:232864/150‑1 (MQ=255)
        caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaagt                                                                                                                                            <  2:1629608/150‑5 (MQ=255)
        caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaagt                                                                                                                                            <  1:1493806/150‑5 (MQ=255)
        caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaagt                                                                                                                                            <  1:483495/150‑5 (MQ=255)
         aaGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGgag                                                                                                                                           <  2:1288421/150‑1 (MQ=255)
          aGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGgaga                                                                                                                                          >  2:228483/1‑150 (MQ=255)
             aaaGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                       >  1:1138052/1‑150 (MQ=255)
                 aaTTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGa                                                                                                                                   <  1:229665/150‑1 (MQ=255)
                 aaTTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTAGGGAGAATTTTGGCTTTCAGGGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGa                                                                                                                                   <  2:807191/150‑1 (MQ=255)
                   ttCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATc                                                                                                                                 >  1:317846/1‑150 (MQ=255)
                    tCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCa                                                                                                                                <  1:1262789/150‑1 (MQ=255)
                    tCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCa                                                                                                                                >  2:1262789/1‑150 (MQ=255)
                      aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAAc                                                                                                                              <  1:1477704/150‑1 (MQ=255)
                        cgcgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtg                                                                                                                            <  1:1016008/150‑1 (MQ=255)
                         gcgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtgt                                                                                                                           >  1:1218353/1‑150 (MQ=255)
                           gTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGc                                                                                                                         <  2:1345489/150‑1 (MQ=255)
                             cAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtg                                                                                                                            >  1:1400062/1‑145 (MQ=255)
                             cAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtg                                                                                                                            <  2:1400062/145‑1 (MQ=255)
                             cAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGa                                                                                                                       >  1:1065040/1‑150 (MQ=255)
                              aGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGAt                                                                                                                      >  2:449259/1‑150 (MQ=255)
                               gAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATc                                                                                                                     >  2:451486/1‑150 (MQ=255)
                                 aaCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAg                                                                                                                   >  2:1363650/1‑150 (MQ=255)
                                   ccAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCTCTGGAGAAGGGCGATCCACTGTGCGATCcgcg                                                                                                                 >  2:192759/1‑150 (MQ=255)
                                   ccAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCg                                                                                                                 <  1:398741/150‑1 (MQ=255)
                                    cAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGt                                                                                                                >  1:1579293/1‑150 (MQ=255)
                                     aaCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTa                                                                                                               <  1:1202499/150‑1 (MQ=255)
                                         ccGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCca                                                                                                           <  2:330190/150‑1 (MQ=255)
                                         ccGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCca                                                                                                           <  2:5786/150‑1 (MQ=255)
                                            cgcTATGTTCTCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagca                                                                                                        <  1:436790/150‑1 (MQ=255)
                                              cTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCagcagc                                                                                                      <  1:1012466/150‑1 (MQ=255)
                                                  gTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCt                                                                                                  >  2:420191/1‑150 (MQ=255)
                                                   ttCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtg                                                                                                 <  1:1078541/150‑1 (MQ=255)
                                                     cGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCCCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgat                                                                                               <  1:1023900/150‑1 (MQ=255)
                                                      gCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                              <  2:1431472/150‑1 (MQ=255)
                                                      gCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                              <  1:536365/150‑1 (MQ=255)
                                                      gCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                              >  1:1134998/1‑150 (MQ=255)
                                                      gCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                              <  1:235992/150‑1 (MQ=255)
                                                      gCCCAACCATTCGGGGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCGCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                              <  2:939191/150‑1 (MQ=255)
                                                        ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGAc                                                                                            <  1:1569628/150‑1 (MQ=255)
                                                             cATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCg                                                                                       <  1:201337/150‑1 (MQ=255)
                                                                    tGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTgg                                                                                <  1:449259/150‑1 (MQ=255)
                                                                     ggggAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGt                                                                               <  2:416680/150‑1 (MQ=255)
                                                                         agaATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCAAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATa                                                                           <  2:732604/150‑1 (MQ=255)
                                                                          gaATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATAc                                                                          >  1:1370000/1‑150 (MQ=255)
                                                                            aTTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATActc                                                                        >  1:986276/1‑150 (MQ=255)
                                                                            aTTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATActc                                                                        >  1:131379/1‑150 (MQ=255)
                                                                               ttGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAggg                                                                                                                                      <  1:715266/85‑1 (MQ=255)
                                                                               ttGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAggg                                                                                                                                      >  2:715266/1‑85 (MQ=255)
                                                                                   cTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCa                                                                 >  2:1180532/1‑150 (MQ=255)
                                                                                   cTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCa                                                                 >  2:1090641/1‑150 (MQ=255)
                                                                                   cTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCa                                                                 <  2:1065040/150‑1 (MQ=255)
                                                                                           ggATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCa                                                         >  1:62305/1‑150 (MQ=255)
                                                                                               gctcgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTgg                                                     <  2:692152/150‑1 (MQ=255)
                                                                                                        ccagccagCGAAGAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGAGATCAACTGTGCGATCAGCGTATCCAGCAGCAGCTGATGACGATAGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGaaa                                            >  2:728882/1‑150 (MQ=255)
                                                                                                         cagccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAg                                           <  1:775390/150‑1 (MQ=255)
                                                                                                         cagccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAg                                           >  1:1040343/1‑150 (MQ=255)
                                                                                                          agccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGa                                          <  1:451486/150‑1 (MQ=255)
                                                                                                            ccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATc                                        <  2:170302/150‑1 (MQ=255)
                                                                                                            ccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATc                                        <  1:602135/150‑1 (MQ=255)
                                                                                                             cagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCg                                       >  2:1243356/1‑150 (MQ=255)
                                                                                                               gcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAc                                     >  2:748933/1‑150 (MQ=255)
                                                                                                                     aaTTTGCTGTGCTTTAAACTTGCTCTATTTTTCGCGCTGGAGATGGGCGAGCAACTGTGCGCTCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCACCGCCCTGGGGCCGGAAAGATCGACCTTATc                               <  2:1513530/150‑1 (MQ=255)
                                                                                                                      aTTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCa                              >  1:1277623/1‑150 (MQ=255)
                                                                                                                         tGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATc                           >  1:447333/1‑150 (MQ=255)
                                                                                                                         tGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATc                           >  2:709687/1‑150 (MQ=255)
                                                                                                                          gCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCg                          <  1:843514/150‑1 (MQ=255)
                                                                                                                           cTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGc                         >  1:1339695/1‑150 (MQ=255)
                                                                                                                           cTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGc                         >  1:1318633/1‑150 (MQ=255)
                                                                                                                             gggCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGccc                       >  2:1000938/1‑150 (MQ=255)
                                                                                                                                cGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGAt                    >  2:209332/1‑150 (MQ=255)
                                                                                                                                  aTAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATgg                  >  2:97687/1‑150 (MQ=255)
                                                                                                                                     aaCTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCTACCTTATCAATCGCCCGATGGAcc               >  1:930599/1‑150 (MQ=255)
                                                                                                                                      aCTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCCCTGGTGCCGGAAAGATCGACCTTATCACTCGCCCGATGGCCcg              >  1:1383318/1‑150 (MQ=255)
                                                                                                                                       ccagcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgc             >  1:292106/4‑150 (MQ=255)
                                                                                                                                       ccagcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgc             >  1:1404508/4‑150 (MQ=255)
                                                                                                                                       cTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgc             >  1:1061379/1‑150 (MQ=255)
                                                                                                                                         tGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgcgc           >  1:101143/1‑150 (MQ=255)
                                                                                                                                                tttATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAg    <  2:1339695/150‑1 (MQ=255)
                                                                                                                                                tttATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAg    <  2:1318633/150‑1 (MQ=255)
                                                                                                                                                 ttATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGc   <  2:1061379/150‑1 (MQ=255)
                                                                                                                                                  tATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCa  <  1:420191/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCA  >  NC_000913/4566500‑4566795

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