Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,839 T→C 7.6% Q10Q (CAA→CAG dppF ← dipeptide ABC transporter ATP binding subunit DppF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8390TC 7.6% 224.1 / 7.8 79Q10Q (CAA→CAGdppFdipeptide ABC transporter ATP binding subunit DppF
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base C (4/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

CCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGC  >  NC_000913/3702690‑3702985
                                                                                                                                                     |                                                                                                                                                  
ccGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                                                                    >  2:1206790/1‑150 (MQ=255)
ccGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                                                                    <  2:1250984/150‑1 (MQ=255)
   cAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                                                                 <  2:1163924/150‑1 (MQ=255)
     gCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGttgttg                                                                                                                                               <  1:1197925/150‑1 (MQ=255)
        tgATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCATCGCTTTAACCAGACTTTCCGGCTCGAACATGCCTTTATTCACCTTATAATGTTTTTTCAGGTCGATATCCTGCAACAGCGGTTGTTGCAg                                                                                                                                            <  2:342048/149‑1 (MQ=255)
         gATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAgg                                                                                                                                           >  1:133790/1‑150 (MQ=255)
         gATTAGCCCACTACTGCCAGCGTTTTGCAACGTTCAAGGTTAAACGAAAAGCCATCCAGCTCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                                                                           <  1:868001/150‑1 (MQ=255)
          aTTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                                                                          <  1:1209196/150‑2 (MQ=255)
            tCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTg                                                                                                                                        <  1:979291/150‑1 (MQ=255)
             cGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTgg                                                                                                                                       <  2:391197/150‑1 (MQ=255)
             cGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTgg                                                                                                                                       <  1:1332730/150‑1 (MQ=255)
              gCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGc                                                                                                                                      >  1:593252/1‑150 (MQ=255)
               cccACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                                                                     >  1:1013072/1‑150 (MQ=255)
                       tGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGt                                                                                                                             <  1:713492/150‑1 (MQ=255)
                        gCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                                                                            >  1:23534/1‑150 (MQ=255)
                        gCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                                                                            >  2:923810/1‑150 (MQ=255)
                         ccAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                                                                           >  2:245181/1‑150 (MQ=255)
                         ccAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                                                                           >  1:375072/1‑150 (MQ=255)
                         ccAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                                                                           <  1:376334/150‑1 (MQ=255)
                           aGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTc                                                                                                                         <  1:1081754/150‑1 (MQ=255)
                            gCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                                                                        <  1:1429334/150‑1 (MQ=255)
                              gTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATa                                                                                                                      >  1:1073826/1‑150 (MQ=255)
                               ttttGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                                                                     <  2:133790/150‑1 (MQ=255)
                                    ccACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                                                                                <  1:605920/150‑1 (MQ=255)
                                          tCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTcc                                                                                                          >  1:1260570/1‑150 (MQ=255)
                                              ggTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                                                                      >  2:492165/1‑150 (MQ=255)
                                               gTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                                                                     <  1:473481/150‑1 (MQ=255)
                                               gTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                                                                     <  2:1109687/150‑1 (MQ=255)
                                                ttAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGCTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                                                                    >  1:945999/1‑150 (MQ=255)
                                                  aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                                                                  <  1:225749/150‑1 (MQ=255)
                                                  aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                                                                  <  1:59140/150‑1 (MQ=255)
                                                   aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatcat                                                                                                 >  2:165843/1‑150 (MQ=255)
                                                   aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatcat                                                                                                 <  1:923810/150‑1 (MQ=255)
                                                    acgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                                                                                <  1:992009/150‑1 (MQ=255)
                                                       aaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAg                                                                                             <  2:304152/150‑1 (MQ=255)
                                                         acgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTg                                                                                           <  1:5547/150‑1 (MQ=255)
                                                          cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTCGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                          >  1:1322415/1‑150 (MQ=255)
                                                           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTggg                                                                                         <  1:1559145/150‑1 (MQ=255)
                                                           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTggg                                                                                         >  2:81638/1‑150 (MQ=255)
                                                            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgggg                                                                                        >  1:837548/1‑149 (MQ=255)
                                                            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGt                                                                                        <  1:1354379/150‑1 (MQ=255)
                                                                 cAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATgg                                                                                   >  2:1279205/1‑150 (MQ=255)
                                                                       tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                                                                             >  1:1035285/1‑150 (MQ=255)
                                                                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTgg                                                                            >  1:403238/1‑150 (MQ=255)
                                                                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTcc                                                                                                          <  2:1100257/118‑1 (MQ=255)
                                                                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTcc                                                                                                          >  1:1100257/1‑118 (MQ=255)
                                                                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGAc                                                                          <  1:1424580/150‑1 (MQ=255)
                                                                           aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                                                                         <  2:375072/150‑1 (MQ=255)
                                                                            ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTg                                                                        >  1:414039/1‑150 (MQ=255)
                                                                                 cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                                                                   >  1:306136/1‑150 (MQ=255)
                                                                                          gcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGa                                                          >  1:1118749/1‑150 (MQ=255)
                                                                                            gAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGc                                                        <  2:593252/150‑1 (MQ=255)
                                                                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt                                                   <  2:1415448/150‑1 (MQ=255)
                                                                                                  gCCTTTCTTCACCGGATAATGTTTTTTCAGGGCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGCAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTc                                                  <  1:1038240/150‑1 (MQ=255)
                                                                                                        cttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGcc                                            <  2:886180/150‑1 (MQ=255)
                                                                                                        cttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGcc                                            >  1:1208829/1‑150 (MQ=255)
                                                                                                           cACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGt                                         <  2:23534/150‑1 (MQ=255)
                                                                                                            aCCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGtt                                        <  2:1013072/150‑1 (MQ=255)
                                                                                                            aCCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGtt                                        >  1:1256555/1‑150 (MQ=255)
                                                                                                               ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                                                                              <  1:934886/77‑1 (MQ=255)
                                                                                                               ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                                                                              >  2:934886/1‑77 (MQ=255)
                                                                                                               ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttctt                                     >  2:1424145/1‑150 (MQ=255)
                                                                                                                gATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttcttc                                    <  2:74453/150‑1 (MQ=255)
                                                                                                                gATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttcttc                                    <  1:1284180/150‑1 (MQ=255)
                                                                                                                      gTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGcgcg                              >  2:343708/1‑150 (MQ=255)
                                                                                                                      gTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGcgcg                              >  1:1114466/1‑150 (MQ=255)
                                                                                                                          tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACAt                          <  2:403238/150‑1 (MQ=255)
                                                                                                                           ttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATc                         >  2:1203499/1‑150 (MQ=255)
                                                                                                                               ggTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTc                     <  2:1118749/150‑1 (MQ=255)
                                                                                                                               ggTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTc                     <  2:414039/150‑1 (MQ=255)
                                                                                                                               ggTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTc                     >  2:146609/1‑150 (MQ=255)
                                                                                                                                gTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCa                    >  2:639850/1‑150 (MQ=255)
                                                                                                                                  ggccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGt                  >  1:1055180/7‑150 (MQ=255)
                                                                                                                                  cGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGt                  <  1:245181/150‑1 (MQ=255)
                                                                                                                                   gATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg                                                                 >  2:468085/1‑102 (MQ=255)
                                                                                                                                   gATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg                                                                 <  1:468085/102‑1 (MQ=255)
                                                                                                                                   gATAGCCTGCAACAGCGGTTGTCGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTg                 <  2:1322415/150‑1 (MQ=255)
                                                                                                                                    ccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTgg                >  2:3769/5‑150 (MQ=255)
                                                                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCAt             >  1:417394/2‑150 (MQ=255)
                                                                                                                                        ccTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATa            >  2:95177/1‑150 (MQ=255)
                                                                                                                                        ccTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATa            >  1:713102/1‑150 (MQ=255)
                                                                                                                                         cTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAg           >  2:449885/1‑150 (MQ=255)
                                                                                                                                                  cGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAgcgc  >  2:876236/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGC  >  NC_000913/3702690‑3702985

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: