Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 966,301 T→C 5.9% V661A (GTA→GCA)  ycaI → conserved inner membrane protein YcaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913966,3010TC 5.9% 211.3 / 2.3 72V661A (GTA→GCA) ycaIconserved inner membrane protein YcaI
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/2);  total (34/38)
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.98e-01

TAGCAACGATAAAGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAAT  >  NC_000913/966152‑966450
                                                                                                                                                     |                                                                                                                                                     
tAGCAAGGATAAAGGAAATAACCATTCCTGGGGGGTTAAGGTTGATGACGGGACGAATAGCATGCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGGTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGt                                                                                                                                                       <  2:1567215/150‑1 (MQ=255)
   cAACGATAAAGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTAcc                                                                                                                                                    >  2:1073716/1‑150 (MQ=255)
         tAAAGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCa                                                                                                                                              >  2:1497551/1‑150 (MQ=255)
           aaGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATg                                                                                                                                            >  2:1272175/1‑150 (MQ=255)
            aGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATgg                                                                                                                                           >  1:1145227/1‑150 (MQ=255)
            aGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCAGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGCACCTCACCATgg                                                                                                                                           >  1:1566267/1‑150 (MQ=255)
             ggAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGc                                                                                                                                          <  2:529808/150‑1 (MQ=255)
               aaaTAACCATTCCTGTGTGGTTAAGGTTGATGTCGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAg                                                                                                                                        <  2:290019/150‑1 (MQ=255)
                aaTAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGt                                                                                                                                       <  1:822058/150‑1 (MQ=255)
                 aTAACCATTCCTGTGTGGTTAAGGTTGATGTCGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTa                                                                                                                                      <  2:1041799/150‑1 (MQ=255)
                  tAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGAATTCTTCTAACCGGTGATATTGAAGCCCAAGCTGAACAAAAGATGCTAAGCCGTTTCTGGCAGCAAGTGCAGGCAACATTGCTTCAGGAACCTCACCATGGCAGTaa                                                                                                                                     >  1:530365/1‑150 (MQ=255)
                    aCCATTCCTGTGTGGGTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATa                                                                                                                                   <  2:605018/150‑1 (MQ=255)
                         tCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCtca                                                                                                                              <  1:253455/150‑1 (MQ=255)
                                  gTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTg                                                                                                                     >  2:3085/1‑150 (MQ=255)
                                    tAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGcc                                                                                                                   >  1:463614/1‑150 (MQ=255)
                                       ggTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCAtt                                                                                                                >  1:1274873/1‑150 (MQ=255)
                                          tgatgaCGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAAt                                                                                                             >  1:158355/1‑150 (MQ=255)
                                             tgaCGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCa                                                                                                          >  1:1236582/1‑150 (MQ=255)
                                             tgaAGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCa                                                                                                          >  2:429/1‑150 (MQ=255)
                                                cGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCg                                                                                                       <  2:1184147/150‑1 (MQ=255)
                                                  ggACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGGACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAg                                                                                                     >  1:124193/1‑150 (MQ=255)
                                                       aaTAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAAt                                                                                                >  2:205280/1‑150 (MQ=255)
                                                        aTAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATg                                                                                               <  1:1031075/150‑1 (MQ=255)
                                                          aGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGa                                                                                             <  1:1073716/150‑1 (MQ=255)
                                                            cATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGaaa                                                                                           <  1:1413789/150‑1 (MQ=255)
                                                              ttcttcTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGGACCTCACAATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGGGAGTGAATGGAAAAg                                                                                         <  2:478100/150‑1 (MQ=255)
                                                                  tctAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGc                                                                                     <  1:1561165/150‑1 (MQ=255)
                                                                   ctAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTggaa                                                                                    >  2:712086/1‑148 (MQ=255)
                                                                    tAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCAc                                                                                   >  1:465012/1‑150 (MQ=255)
                                                                       ccGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCg                                                                                <  1:1477585/150‑1 (MQ=255)
                                                                             atatTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcgg                                                                          <  2:1216620/150‑1 (MQ=255)
                                                                              tatTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggc                                                                         <  2:530365/150‑1 (MQ=255)
                                                                                 tGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggcatc                                                                      >  1:322699/1‑150 (MQ=255)
                                                                                   aaGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATcgc                                                                    <  1:1429315/150‑1 (MQ=255)
                                                                                     gCCCCAGCTGAACAAAAGATCCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCAGTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCgcgc                                                                  <  1:878412/150‑1 (MQ=255)
                                                                                       cccAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCta                                                                >  1:713967/1‑150 (MQ=255)
                                                                                       cccAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCta                                                                >  1:811882/1‑150 (MQ=255)
                                                                                       cccAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCta                                                                <  1:1323147/150‑1 (MQ=255)
                                                                                         cAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGCACCTCACCATGGCAGTAATTCCTCATCTTGGTTGCCATTAATTCATCGAGTGAATGGAAAAGTGGCTCTCGGATCGGCATCGCGCtata                                                              <  2:1566267/150‑1 (MQ=255)
                                                                                           gCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAAc                                                            >  2:301027/1‑150 (MQ=255)
                                                                                              gAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTTGCACTCGCATCGGCATCGCGCTATAACGCa                                                         >  2:705619/1‑150 (MQ=255)
                                                                                              gAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCa                                                         >  2:802794/1‑150 (MQ=255)
                                                                                                 cAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATgg                                                      <  2:1018000/150‑1 (MQ=255)
                                                                                                   aaaGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGACACTCGCATCGGCATCGCGCTATAACGCATGGCg                                                    >  1:1082829/1‑150 (MQ=255)
                                                                                                         gCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGcc                                              <  2:1145227/150‑1 (MQ=255)
                                                                                                         gCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGcc                                              >  1:1087733/1‑150 (MQ=255)
                                                                                                          cTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                             <  2:1130367/150‑1 (MQ=255)
                                                                                                           tAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCct                                            >  2:871431/1‑150 (MQ=255)
                                                                                                           tAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCct                                            >  1:1360161/1‑150 (MQ=255)
                                                                                                            aaGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                           <  1:205280/150‑1 (MQ=255)
                                                                                                            aaGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                           >  2:915097/1‑150 (MQ=255)
                                                                                                            aaGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGCACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                           >  2:920602/1‑150 (MQ=255)
                                                                                                             aGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                          <  2:463614/150‑1 (MQ=255)
                                                                                                              gCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACATCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                         <  1:443089/150‑1 (MQ=255)
                                                                                                               ccGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                        <  1:429/150‑1 (MQ=255)
                                                                                                                    aCTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                   <  2:424112/150‑1 (MQ=255)
                                                                                                                     cTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACATCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGt                                  >  1:109822/1‑150 (MQ=255)
                                                                                                                               gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                        <  1:1593503/150‑1 (MQ=255)
                                                                                                                               gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                        <  1:782397/150‑1 (MQ=255)
                                                                                                                               gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                        <  1:871431/150‑1 (MQ=255)
                                                                                                                                    ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGGGACTGCCATCTAACAAAGTTAAGCATCGCTAACa                   >  1:1281266/1‑150 (MQ=255)
                                                                                                                                    ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCa                   >  1:940182/1‑150 (MQ=255)
                                                                                                                                     gCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCaa                  >  2:1505769/1‑150 (MQ=255)
                                                                                                                                       aaCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                <  1:1119112/150‑1 (MQ=255)
                                                                                                                                         cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc              <  1:1497551/150‑1 (MQ=255)
                                                                                                                                            tGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg           <  2:767689/150‑1 (MQ=255)
                                                                                                                                             gCTTCAGGTACCTGTCGATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg          <  2:1194225/150‑1 (MQ=255)
                                                                                                                                             gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg          <  2:713967/150‑1 (MQ=255)
                                                                                                                                             gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg          <  2:322699/150‑1 (MQ=255)
                                                                                                                                             gCTTCAGGCACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg          <  1:920602/150‑1 (MQ=255)
                                                                                                                                                tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGata       <  1:1272175/150‑1 (MQ=255)
                                                                                                                                                 cAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGatat      >  2:1396475/1‑150 (MQ=255)
                                                                                                                                                   ggTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa    >  2:830591/1‑150 (MQ=255)
                                                                                                                                                     tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAAt  >  2:222620/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGCAACGATAAAGGAAATAACCATTCCTGTGTGGTTAAGGTTGATGACGGGACGAATAGCATTCTTCTAACCGGTGATATTGAAGCCCCAGCTGAACAAAAGATGCTAAGCCGTTACTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAAT  >  NC_000913/966152‑966450

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