Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 355,892 G→A 5.5% W324* (TGG→TAG)  codB → cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913355,8920GA 5.5% 233.5 / 2.6 73W324* (TGG→TAG) codBcytosine transporter
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base A (2/2);  total (36/37)
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.03e-01

TATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGG  >  NC_000913/355746‑356035
                                                                                                                                                  |                                                                                                                                               
tatCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGt                                                                                                                                               <  2:277781/149‑1 (MQ=255)
        ccaccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttc                                                                                                                                       >  4:259661/1‑149 (MQ=255)
         caccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcc                                                                                                                                      >  5:321524/1‑149 (MQ=255)
         caccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcc                                                                                                                                      >  3:162425/1‑149 (MQ=255)
         caccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCGTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcc                                                                                                                                      >  5:223858/1‑149 (MQ=255)
         caccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCGTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcc                                                                                                                                      >  5:230314/1‑149 (MQ=255)
           ccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTgg                                                                                                                                                <  4:359493/137‑1 (MQ=255)
           ccaACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTgg                                                                                                                                                >  3:359493/1‑137 (MQ=255)
            caACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCttt                                                                                                                                   >  1:355296/1‑149 (MQ=255)
            caACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCttt                                                                                                                                   >  8:153886/1‑149 (MQ=255)
             aaCTATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTc                                                                                                                                  <  5:505/149‑1 (MQ=255)
                  tAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAg                                                                                                                             >  8:328984/1‑149 (MQ=255)
                        aCTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATtc                                                                                                                       >  8:50514/1‑149 (MQ=255)
                         ctctATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATtcc                                                                                                                      <  3:337269/149‑1 (MQ=255)
                          tctATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcct                                                                                                                                     >  8:359075/1‑133 (MQ=255)
                          tctATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGAccttcct                                                                                                                                     <  7:359075/133‑1 (MQ=255)
                              ttcgtatggTTTAGGTTTCGACAACATTACAGGGATGGCGAGCAAAACCCTTTCGGTAATCAACGTTATTATAGGTAATTTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTTACCTTCCGTTCGGCAGCTATTCCTCCa                                                                                                                  <  8:143158/141‑1 (MQ=255)
                              tGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCa                                                                                                                  <  5:2093/148‑1 (MQ=255)
                                 gTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgg                                                                                                              <  3:293358/149‑1 (MQ=255)
                                  tCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtggg                                                                                                             <  7:57729/149‑1 (MQ=255)
                                    gggTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgggtg                                                                                                           >  3:307516/1‑149 (MQ=255)
                                    gggTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgggtg                                                                                                           <  4:391824/149‑1 (MQ=255)
                                    gggTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgggt                                                                                                            <  5:400541/148‑1 (MQ=255)
                                      gTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCa                                                                                                                  >  4:413096/1‑140 (MQ=255)
                                       tttAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCg                                                                                                        >  3:40938/1‑149 (MQ=255)
                                           ggTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGat                                                                                                    <  4:307516/149‑1 (MQ=255)
                                            gTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGatc                                                                                                   <  6:200415/149‑1 (MQ=255)
                                                  ccAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGcc                                                                                             >  2:387737/1‑149 (MQ=255)
                                                   cAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGcc                                                                                             >  7:121978/1‑148 (MQ=255)
                                                    aaCATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGGCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGa                                                                                           <  7:328984/149‑1 (MQ=255)
                                                    aaCATTACCGGGATGTAGAGCAAAACCCTTTAGGTTATCAACGTTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGGCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGa                                                                                           <  6:223858/149‑1 (MQ=255)
                                                    aaCATTACCGGGAGGTCGAGCAAAACCCTTTCGGTAATCCACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGGCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGa                                                                                           <  6:230314/149‑1 (MQ=255)
                                                            cGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTAGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                   >  8:255595/1‑149 (MQ=255)
                                                            cGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTAGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                   >  8:262351/1‑149 (MQ=255)
                                                                      aGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTcg                                                                         >  3:168683/1‑149 (MQ=255)
                                                                         aaaaCCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTcgccg                                                                      <  7:50514/149‑1 (MQ=255)
                                                                         aaaaCCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTcgccg                                                                      >  7:366831/1‑149 (MQ=255)
                                                                         aaaaCCCTTTCGGTAATAAACAGTATTATCGGTACGGTCTGCGCATTATCTCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATACCTCCAGTGGGTGGCGTGATCAACGCCGACTTTCTGATGAACCGAcgacc                                                                      >  4:255722/1‑146 (MQ=255)
                                                                             cccTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTAt                                                                  <  6:321524/149‑1 (MQ=255)
                                                                             cccTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTAt                                                                  <  4:40938/149‑1 (MQ=255)
                                                                                tttCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAg                                                               >  3:165767/1‑149 (MQ=255)
                                                                                tttCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAg                                                               >  2:308166/1‑149 (MQ=255)
                                                                                  tCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAg                                                               <  1:387737/147‑1 (MQ=255)
                                                                                        aTCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATCACAATTTTGTCGGCTGGTGGACCTTCCTTTCGGCAGCTATTCCTCCATTTGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGc                                                       <  5:375963/149‑1 (MQ=255)
                                                                                        aTCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGc                                                       <  3:259661/149‑1 (MQ=255)
                                                                                           aaCGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGa                                                     <  1:308166/148‑1 (MQ=255)
                                                                                            aCGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGAc                                                    >  5:375676/1‑148 (MQ=255)
                                                                                                tattatCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCAcgc                                               >  8:74536/1‑149 (MQ=255)
                                                                                                tattatCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCAcgc                                               <  2:355296/149‑1 (MQ=255)
                                                                                                 attatCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCttt                                                                                                                                   >  1:281875/1‑64 (MQ=255)
                                                                                                 attatCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCttt                                                                                                                                   <  2:281875/64‑1 (MQ=255)
                                                                                                    atCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTa                                            <  3:167852/148‑1 (MQ=255)
                                                                                                     tCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTat                                           >  8:235118/1‑148 (MQ=255)
                                                                                                       ggTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTGCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatgat                                        <  3:255722/149‑1 (MQ=255)
                                                                                                       ggTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatga                                         <  4:162425/148‑1 (MQ=255)
                                                                                                        gTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCAcgc                                               >  5:112879/1‑141 (MQ=255)
                                                                                                        gTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCAcgc                                               <  6:112879/141‑1 (MQ=255)
                                                                                                        gTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatgat                                        >  4:123915/1‑148 (MQ=255)
                                                                                                         tACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgggt                                                                                                            >  6:100364/1‑79 (MQ=255)
                                                                                                         tACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAgtgggt                                                                                                            <  5:100364/79‑1 (MQ=255)
                                                                                                         tACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTAtgatga                                      <  7:153886/149‑1 (MQ=255)
                                                                                                         tACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTAtgatga                                      <  7:235118/149‑1 (MQ=255)
                                                                                                             gTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAgtgt                                  <  8:366831/149‑1 (MQ=255)
                                                                                                                   gcATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTAATCCAGTTGGTTGCGTGATCATCGCCGACTATCTCATGAACAGTCGCCGCTATGATCAATTTGCGACCAAGCGTATGATGAGTGTCAAtag                            >  3:257839/1‑147 (MQ=255)
                                                                                                                      ttATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGt                         >  8:417247/1‑149 (MQ=255)
                                                                                                                                 taACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTg               <  3:123915/148‑1 (MQ=255)
                                                                                                                                  aaCAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGc             >  4:96373/1‑149 (MQ=255)
                                                                                                                                  aaCAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGc             >  6:62244/1‑149 (MQ=255)
                                                                                                                                     aaTTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGt          >  2:246383/1‑149 (MQ=255)
                                                                                                                                     aaTTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCGGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGt          >  4:298761/1‑149 (MQ=255)
                                                                                                                                            tCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTg   <  8:121978/149‑1 (MQ=255)
                                                                                                                                             cGGCTAGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgg  <  7:255595/149‑1 (MQ=255)
                                                                                                                                              ggCTAGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgg  <  7:262351/148‑1 (MQ=255)
                                                                                                                                               gCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGcc      >  3:407977/1‑143 (MQ=255)
                                                                                                                                               gCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGcc      <  4:407977/143‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGG  >  NC_000913/355746‑356035

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: ‑