Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,857 C→T 14.0% intergenic (‑99/+16) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8570CT14.0% 164.7 / 13.4 64intergenic (‑99/+16)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (4/5);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCT  >  NC_000913/1987707‑1988006
                                                                                                                                                      |                                                                                                                                                     
atcaaacaacATAGGAACCCAATACTTGTTCGCCATGCTGCGGTTGGTACACTGTAACGCCCTTATACCATTTTCAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGcc                                                                                                                                                       <  2:512471/146‑1 (MQ=255)
      ccACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                                  <  1:1273597/150‑1 (MQ=255)
       cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACACCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCa                                                                                                                                                 <  1:424789/150‑1 (MQ=255)
        acaTCAGAAACCAATTCCTGATCTGCATTCAGCGGCTTGGACATCGTAAACTCCTGAAAGCATTTTGAAAGCCGTTATACTAGCGACTGCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAg                                                                                                                                                <  1:24664/150‑1 (MQ=255)
            cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGttt                                                                                                                                            <  1:925846/150‑1 (MQ=255)
               aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCTCCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttt                                                                                                                                         <  1:860084/150‑1 (MQ=255)
               aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttt                                                                                                                                         <  1:1027502/150‑1 (MQ=255)
                acccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTAATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTg                                                                                                                                        <  1:668383/148‑1 (MQ=255)
                 aCCAATTCCTGATCGGCACTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCATGACACAAAAAAACCTGCTCCGGCAGGtttttttg                                                                                                                                       <  1:659288/150‑3 (MQ=255)
                    aaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCAACCGTCATCAGGACACAAAAAAACCTGCCGGAGCTGCTTTTTTGTTAt                                                                                                                                    >  2:1086951/1‑150 (MQ=255)
                      ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCg                                                                                                                                  <  2:745614/150‑1 (MQ=255)
                      ttCCTGATCGGCATGCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttttgtgtcc                                                                                                                                  <  1:668349/150‑8 (MQ=255)
                       tCCTTATCTGCATTCATCGGCTGGTATATCGTAACCTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGGCACATCCACCGTCATCAGGACACAAAAAAACCTTCCGGAGCAGGTTTTTTGTTATCgg                                                                                                                                 <  2:61759/150‑1 (MQ=255)
                        ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTAACATCCACCTACATCAGGACACAAAAAAACCTGCCGGATCAGGTTTTTAGTGATCGGa                                                                                                                                >  2:556691/1‑150 (MQ=255)
                          tGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                              >  1:218636/1‑150 (MQ=255)
                          tGATCGGCAATCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                              <  1:1049477/150‑1 (MQ=255)
                           gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACa                                                                                                                             <  2:28057/150‑1 (MQ=255)
                           gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACa                                                                                                                             >  1:1262830/1‑150 (MQ=255)
                            aTCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACat                                                                                                                            >  1:546033/1‑150 (MQ=255)
                                    tCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTg                                                                                                                    <  1:396198/150‑1 (MQ=255)
                                     cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCCCCGTCATCAGGACACAAAAAAACCTGCCGGCGCAGGTGTTTTGTTATCGGAACATATTGCCTgg                                                                                                                   >  2:588891/1‑150 (MQ=255)
                                     cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTgg                                                                                                                   >  1:590671/1‑150 (MQ=255)
                                       gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCg                                                                                                                 <  2:289053/150‑1 (MQ=255)
                                           cTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGCCTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                             >  1:1108907/1‑150 (MQ=255)
                                              ggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTc                                                                                                          <  2:894478/150‑1 (MQ=255)
                                                aCATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                        >  1:601324/1‑150 (MQ=255)
                                                 cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                                       <  1:245946/150‑1 (MQ=255)
                                                  aTCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTg                                                                                                      >  2:1054444/1‑150 (MQ=255)
                                                       aaaCTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGt                                                                                                 >  2:971529/1‑150 (MQ=255)
                                                           tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGTTACGTCTTTGAACGTCTTg                                                                                             >  1:418684/1‑150 (MQ=255)
                                                            ccTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                            >  1:1107726/1‑150 (MQ=255)
                                                            ccTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                            <  2:998513/150‑1 (MQ=255)
                                                             cTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCa                                                                                           <  1:10736/150‑1 (MQ=255)
                                                                  aGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGt                                                                                      <  1:1044658/150‑1 (MQ=255)
                                                                       tttGATAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGTAACATATTGCATGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTg                                                                                 <  1:853378/150‑1 (MQ=255)
                                                                         tGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGaa                                                                               <  1:746682/150‑1 (MQ=255)
                                                                           aaaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATTGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACa                                                                             <  1:835551/150‑1 (MQ=255)
                                                                            aatccGTTATAGTAGCTACTTCATATCTTCAGCGATAGTCACATGCTCCGTCATAAGGACACAAATAAACCTTCCTGAGCAGGTTTTTTGTTATCGGAACATAATGCGTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACAt                                                                            <  1:179342/147‑1 (MQ=255)
                                                                             aGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATa                                                                           >  2:291705/1‑150 (MQ=255)
                                                                                cGTTATAGTAGCGAATTCACATATTCAGAGATAGTCACATCCACCGTAATCAGTACACAAAAAAACCTGCCTGAGCAGGTTTTTTGTTATCGTAACATATTGCCTGGCGGTACGTCTTTGATCGTCTTGCAATAGTTATTGAACATACtt                                                                        <  1:400402/150‑1 (MQ=255)
                                                                                    ataGTAGCGACTTCACATCTTCATCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                    <  1:724060/150‑1 (MQ=255)
                                                                                    ataGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTTCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                    >  2:855422/1‑150 (MQ=255)
                                                                                       gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                                 >  2:821981/1‑150 (MQ=255)
                                                                                            gACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTg                                                            >  2:103018/1‑150 (MQ=255)
                                                                                               ttCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgcgc                                                         <  1:971529/150‑1 (MQ=255)
                                                                                               ttCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgcgc                                                         <  1:399062/150‑1 (MQ=255)
                                                                                                     tCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTc                                                   <  2:509997/150‑1 (MQ=255)
                                                                                                                   cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTAACGGAACATATTCCCTGGCGGTACGTCTTTGATCGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTccc                                     >  2:93645/1‑149 (MQ=255)
                                                                                                                     caTCCACCGTCATCAGGACCCAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGGACGTCGTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                   <  1:824484/150‑1 (MQ=255)
                                                                                                                     caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                   <  1:884723/150‑1 (MQ=255)
                                                                                                                     caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGATATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTAATCGCGGCACTCCGAc                                   >  2:496907/1‑150 (MQ=255)
                                                                                                                        ccACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAt                                >  1:412835/1‑150 (MQ=255)
                                                                                                                         cACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAtt                               >  2:1032222/1‑150 (MQ=255)
                                                                                                                                 tcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTCCGGCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATAc                       <  1:1108323/150‑1 (MQ=255)
                                                                                                                                    ggACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGTATTTTGCGCAGTTGCATCGCTGCACTCCGACCATTTGTTATACAgg                    <  1:777025/150‑1 (MQ=255)
                                                                                                                                       cacaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                 <  2:418684/150‑1 (MQ=255)
                                                                                                                                          caaaaaaCCTGCTCCGGCAGGTTTTTTGTTATGGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCTGCACTCCGACCATTTGTTATACTGGTGTTAt              >  2:202032/2‑150 (MQ=255)
                                                                                                                                          caaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCATTAGTTATTGAACATACTTTTCAGGATTTTGCGCGGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt              >  1:31514/2‑150 (MQ=255)
                                                                                                                                           aaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt             >  2:910430/1‑150 (MQ=255)
                                                                                                                                           aaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCCCTCCGACCATTTGTTATACAGGTGTTAtt             >  1:1189832/1‑150 (MQ=255)
                                                                                                                                           aaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt             >  2:884006/1‑150 (MQ=255)
                                                                                                                                            aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg            >  1:577992/1‑150 (MQ=255)
                                                                                                                                              aaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTc          >  1:182117/1‑150 (MQ=255)
                                                                                                                                                   tGCCGGAGCAGGTTTTTTGTTATTGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGCCCATTTGTTATACAGGTGTTTTTGTCTTTgc     >  2:342222/1‑150 (MQ=255)
                                                                                                                                                    gCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTgcg    >  2:430627/1‑150 (MQ=255)
                                                                                                                                                      cGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCt  >  1:920593/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCT  >  NC_000913/1987707‑1988006

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