Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,407,982 T→C R443R (CGT→CGC rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,9820TC98.9% 356.3 / ‑3.8 94R443R (CGT→CGCrnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (47/46);  total (47/47)
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.13e-01

GGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAA  >  NC_000913/4407835‑4408129
                                                                                                                                                   |                                                                                                                                                   
ggATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCAtt                                                                                                                                                  >  1:457045/1‑151 (MQ=255)
     tGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACg                                                                                                                                             >  2:695513/1‑151 (MQ=255)
      gcgcGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTg                                                                                                                                                 <  1:139903/146‑1 (MQ=255)
          gAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCg                                                                                                                                        <  3:1093945/151‑1 (MQ=255)
          gAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCg                                                                                                                                        >  4:171723/1‑151 (MQ=255)
               gTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAAc                                                                                                                                     >  2:781353/1‑149 (MQ=255)
                  cGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAcc                                                                                                                                >  4:298464/1‑151 (MQ=255)
                   gCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAccc                                                                                                                               <  2:339439/151‑1 (MQ=255)
                   gCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAccc                                                                                                                               >  3:943643/1‑151 (MQ=255)
                    ccccGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCa                                                                                                                              <  3:383034/151‑1 (MQ=255)
                     cccGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAg                                                                                                                             <  1:161220/151‑1 (MQ=255)
                        gCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGt                                                                                                                          <  3:1087334/151‑1 (MQ=255)
                           ggTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTa                                                                                                                         <  4:334555/149‑1 (MQ=255)
                               aaGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAacgac                                                                                                                    <  3:1303438/150‑1 (MQ=255)
                                 gCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGAcgcg                                                                                                                 <  3:818555/151‑1 (MQ=255)
                                  cATCTCTAAGAGTTGCATATCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGAGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTATCAACGCTCAACGCCGCATTGAACGTATCGAACAGTCCCAGCGTAACGAcgcg                                                                                                                 <  3:1244858/150‑1 (MQ=255)
                                       cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGca                                                                                                               <  3:517609/147‑1 (MQ=255)
                                       cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGca                                                                                                               <  1:237411/147‑1 (MQ=255)
                                        gAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaat                                                                                                          >  4:233812/1‑151 (MQ=255)
                                          agagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaatt                                                                                                         >  1:242374/1‑150 (MQ=255)
                                           gagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaatta                                                                                                        <  3:708835/150‑1 (MQ=255)
                                           gagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaatt                                                                                                         >  4:26586/1‑149 (MQ=255)
                                               tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                                                                    >  3:243193/1‑150 (MQ=255)
                                                gCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGaa                                                                                                  <  3:924594/151‑1 (MQ=255)
                                                  aTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGa                                                                                                   >  3:832772/1‑148 (MQ=255)
                                                  aTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGa                                                                                                   >  3:338775/1‑148 (MQ=255)
                                                   tAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAaga                                                                                                <  3:849367/150‑1 (MQ=255)
                                                    aaCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGt                                                                                              >  1:238958/1‑151 (MQ=255)
                                                      ccTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGc                                                                                            <  2:353036/151‑1 (MQ=255)
                                                      ccTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGc                                                                                            <  4:4794/151‑1 (MQ=255)
                                                       ctctATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGc                                                                                            >  2:164314/1‑150 (MQ=255)
                                                             aaaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTc                                                                                     >  2:26443/1‑151 (MQ=255)
                                                              aaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCt                                                                                    <  2:608702/151‑1 (MQ=255)
                                                                gTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTgg                                                                                  >  3:1308275/1‑151 (MQ=255)
                                                                 tGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                                                                 >  2:1120528/1‑151 (MQ=255)
                                                                  gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                                                                                <  2:1402892/151‑1 (MQ=255)
                                                                    tGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGa                                                                               >  2:209293/1‑150 (MQ=255)
                                                                    tGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGa                                                                               >  4:175365/1‑150 (MQ=255)
                                                                       aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGaa                                                                              >  4:186991/1‑148 (MQ=255)
                                                                       aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGaa                                                                              <  1:186993/148‑1 (MQ=255)
                                                                        tAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGACTCTGGCGAATAtc                                                                          >  2:1311520/1‑151 (MQ=255)
                                                                         aaaGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATAtct                                                                         <  3:200843/151‑1 (MQ=255)
                                                                         aaaGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATAtc                                                                          <  2:332996/150‑1 (MQ=255)
                                                                          aaGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATAtctc                                                                        <  1:352921/151‑1 (MQ=255)
                                                                           aGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcg                                                                       <  4:474191/151‑1 (MQ=255)
                                                                             cccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggc                                                                     >  1:118814/1‑151 (MQ=255)
                                                                             cccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggc                                                                     <  3:1071223/151‑1 (MQ=255)
                                                                                gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCggcggc                                                                  >  4:305349/1‑151 (MQ=255)
                                                                                gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCggcggc                                                                  <  1:85327/151‑1 (MQ=255)
                                                                                gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCggcggc                                                                  >  2:660173/1‑151 (MQ=255)
                                                                                 tGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGgcggcg                                                                 >  2:854230/1‑151 (MQ=255)
                                                                                  gaagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGgcggcg                                                                 <  3:97058/150‑1 (MQ=255)
                                                                                   aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGcgcg                                                               >  3:265990/1‑151 (MQ=255)
                                                                                      aaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGttt                                                            >  3:1337145/1‑151 (MQ=255)
                                                                                      aaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGttt                                                            >  2:53515/1‑151 (MQ=255)
                                                                                      aaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGttt                                                            >  2:475757/1‑151 (MQ=255)
                                                                                       aCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTc                                                           <  3:928296/151‑1 (MQ=255)
                                                                                        cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCg                                                          >  3:1067164/1‑151 (MQ=255)
                                                                                          cgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGtt                                                        >  1:298716/1‑151 (MQ=255)
                                                                                              gggATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTgag                                                     <  2:1044379/150‑1 (MQ=255)
                                                                                                gATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaa                                                   >  2:277691/1‑150 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaaa                                                  <  3:221604/146‑1 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAg                                                 <  1:233814/147‑1 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAg                                                 <  3:728746/147‑1 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGc                                                <  2:124238/148‑1 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGa                                              <  4:92203/150‑1 (MQ=255)
                                                                                                      aTTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGaaa                                            <  2:474125/151‑1 (MQ=255)
                                                                                                         tGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGa                                          >  2:21110/1‑150 (MQ=255)
                                                                                                         tGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGa                                          >  4:343972/1‑150 (MQ=255)
                                                                                                          gagaGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAAc                                        >  2:1166235/1‑151 (MQ=255)
                                                                                                            gagCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCg                                      <  4:42493/151‑1 (MQ=255)
                                                                                                             agCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCgg                                     >  2:941750/1‑151 (MQ=255)
                                                                                                              gCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGc                                    >  1:294365/1‑151 (MQ=255)
                                                                                                                   gaagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTg                                >  3:809625/1‑150 (MQ=255)
                                                                                                                   gaagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGt                               <  2:943606/151‑1 (MQ=255)
                                                                                                                   gaagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGt                               <  3:695537/151‑1 (MQ=255)
                                                                                                                    aagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGtt                              <  3:26449/151‑1 (MQ=255)
                                                                                                                     agCGAAGTTCATTTTCCACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAAATAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGtt                              >  4:198720/1‑150 (MQ=255)
                                                                                                                     agCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTc                             >  3:188345/1‑151 (MQ=255)
                                                                                                                      gCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTcc                            >  1:407541/1‑151 (MQ=255)
                                                                                                                         aaGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGt                          >  1:103830/1‑150 (MQ=255)
                                                                                                                         aaGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTa                         <  2:109285/151‑1 (MQ=255)
                                                                                                                          aGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTa                         <  3:164337/150‑1 (MQ=255)
                                                                                                                             tCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTc                      <  3:138034/150‑1 (MQ=255)
                                                                                                                              cATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCa                     <  4:78329/150‑1 (MQ=255)
                                                                                                                               aTTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACg                   >  1:371690/1‑151 (MQ=255)
                                                                                                                               aTTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACg                   <  4:241348/151‑1 (MQ=255)
                                                                                                                                ttttCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACg                   >  1:156798/1‑150 (MQ=255)
                                                                                                                                    cAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACaa               <  4:236490/150‑1 (MQ=255)
                                                                                                                                      aCGCTGAACGCCGCATTGAAAGTATCGAACAGAACCAGCGAAAAGACGCGCACAAATTAATTGGCGAGTGCATGATTCTGGCGAATATCTCCGCCGCGCGCTTCGTTGAGACCACGAAAGAACCGGCACTGTTCCGCATTCACGACGAGcc            >  4:393749/1‑151 (MQ=255)
                                                                                                                                          tGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAg         >  4:39569/1‑150 (MQ=255)
                                                                                                                                            aaCGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCAc      <  4:145229/151‑1 (MQ=255)
                                                                                                                                            aaCGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCAc      >  2:118720/1‑151 (MQ=255)
                                                                                                                                              cgccgcATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCg    >  2:167441/1‑151 (MQ=255)
                                                                                                                                                ccGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGaa  <  4:148627/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAA  >  NC_000913/4407835‑4408129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: