Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,933 T→A 5.2% I251N (ATC→AAC)  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,9330TA 5.2% 309.3 / 4.8 95I251N (ATC→AAC) nikCNi(2(+)) ABC transporter membrane subunit NikC
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base A (2/3);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGTATC  >  NC_000913/3616784‑3617082
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                                      <  2:1587852/138‑1 (MQ=255)
             cACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATc                                                                                                                                                      >  1:1587852/1‑138 (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                                                                         <  1:454037/150‑1 (MQ=255)
                                                                              gatTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGcaca                                                                         >  2:360308/1‑150 (MQ=255)
                                                                               atTAAAGAAGAGCGCAAGTATATCTGGACCCAGCAGCTGCAAATGTTCTGGACGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCacac                                                                        >  2:750793/1‑150 (MQ=255)
                                                                                ttAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                       >  1:998867/1‑150 (MQ=255)
                                                                                ttAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                       <  1:1788957/150‑1 (MQ=255)
                                                                                ttAACGACGCGCGCCAGTATATCTGGACACAGCCGCTGCAAATGTTCTGGACGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACt                                                                       >  2:1627671/1‑150 (MQ=255)
                                                                                     gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                                  >  1:1611526/1‑150 (MQ=255)
                                                                                     gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                                  <  1:803482/150‑1 (MQ=255)
                                                                                     gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                                  >  1:1279933/1‑150 (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)
                                                                                                                                                    aTCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGTAt   >  2:253807/1‑150 (MQ=255)
                                                                                                                                                     tCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGTATc  >  1:195584/1‑150 (MQ=255)
                                                                                                                                                     tCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGTATc  <  1:1728003/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTACACGGTGTATC  >  NC_000913/3616784‑3617082

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: