Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,213 A→T 9.8% intergenic (‑94/+45) glnP ← / ← glnH glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,2130AT 9.8% 327.8 / 23.2 104intergenic (‑94/+45)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base T (3/7);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 105 total aligned reads displayed.

CCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCG  >  NC_000913/847067‑847360
                                                                                                                                                  |                                                                                                                                                    
ccTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAt                                                                                                                                                     >  1:787385/1‑148 (MQ=255)
ccTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAt                                                                                                                                                     >  1:760702/1‑148 (MQ=255)
 cTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                    <  1:387879/148‑1 (MQ=255)
  ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                 >  4:673288/1‑149 (MQ=255)
  ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                 >  3:739538/1‑149 (MQ=255)
  ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                 >  4:352985/1‑149 (MQ=255)
      aTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                             >  4:806408/1‑149 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                            <  4:791648/149‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                             <  2:522987/148‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                             <  1:776790/148‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                             <  1:243993/148‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACCGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                            <  4:739538/149‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGagag                                                                                                                                           <  4:56580/149‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGagag                                                                                                                                           <  4:298961/149‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGAGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGgaga                                                                                                                                            <  1:462673/148‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTCGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGagag                                                                                                                                           <  4:194816/149‑1 (MQ=255)
            aGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGagag                                                                                                                                           <  2:701634/145‑1 (MQ=255)
            aGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGagag                                                                                                                                           >  1:701634/1‑145 (MQ=255)
             gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAAc                                                                                                                                      <  3:352985/149‑1 (MQ=255)
              cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCCTAT‑GGGGAGAggattag                                                                                                                                     <  3:739655/149‑6 (MQ=255)
               ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGAGAggattag                                                                                                                                     <  2:854/148‑6 (MQ=255)
               ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCC‑TATGGGGAGAggattagg                                                                                                                                    <  3:110593/149‑7 (MQ=255)
                gAATGGCAGGCCAGATGGCACTCCAGTCAGACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGc                                                                                                                                   >  3:677918/1‑149 (MQ=255)
                gAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGc                                                                                                                                   <  3:805261/149‑1 (MQ=255)
                        ggCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCg                                                                                                                            >  1:352630/1‑148 (MQ=255)
                         gCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                          <  4:780860/149‑1 (MQ=255)
                         gCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGAGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                          >  4:63810/1‑149 (MQ=255)
                              aTGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                      <  2:896093/148‑1 (MQ=255)
                                    cTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGgtgt                                                                                                                >  1:449280/1‑148 (MQ=255)
                                     tCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTa                                                                                                               <  1:897429/148‑1 (MQ=255)
                                         gTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAAGCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtt                                                                                                          <  3:63810/149‑1 (MQ=255)
                                              aCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                      >  2:860024/1‑148 (MQ=255)
                                               cTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTtattat                                                                                                     <  1:300243/148‑1 (MQ=255)
                                                 gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                   <  2:787385/148‑1 (MQ=255)
                                                   atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                >  3:183782/1‑149 (MQ=255)
                                                     atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                  >  1:67265/1‑145 (MQ=255)
                                                     atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                  <  2:67265/145‑1 (MQ=255)
                                                     atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGtt                                                                                              <  3:174887/149‑1 (MQ=255)
                                                     atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGtt                                                                                              <  4:323489/149‑1 (MQ=255)
                                                       gttgttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGAGAgga                                                                                                                                         >  4:413902/1‑103 (MQ=255)
                                                       gttgttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                           >  4:430755/1‑134 (MQ=255)
                                                       gttgttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCC‑TATGGGGAGAgga                                                                                                                                         <  3:413902/104‑2 (MQ=255)
                                                              cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCg                                                                                      <  4:703353/148‑1 (MQ=255)
                                                              cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCg                                                                                      >  3:703353/1‑148 (MQ=255)
                                                               tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                    >  3:771129/1‑149 (MQ=255)
                                                                 ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAc                                                                                   >  2:335753/1‑148 (MQ=255)
                                                                 ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAc                                                                                   <  2:40812/148‑1 (MQ=255)
                                                                       tggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAttttt                                                                            <  3:687908/149‑1 (MQ=255)
                                                                        ggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                           <  3:492701/149‑1 (MQ=255)
                                                                         gtgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                           <  2:708646/148‑1 (MQ=255)
                                                                          tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                          >  1:50825/1‑148 (MQ=255)
                                                                           gtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTa                                                                        <  4:183782/149‑1 (MQ=255)
                                                                              aataatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAt                                                                     <  4:771129/149‑1 (MQ=255)
                                                                                 aatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAttt                                                                   >  1:127608/1‑148 (MQ=255)
                                                                                  atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTc                                                                  <  2:50825/148‑1 (MQ=255)
                                                                                   tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                 >  2:834815/1‑148 (MQ=255)
                                                                                   tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                 >  1:139099/1‑148 (MQ=255)
                                                                                   tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                >  3:206852/1‑149 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                             <  2:712410/71‑1 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                <  2:503786/148‑1 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                >  1:788862/1‑148 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                <  2:180567/148‑1 (MQ=255)
                                                                                     aaaTTCAGGGGAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                              >  3:57166/1‑149 (MQ=255)
                                                                                        ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                            <  2:352630/148‑1 (MQ=255)
                                                                                        ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                            <  1:170644/148‑1 (MQ=255)
                                                                                             ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                <  4:583807/139‑1 (MQ=255)
                                                                                             ggTAATTAATGGCCCCGGAGCGGGGTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCTGTGCGTTTTACAGGTGTAGCGTTATGATTTCGGTTCAGTACCGAACGATTTTTTGTAGATTTCGTAGTAAGTTct                                                      >  3:764888/1‑148 (MQ=255)
                                                                                             ggTAATTAATGACCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                >  3:583807/1‑139 (MQ=255)
                                                                                                aattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTattatt                                                                                                    >  1:383232/1‑100 (MQ=255)
                                                                                                aattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTattatt                                                                                                    <  2:383232/100‑1 (MQ=255)
                                                                                                  ttaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtc                                                  >  2:330696/1‑148 (MQ=255)
                                                                                                   taatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtct                                                 >  2:542365/1‑148 (MQ=255)
                                                                                                       ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCgcgc                                            >  4:499777/1‑149 (MQ=255)
                                                                                                           ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAggg                                        >  4:46463/1‑149 (MQ=255)
                                                                                                           ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                         >  1:493015/1‑148 (MQ=255)
                                                                                                           ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                         <  2:760702/148‑1 (MQ=255)
                                                                                                               aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGttt                                     <  2:773148/148‑1 (MQ=255)
                                                                                                               aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGttt                                     >  1:196659/1‑148 (MQ=255)
                                                                                                                   ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                <  4:734569/149‑1 (MQ=255)
                                                                                                                            cccTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGtt                        >  2:854569/1‑148 (MQ=255)
                                                                                                                             ccTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                       <  1:516383/148‑1 (MQ=255)
                                                                                                                             ccTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                      <  4:804868/149‑1 (MQ=255)
                                                                                                                             ccTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                      <  4:206852/149‑1 (MQ=255)
                                                                                                                                cACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                    <  1:834815/148‑1 (MQ=255)
                                                                                                                                 aCCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACtt                   <  1:335753/148‑1 (MQ=255)
                                                                                                                                    gatctCCTCTCCCCATA‑GGGGAGAGATAAGTGCCAGTGCGTTTTACAGGTGTAGCGTTAGTATATCGGTACAGTGCCGAAACATTTTTTGTAGATTTCGTTGTAAGTTCCGTTATCGCTGAGGGATTGCAACGTGCCGTTGACTTTGTc               <  4:559968/145‑1 (MQ=255)
                                                                                                                                    cTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                >  1:117692/1‑148 (MQ=255)
                                                                                                                                    cTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                >  2:178238/1‑148 (MQ=255)
                                                                                                                                        tCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                    <  3:358881/76‑1 (MQ=255)
                                                                                                                                        tCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACg            <  1:860024/148‑1 (MQ=255)
                                                                                                                                         ccTCTCCCCTAT‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc           >  2:480484/1‑148 (MQ=255)
                                                                                                                                          ctctCCCCATA‑GGGGAGAGGGAACTGACAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAACATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTCTCaagatc         >  3:559968/1‑144 (MQ=255)
                                                                                                                                           tctcCCCTAT‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAggg                                        <  2:885094/117‑1 (MQ=255)
                                                                                                                                           tctcCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg         >  1:120830/1‑148 (MQ=255)
                                                                                                                                             tcCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtc      >  4:679237/1‑149 (MQ=255)
                                                                                                                                             tcCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTCCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTATCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtc      >  4:593360/1‑149 (MQ=255)
                                                                                                                                              ccccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtc      >  1:189869/1‑148 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTc                                   >  4:787219/1‑117 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTc                                   <  3:787219/117‑1 (MQ=255)
                                                                                                                                                  aTA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg  >  2:220778/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCG  >  NC_000913/847067‑847360

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