Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,931 T→A 5.3% F250L (TTT→TTA) ‡ nikC → Ni(2(+)) ABC transporter membrane subunit NikC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,616,9310TA 5.3% 304.3 / 5.1 95F250L (TTT→TTA) ‡nikCNi(2(+)) ABC transporter membrane subunit NikC
Reads supporting (aligned to +/- strand):  ref base T (45/45);  new base A (2/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGT  >  NC_000913/3616783‑3617079
                                                                                                                                                    |                                                                                                                                                    
gCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTa                                                                                                                                                     >  2:1587895/1‑150 (MQ=255)
 ccATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTAt                                                                                                                                                    >  2:1947103/1‑150 (MQ=255)
  cATATGATGCTGCACGTTGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATc                                                                                                                                                   <  1:1852778/150‑1 (MQ=255)
    tatGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAg                                                                                                                                                 <  2:1069560/150‑1 (MQ=255)
      tgatgCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCg                                                                                                                                               >  1:1802340/1‑150 (MQ=255)
       gatgCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGt                                                                                                                                              <  2:31574/150‑1 (MQ=255)
       gatgCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGt                                                                                                                                              >  1:1374622/1‑150 (MQ=255)
        atgCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTg                                                                                                                                             <  1:1587895/150‑1 (MQ=255)
         tgctgcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGa                                                                                                                                            >  2:1412065/1‑150 (MQ=255)
           ctgcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGATAAACAGCgccagc                                                                                                                                          <  2:133550/150‑6 (MQ=255)
           ctgcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGATAAACAGCgccagc                                                                                                                                          <  2:1215268/150‑6 (MQ=255)
           ctgcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGATAAACAGCgccagc                                                                                                                                          <  2:1207799/150‑6 (MQ=255)
             gcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACAGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGc                                                                                                                                        >  2:1287578/1‑150 (MQ=255)
              cACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATc                                                                                                                                                   >  1:1587852/1‑138 (MQ=255)
              cACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATc                                                                                                                                                   <  2:1587852/138‑1 (MQ=255)
               aCGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCt                                                                                                                                      <  2:53859/150‑1 (MQ=255)
                cGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCACCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCtt                                                                                                                                     >  2:1335597/1‑150 (MQ=255)
                 gTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCttt                                                                                                                                    >  1:1883963/1‑150 (MQ=255)
                   cgcgGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTaa                                                                                                                                  <  2:32151/150‑1 (MQ=255)
                    gcgGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAAc                                                                                                                                 <  1:1826908/150‑1 (MQ=255)
                      ggggATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCt                                                                                                                               <  2:1815300/150‑1 (MQ=255)
                         gATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCtggt                                                                                                                            <  2:402648/150‑1 (MQ=255)
                                  ccTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGc                                                                                                                   >  1:1264395/1‑150 (MQ=255)
                                      ggCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTg                                                                                                               <  1:1335597/150‑1 (MQ=255)
                                        ccTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTgcg                                                                                                             >  1:1140043/1‑150 (MQ=255)
                                         cTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTgcgc                                                                                                            <  2:1883963/150‑1 (MQ=255)
                                         cTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTgcgc                                                                                                            >  1:454305/1‑150 (MQ=255)
                                          tCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGcgcg                                                                                                           >  2:590063/1‑150 (MQ=255)
                                          tCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGCGCTGGCGCTGTTTATCAGCCTGATGGCCTTTAACCTGGTGGGTGACGCACTGcgcg                                                                                                           >  2:1503459/1‑150 (MQ=255)
                                           cGGTGTGACCGCGCCGACGGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATATGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGa                                                                                                          <  1:600898/150‑1 (MQ=255)
                                           cGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGa                                                                                                          <  1:1337666/150‑1 (MQ=255)
                                              tgtgACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGatca                                                                                                       <  1:1320401/150‑1 (MQ=255)
                                               gtgACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGatcat                                                                                                      >  1:1007256/1‑150 (MQ=255)
                                                 gACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCt                                                                                                    <  2:1264395/150‑1 (MQ=255)
                                                 gACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCt                                                                                                    >  1:1805263/1‑150 (MQ=255)
                                                   ccGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCAGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTgg                                                                                                  >  1:743582/1‑150 (MQ=255)
                                                     gcgcCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGAc                                                                                                >  1:196590/1‑150 (MQ=255)
                                                     gcgcCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGAc                                                                                                <  1:282595/150‑1 (MQ=255)
                                                          gACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCa                                                                                           >  1:786891/1‑150 (MQ=255)
                                                            ccgccgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGTGCTGTTTATCATCGGGATTGCCTTTATCCTTTTGGGTGTCGAAATGCGCGATCATCTGGACCCTCATc                                                                                         <  2:1824111/150‑1 (MQ=255)
                                                            ccgccgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATc                                                                                         <  2:786891/150‑1 (MQ=255)
                                                             cgccgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCt                                                                                        <  2:1374622/150‑1 (MQ=255)
                                                             cgccgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGATGACGCACTGCGCGATCATCTGGACCCTCATCt                                                                                        >  1:1501748/1‑150 (MQ=255)
                                                              gccgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTg                                                                                       >  1:392019/1‑150 (MQ=255)
                                                                cgAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGt                                                                                     <  1:206503/150‑1 (MQ=255)
                                                                 gAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGGGACTCACTGCGCGATCATCTGGACCCTCATCTGGTg                                                                                    <  2:1828798/150‑1 (MQ=255)
                                                                  aaTGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGa                                                                                   <  1:1900077/150‑1 (MQ=255)
                                                                   aTGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGAc                                                                                  <  2:1533258/150‑1 (MQ=255)
                                                                     ggggCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTAACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACgg                                                                                <  1:17205/150‑1 (MQ=255)
                                                                       ggCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAg                                                                              <  2:1408369/150‑1 (MQ=255)
                                                                           tgatgaTTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAgcacg                                                                          >  2:1114847/1‑150 (MQ=255)
                                                                               gatTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGcaca                                                                      >  2:360308/1‑150 (MQ=255)
                                                                               gatTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGcaca                                                                      <  1:454037/150‑1 (MQ=255)
                                                                                atTAAAGAAGAGCGCAAGTATATCTGGACCCAGCAGCTGCAAATGTTCTGGACGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCacac                                                                     >  2:750793/1‑150 (MQ=255)
                                                                                 ttAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                    <  1:1788957/150‑1 (MQ=255)
                                                                                 ttAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                    >  1:998867/1‑150 (MQ=255)
                                                                                 ttAACGACGCGCGCCAGTATATCTGGACACAGCCGCTGCAAATGTTCTGGACGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                    >  2:1627671/1‑150 (MQ=255)
                                                                                      gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                               >  1:1279933/1‑150 (MQ=255)
                                                                                      gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                               >  1:1611526/1‑150 (MQ=255)
                                                                                      gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                               <  1:803482/150‑1 (MQ=255)
                                                                                       acgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGcc                                                              >  2:1159186/1‑150 (MQ=255)
                                                                                            cgcCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGATGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAAc                                                         <  2:1501748/150‑1 (MQ=255)
                                                                                             gcCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACa                                                        <  1:1412065/150‑1 (MQ=255)
                                                                                                aGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGAt                                                     >  1:659990/1‑150 (MQ=255)
                                                                                                  tataTCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTg                                                   >  2:53641/1‑150 (MQ=255)
                                                                                                    tatCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGaa                                                 >  1:1939561/1‑150 (MQ=255)
                                                                                                     atCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAAc                                                >  2:1907396/1‑150 (MQ=255)
                                                                                                       cTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTa                                              >  2:921537/1‑150 (MQ=255)
                                                                                                       cTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGACGCAACAGATTGaaaaa                                              >  2:849710/1‑147 (MQ=255)
                                                                                                       cTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCAAACTAATGCCGCAACAGATTGAACTa                                              <  1:1947103/150‑1 (MQ=255)
                                                                                                         ggACCCAGCTGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACg                                            >  2:222577/1‑150 (MQ=255)
                                                                                                         ggACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACg                                            <  2:607895/150‑1 (MQ=255)
                                                                                                               aGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAatat                                      >  2:695746/1‑150 (MQ=255)
                                                                                                                  cGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATcgc                                   <  2:1007256/150‑1 (MQ=255)
                                                                                                                    cTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCgcgc                                 <  2:1311776/150‑1 (MQ=255)
                                                                                                                      gCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTa                               >  1:1793948/1‑150 (MQ=255)
                                                                                                                        aaaTGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGCCGGAGCACGCACACGAATGCAGCAACAGATTGAACTACGGAATATCGCGCTACa                             >  2:819223/1‑150 (MQ=255)
                                                                                                                        aaaTGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACa                             >  2:1400965/1‑150 (MQ=255)
                                                                                                                         aaTGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAg                            <  1:1287578/150‑1 (MQ=255)
                                                                                                                           tGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGc                          <  1:142260/150‑1 (MQ=255)
                                                                                                                              tCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCcgc                       >  1:629726/1‑150 (MQ=255)
                                                                                                                               cTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCcgcg                      >  2:373359/1‑150 (MQ=255)
                                                                                                                                   ccGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGt                                           >  1:1032266/1‑125 (MQ=255)
                                                                                                                                   ccGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGt                                           <  2:1032266/125‑1 (MQ=255)
                                                                                                                                    cGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGcc                 <  2:629726/150‑1 (MQ=255)
                                                                                                                                     ggggCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCg                <  2:220837/150‑1 (MQ=255)
                                                                                                                                      gccatcacgcTGATAAACAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGc               >  1:1757606/8‑150 (MQ=255)
                                                                                                                                       ccatcacgcTGATAAACAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCt              >  1:1491483/7‑150 (MQ=255)
                                                                                                                                         cTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTgg            <  1:966435/150‑1 (MQ=255)
                                                                                                                                         cTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTgg            <  1:590063/150‑1 (MQ=255)
                                                                                                                                          tGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGt           <  2:998867/150‑1 (MQ=255)
                                                                                                                                           ggCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGc                          <  2:1814455/134‑1 (MQ=255)
                                                                                                                                           ggCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGc                          >  1:1814455/1‑134 (MQ=255)
                                                                                                                                            gcgcTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTac         <  2:1411977/150‑1 (MQ=255)
                                                                                                                                            gcgcTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTac         >  2:1264242/1‑150 (MQ=255)
                                                                                                                                                   ttATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGgttt  <  1:1264242/150‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGT  >  NC_000913/3616783‑3617079

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: