Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,929 T→A 5.3% F250I (TTT→ATT) ‡ 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,9290TA 5.3% 305.5 / 5.1 95F250I (TTT→ATT) ‡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 = 5.68e-01

TCGGCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTAC  >  NC_000913/3616780‑3617072
                                                                                                                                                     |                                                                                                                                               
tCGGCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGt                                                                                                                                                 <  1:1505113/150‑1 (MQ=255)
 cGGCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGtt                                                                                                                                                >  1:220837/1‑150 (MQ=255)
 cGGCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGtt                                                                                                                                                >  2:1728003/1‑150 (MQ=255)
   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:1215268/150‑6 (MQ=255)
              ctgcACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGATAAACAGCgccagc                                                                                                                                   <  2:133550/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                                                                                             >  1:1805263/1‑150 (MQ=255)
                                                    gACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCt                                                                                             <  2:1264395/150‑1 (MQ=255)
                                                      ccGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCAGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTgg                                                                                           >  1:743582/1‑150 (MQ=255)
                                                        gcgcCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGAc                                                                                         <  1:282595/150‑1 (MQ=255)
                                                        gcgcCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGAc                                                                                         >  1:196590/1‑150 (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:803482/150‑1 (MQ=255)
                                                                                         gacgCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGc                                                        >  1:1611526/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)
                                                                                                                                                     |                                                                                                                                               
TCGGCCATATGATGCTGCACGTCGCGGGGATGTCTTTCCTTGGCCTCGGTGTGACCGCGCCGACCGCCGAATGGGGCGTGATGATTAACGACGCGCGCCAGTATATCTGGACCCAGCCGCTGCAAATGTTCTGGCCGGGGCTGGCGCTGTTTATCAGCGTGATGGCCTTTAACCTGGTGGGTGACGCACTGCGCGATCATCTGGACCCTCATCTGGTGACGGAGCACGCACACTAATGCCGCAACAGATTGAACTACGTAATATCGCGCTACAGGCCGCGCAGCCGCTGGTAC  >  NC_000913/3616780‑3617072

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: