Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 916,431 A→T 7.1% intergenic (‑384/‑42) aqpZ ← / → ybjD aquaporin Z/putative OLD family ATP‑dependent endonuclease; DUF2813 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913916,4310AT 7.1% 213.5 / 4.3 69intergenic (‑384/‑42)aqpZ/ybjDaquaporin Z/putative OLD family ATP‑dependent endonuclease; DUF2813 family protein
Reads supporting (aligned to +/- strand):  ref base A (25/39);  new base T (3/2);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACGCGTGG  >  NC_000913/916283‑916571
                                                                                                                                                    |                                                                                                                                            
aataatTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTtata                                                                                                                                              <  1:144026/149‑1 (MQ=255)
 ataatTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATAc                                                                                                                                             <  6:345216/149‑1 (MQ=255)
  taatTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACt                                                                                                                                            <  8:247439/149‑1 (MQ=255)
       tATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCa                                                                                                                                        <  6:192276/148‑1 (MQ=255)
           aTGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATTCTGCCAGCGt                                                                                                                                    >  4:87376/1‑148 (MQ=255)
           aTGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATTCTGCCAGCGt                                                                                                                                    >  4:92346/1‑148 (MQ=255)
           aTGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGt                                                                                                                                    >  5:205238/1‑148 (MQ=255)
            tGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTa                                                                                                                                   >  7:99977/1‑148 (MQ=255)
              cTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAg                                                                                                                                <  6:23487/149‑1 (MQ=255)
                aaTTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGa                                                                                                                              >  8:210772/1‑149 (MQ=255)
                 attattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGAt                                                                                                                             <  6:205238/149‑1 (MQ=255)
                   tattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATa                                                                                                                            <  4:297212/148‑1 (MQ=255)
                   tattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATa                                                                                                                            >  3:297212/1‑148 (MQ=255)
                    attTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATaa                                                                                                                           <  1:300578/148‑1 (MQ=255)
                      ttGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGt                                                                                                                         >  3:132700/1‑148 (MQ=255)
                          ttttGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTcaca                                                                                                                     >  1:19453/1‑148 (MQ=255)
                             tgttgtTTTTGCAGACTTATCAGCAAGAGGGAGTATAACTCGATAATTCGCTCATTTTTCAGACATTTGACATGCTTAAATGTGAAGTCATAACGTATTAGCAAGGCCTTTCCCGTTATACTGCAAGCGTAAAGGATAAGTAACATAtt                                                                                                                 >  5:235383/1‑149 (MQ=255)
                               ttgttTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATAttt                                                                                                                <  4:90070/148‑1 (MQ=255)
                                    tttGCAGACTTATCATCAACACTGAGTATAACCCGATTATTCTCTCATTATTCAGACATTTCCCATGCTTAAATCTGATGTCATCAAGAATTAGCAAGGCCTTTCCCGTTATACTGCCATCGTAAAGCATAAGTCACATATTTCTGCAg                                                                                                          <  2:174029/149‑1 (MQ=255)
                                      tGCAGACTTATCAGCAAGAGGTAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAgg                                                                                                         <  1:321907/148‑1 (MQ=255)
                                      tGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAgg                                                                                                         <  7:379376/148‑1 (MQ=255)
                                      tGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAgg                                                                                                         <  3:8698/148‑1 (MQ=255)
                                          gACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGata                                                                                                    <  8:99977/149‑1 (MQ=255)
                                           aCTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGatat                                                                                                   <  4:166117/149‑1 (MQ=255)
                                               aTCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGa                                                                                                 <  7:327059/147‑1 (MQ=255)
                                                tCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGAtt                                                                                               <  7:111675/148‑1 (MQ=255)
                                                tCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTAGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCACGTTATTCTGCAAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTc                                                                                              >  8:375210/1‑149 (MQ=255)
                                                  aGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCtt                                                                                            <  7:210772/149‑1 (MQ=255)
                                                  aGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCtt                                                                                            <  8:211886/149‑1 (MQ=255)
                                                      agatGGAGTATATCGCGATTATTCGCTCATTTTTCAGACTTTTTCCATGCTTATATGTGATGTCATCTCGTATTAGCATGGCGTTTCCCGTTATACTGACAGCGTAAATTATAAGTCACATATTTCTGGAGGTTATATGATTCTTTAgc                                                                                        <  6:235383/149‑1 (MQ=255)
                                                      agagGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAgc                                                                                        <  3:128916/149‑1 (MQ=255)
                                                       gagGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAgcg                                                                                       >  8:212967/1‑149 (MQ=255)
                                                                    gcgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATtgtg                                                                          >  2:153189/1‑149 (MQ=255)
                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATTCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                                                         <  3:87376/149‑1 (MQ=255)
                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATTCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                                                         <  3:92346/149‑1 (MQ=255)
                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                                                         <  1:173781/149‑1 (MQ=255)
                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                                                         <  7:293538/149‑1 (MQ=255)
                                                                         tattCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTg                                                                                           <  8:387864/127‑1 (MQ=255)
                                                                                     tttCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTc                                                                                                               <  3:31547/95‑1 (MQ=255)
                                                                                     tttCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTc                                                                                                               >  4:31547/1‑95 (MQ=255)
                                                                                      ttCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAAc                                                        >  5:288915/1‑149 (MQ=255)
                                                                                         aGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGt                                                     <  8:188596/149‑1 (MQ=255)
                                                                                         aGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGGTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGt                                                     >  5:56483/1‑149 (MQ=255)
                                                                                           aCATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTTAAATTGTGGGTTTTCGCGGTATCAACCGttt                                                   <  7:381681/149‑1 (MQ=255)
                                                                                           aCATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt                                                   <  5:370317/149‑1 (MQ=255)
                                                                                           aCATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt                                                   <  4:132700/149‑1 (MQ=255)
                                                                                               ttGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCg                                               >  3:101530/1‑149 (MQ=255)
                                                                                                tGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCg                                               <  4:101530/148‑1 (MQ=255)
                                                                                                   cATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                                                         <  6:45329/119‑1 (MQ=255)
                                                                                                       cTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCt                                       >  2:376177/1‑149 (MQ=255)
                                                                                                       cTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCt                                       >  3:141365/1‑149 (MQ=255)
                                                                                                        ttAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGt                                                                                                                         >  5:150486/1‑66 (MQ=255)
                                                                                                        ttAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGt                                                                                                                         <  6:150486/66‑1 (MQ=255)
                                                                                                        ttAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                      >  6:234360/1‑149 (MQ=255)
                                                                                                        ttAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                      >  3:246787/1‑149 (MQ=255)
                                                                                                        ttAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                      >  6:40940/1‑149 (MQ=255)
                                                                                                                 aTGTCATCACGTATTAGCAAGGCCTGTCCCGTTGTACTGCCAGCGTAAAGGATAAGTCACTGATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAACTTGTGGGTGTTCGCGGTATCAACCGTTGGTCGTTGATGCTGGaacaatac                             >  2:301266/1‑149 (MQ=255)
                                                                                                                  tGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaacaaaaca                            <  1:153189/149‑1 (MQ=255)
                                                                                                                      atcaCGTATTAGCAAGGCCTTTCCCTTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACg                         <  3:29232/148‑1 (MQ=255)
                                                                                                                       tcaCGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTGTCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTc                       >  7:273284/1‑149 (MQ=255)
                                                                                                                              ttAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTg                <  4:141365/149‑1 (MQ=255)
                                                                                                                                 gCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTggg              <  1:301266/148‑1 (MQ=255)
                                                                                                                                 gCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTggg              <  2:19453/148‑1 (MQ=255)
                                                                                                                                 gCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATATCTGGAGGGGATATGATTCTTGAGCGCGGTGAAATTGTGGGTTTTCGCGGTATCACCCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTtgg              >  3:346781/1‑148 (MQ=255)
                                                                                                                                      gCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAAc        >  2:300757/1‑149 (MQ=255)
                                                                                                                                        cTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGAGGGTTTTCGCGGTATCAACCGTTAGTCGTTGATGATGGAACAAACCAACGTCCTGATTGGGGAGAAcgc      >  6:135071/1‑149 (MQ=255)
                                                                                                                                         tttCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGgag           >  7:5471/1‑143 (MQ=255)
                                                                                                                                         tttCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGgag           <  8:5471/143‑1 (MQ=255)
                                                                                                                                            cccGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACGCGTg   <  6:288915/148‑1 (MQ=255)
                                                                                                                                             ccGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACGCGTgg  <  5:40940/148‑1 (MQ=255)
                                                                                                                                                  ataCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCg                                               <  3:174630/98‑1 (MQ=255)
                                                                                                                                                  ataCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCg                                               >  4:174630/1‑98 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACGCGTGG  >  NC_000913/916283‑916571

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