Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,673 A→C 6.0% E51A (GAG→GCG) ‡ rffG → dTDP‑glucose 4,6‑dehydratase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6730AC 6.0% 238.6 / 5.2 83E51A (GAG→GCG) ‡rffGdTDP‑glucose 4,6‑dehydratase 2
Reads supporting (aligned to +/- strand):  ref base A (32/44);  new base C (3/2);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972526‑3972821
                                                                                                                                                   |                                                                                                                                                    
gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  1:970151/1‑150 (MQ=255)
gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  1:1531464/1‑150 (MQ=255)
  aaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGCAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAgcgc                                                                                                                                                  >  1:1392650/1‑150 (MQ=255)
  aaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGATGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAgcgc                                                                                                                                                  >  2:411698/1‑150 (MQ=255)
   aaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCt                                                                                                                                                 <  2:1509236/150‑1 (MQ=255)
      ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                                                                                              <  2:1298084/150‑1 (MQ=255)
        cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                            >  2:1442145/1‑150 (MQ=255)
        cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                            <  1:746367/150‑1 (MQ=255)
         tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcgc                                                                                                                                           <  1:509214/150‑3 (MQ=255)
         tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcgc                                                                                                                                           <  1:1423412/150‑3 (MQ=255)
           aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                         >  2:558702/1‑150 (MQ=255)
            tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                        <  2:97233/150‑1 (MQ=255)
              aCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgag                                                                                                                                      <  2:970151/150‑1 (MQ=255)
                  gtggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAg                                                                                                                                  <  2:1155795/150‑1 (MQ=255)
                  gtggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAg                                                                                                                                  <  2:474898/150‑1 (MQ=255)
                   tggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCGTCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGt                                                                                                                                 <  2:866892/150‑1 (MQ=255)
                   tggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGt                                                                                                                                 <  1:16148/150‑1 (MQ=255)
                       gCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTTACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGat                                                                                                                             <  2:108976/150‑1 (MQ=255)
                       gCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCGCCGGTCGCGCGAAGGGGGCGCTTTGCCGTTGGGAAGGTTggt                                                                                                                             >  2:1548487/1‑148 (MQ=255)
                       gCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGat                                                                                                                             <  2:1355901/150‑1 (MQ=255)
                        ccGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGata                                                                                                                            >  1:1472729/1‑150 (MQ=255)
                         cGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGatat                                                                                                                           >  1:1177005/1‑150 (MQ=255)
                              ttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCg                                                                                                                      <  1:95535/150‑1 (MQ=255)
                                aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGCACCGGTCGCGCAAAGCCACCGCATTTCCTTTGAGAAAGTTGATATATGCCATc                                                                                                                   >  2:394984/1‑150 (MQ=255)
                                     cTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGc                                                                                                               >  1:1525793/1‑150 (MQ=255)
                                       cGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                                                                                             <  2:1218432/150‑1 (MQ=255)
                                          cgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAAc                                                                                                          <  2:303174/150‑1 (MQ=255)
                                            cTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTg                                                                                                        >  1:442399/1‑150 (MQ=255)
                                               gTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCa                                                                                                     <  2:511890/150‑1 (MQ=255)
                                                tGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       <  1:916437/115‑1 (MQ=255)
                                                tGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       >  2:916437/1‑115 (MQ=255)
                                                   gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                 <  2:536690/150‑1 (MQ=255)
                                                       tatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTAtt                                                                                             >  1:21989/1‑150 (MQ=255)
                                                        atcatgAACGAAACGAGCTACGTGGTGGTAGTGGTCGATAAGCTGACTTACGCCGGAAACCTGATGTAGCTGGCTCCGGTCGCGCAAAGCGTGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                                                                                            <  2:1490741/150‑1 (MQ=255)
                                                          catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc                                                                                          >  2:971770/1‑150 (MQ=255)
                                                            tcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGCACTGGCACGCGTATTCACTg                                                                                        >  1:574308/1‑150 (MQ=255)
                                                            tcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTg                                                                                        <  2:889195/150‑1 (MQ=255)
                                                              aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg                                                                                      >  2:1528083/1‑150 (MQ=255)
                                                              aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg                                                                                      <  1:1497904/150‑1 (MQ=255)
                                                               acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCAGGCACCGGTCGCGCAAAACGAGCGCTTTGCCATTGAGAAATTTGATATATGCGATCGGGCAGAACTGCCACGCCTATTAACTGAGc                                                                                     >  2:879168/1‑150 (MQ=255)
                                                                cgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCa                                                                                    <  1:1041441/150‑1 (MQ=255)
                                                                  aaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGCCAGAACTGGCACGCGTATTCACTGAGCATc                                                                                  <  1:1177077/150‑1 (MQ=255)
                                                                   aacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa                                                                                 <  2:212843/150‑1 (MQ=255)
                                                                    acgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAg                                                                                <  1:971770/150‑1 (MQ=255)
                                                                    acgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGc                                                                               <  1:394984/150‑1 (MQ=255)
                                                                    acgaCCGACGCGGTGGTAGTGGTCGATAATCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCTCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCTCTGAGCATCAg                                                                                <  2:531008/150‑1 (MQ=255)
                                                                          gACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGAc                                                                          >  1:609414/1‑150 (MQ=255)
                                                                           aCGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACt                                                                         >  1:1189736/1‑150 (MQ=255)
                                                                                gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTc                                                                    <  2:1525793/150‑1 (MQ=255)
                                                                                 ttgtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                   >  2:254335/3‑150 (MQ=255)
                                                                                 gggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                   <  2:1496450/149‑1 (MQ=255)
                                                                                  ggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                  <  1:1442145/150‑1 (MQ=255)
                                                                                      gTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCAt                                                              >  1:127552/1‑150 (MQ=255)
                                                                                       gggTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTTCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATc                                                             <  2:825961/149‑1 (MQ=255)
                                                                                           cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGc                                                         <  2:1531464/150‑1 (MQ=255)
                                                                                           cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACGGTGTCATGCATCTGGc                                                         >  1:697825/1‑150 (MQ=255)
                                                                                           cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCAAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGc                                                         >  1:1573346/1‑150 (MQ=255)
                                                                                            gATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGACCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGAGTCATGCATCTGGCa                                                        >  2:1161800/1‑150 (MQ=255)
                                                                                              tAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGc                                                      >  2:767258/1‑150 (MQ=255)
                                                                                              tAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGc                                                      <  2:1182974/150‑1 (MQ=255)
                                                                                                aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                    >  1:1134319/1‑150 (MQ=255)
                                                                                                 gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGGAGAACTGGCACGCGTATTAACAGAGCATCAGCCAGACTGTGTCAAGCATATGGCAGCCGa                                                   >  2:1419361/1‑150 (MQ=255)
                                                                                                       cTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATAGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCa                                             >  1:453744/1‑150 (MQ=255)
                                                                                                        tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                            >  2:674071/1‑150 (MQ=255)
                                                                                                         aCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATg                                           <  2:1392650/150‑1 (MQ=255)
                                                                                                                   ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                 <  1:1316267/150‑1 (MQ=255)
                                                                                                                   ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                 >  1:1126410/1‑150 (MQ=255)
                                                                                                                    cTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCt                                <  1:1046317/150‑1 (MQ=255)
                                                                                                                      gATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAt                              <  2:127552/150‑1 (MQ=255)
                                                                                                                       aTGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAtt                             <  1:1028825/150‑1 (MQ=255)
                                                                                                                          tCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTTCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGAc                          <  2:677409/150‑1 (MQ=255)
                                                                                                                          tCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGAc                          <  1:360174/150‑1 (MQ=255)
                                                                                                                                  aCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGc                  <  1:411698/150‑1 (MQ=255)
                                                                                                                                  aCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGc                  <  1:1184701/150‑1 (MQ=255)
                                                                                                                                    cGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                   <  2:546643/147‑1 (MQ=255)
                                                                                                                                    cGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                   >  1:546643/1‑147 (MQ=255)
                                                                                                                                     tcaaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCGGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGc               >  2:614726/7‑150 (MQ=255)
                                                                                                                                     ggTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGc               <  2:1189736/150‑1 (MQ=255)
                                                                                                                                      gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACAGCCCGGCAGCg              >  1:53729/1‑150 (MQ=255)
                                                                                                                                        cgcgCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt            >  1:1172865/1‑150 (MQ=255)
                                                                                                                                        aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt            >  1:1094523/4‑150 (MQ=255)
                                                                                                                                         aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt           >  1:1107087/3‑150 (MQ=255)
                                                                                                                                             aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTg       >  2:583182/1‑150 (MQ=255)
                                                                                                                                             aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTg       <  1:1416248/150‑1 (MQ=255)
                                                                                                                                              aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGa      <  1:254335/150‑1 (MQ=255)
                                                                                                                                               agcgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGaa     <  1:674071/150‑1 (MQ=255)
                                                                                                                                                gcgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGaaa    >  2:842121/1‑150 (MQ=255)
                                                                                                                                                  gagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAcc  <  2:1472729/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972526‑3972821

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