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.8% 316.0 / NA 81R443R (CGT→CGCrnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/49);  minor base A (1/0);  total (32/49)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCG  >  NC_000913/4407835‑4408127
                                                                                                                                                   |                                                                                                                                                 
ggATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCAtt                                                                                                                                                <  1:221837/151‑1 (MQ=255)
      gcgcGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGaa                                                                                                                                             <  2:1167172/148‑1 (MQ=255)
        gcgAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTAt                                                                                                                                        <  2:1155504/151‑1 (MQ=255)
         cgAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATc                                                                                                                                       >  2:244966/1‑151 (MQ=255)
              aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAAc                                                                                                                                   >  1:1184728/1‑150 (MQ=255)
                 aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAc                                                                                                                               >  1:1120515/1‑151 (MQ=255)
                 aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAc                                                                                                                               >  1:18657/1‑151 (MQ=255)
                   gCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGAccc                                                                                                                             <  1:1116817/151‑1 (MQ=255)
                    ccccGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCa                                                                                                                            <  1:413140/151‑1 (MQ=255)
                        gCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCg                                                                                                                         >  2:432951/1‑150 (MQ=255)
                           ggTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTaa                                                                                                                      <  2:536392/150‑1 (MQ=255)
                           ggTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAac                                                                                                                     >  2:993222/1‑151 (MQ=255)
                            gTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAacg                                                                                                                    <  2:862276/151‑1 (MQ=255)
                              tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAacga                                                                                                                   <  2:423723/150‑1 (MQ=255)
                                aGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGAATTGAACGTATCGAACAGACCCAGCGTAacgacg                                                                                                                 >  1:797618/1‑150 (MQ=255)
                                 gCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGAcgcg                                                                                                               <  2:589796/151‑1 (MQ=255)
                                  cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGCTCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACgcgc                                                                                                              >  2:382201/1‑151 (MQ=255)
                                   aTCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACgcgc                                                                                                              >  1:354688/1‑150 (MQ=255)
                                       cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcaca                                                                                                           <  2:1147529/149‑1 (MQ=255)
                                       cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACaaa                                                                                                         <  1:103089/151‑1 (MQ=255)
                                       cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACaa                                                                                                          >  1:896278/1‑150 (MQ=255)
                                          agagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaatta                                                                                                      <  2:631146/151‑1 (MQ=255)
                                             gTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaatt                                                                                                   <  2:1164397/151‑1 (MQ=255)
                                              ttGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                                                                  >  2:1126106/1‑151 (MQ=255)
                                                gCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                                                                  >  1:245695/1‑149 (MQ=255)
                                                 cATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGaa                                                                                                <  2:806635/150‑1 (MQ=255)
                                                           ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCAt                                                                                        <  1:9393/148‑1 (MQ=255)
                                                           ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGAt                                                                                     >  1:91388/1‑151 (MQ=255)
                                                           ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGAt                                                                                     >  1:625227/1‑151 (MQ=255)
                                                           ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGAt                                                                                     >  2:1023381/1‑151 (MQ=255)
                                                            taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGAtt                                                                                    >  1:832722/1‑151 (MQ=255)
                                                               aGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTg                                                                                 <  2:1107510/151‑1 (MQ=255)
                                                                  gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                                                               <  2:372318/150‑1 (MQ=255)
                                                                  gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                                                                              <  1:2349/151‑1 (MQ=255)
                                                                    tGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGa                                                                             >  2:983990/1‑150 (MQ=255)
                                                                     ggATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGaa                                                                            <  2:778730/150‑1 (MQ=255)
                                                                      gATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAata                                                                          <  2:821892/151‑1 (MQ=255)
                                                                       aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAatat                                                                         <  1:721600/151‑1 (MQ=255)
                                                                       aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAatat                                                                         <  2:1052457/151‑1 (MQ=255)
                                                                         aaaGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATAtct                                                                       <  1:215613/151‑1 (MQ=255)
                                                                            gCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcgg                                                                    >  1:324982/1‑151 (MQ=255)
                                                                            gCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcgg                                                                    <  2:44818/151‑1 (MQ=255)
                                                                              ccGTGGAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcg                                                                  >  1:528965/1‑151 (MQ=255)
                                                                              ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcg                                                                  >  2:447084/1‑151 (MQ=255)
                                                                                gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcgg                                                                 <  2:876382/150‑1 (MQ=255)
                                                                                gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCggcggc                                                                <  1:1181251/151‑1 (MQ=255)
                                                                                     gaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGt                                                            <  1:101342/150‑1 (MQ=255)
                                                                                       aCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTc                                                         >  2:443181/1‑151 (MQ=255)
                                                                                       aCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTc                                                         <  2:28112/151‑1 (MQ=255)
                                                                                        cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGttt                                                          <  2:229262/149‑1 (MQ=255)
                                                                                        cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCGCGGCGGCGCGTTTCg                                                        <  2:985507/151‑1 (MQ=255)
                                                                                          cgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCTCAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGtt                                                      <  1:944921/151‑1 (MQ=255)
                                                                                            gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTga                                                    <  1:559398/151‑1 (MQ=255)
                                                                                            gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTga                                                    <  1:970563/151‑1 (MQ=255)
                                                                                            gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTga                                                    <  2:449639/151‑1 (MQ=255)
                                                                                                   ctcATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCg                                             >  1:571270/1‑151 (MQ=255)
                                                                                                    tcATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCg                                             >  1:823320/1‑150 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaaa                                                <  2:477970/146‑1 (MQ=255)
                                                                                                     cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGaa                                           <  1:606146/151‑1 (MQ=255)
                                                                                                          gagaGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAAc                                      <  1:661304/151‑1 (MQ=255)
                                                                                                             agCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCgg                                   >  1:930919/1‑151 (MQ=255)
                                                                                                              gCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGTCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGc                                  <  2:468397/151‑1 (MQ=255)
                                                                                                               cGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCa                                 <  2:308523/151‑1 (MQ=255)
                                                                                                                 aagaagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACt                               <  2:1120538/151‑1 (MQ=255)
                                                                                                                 aagaagCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACt                               >  2:865478/1‑151 (MQ=255)
                                                                                                                     agCGAAGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTc                           <  2:608324/151‑1 (MQ=255)
                                                                                                                         aaGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGt                        >  2:619974/1‑150 (MQ=255)
                                                                                                                         aaGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTa                       <  1:1052275/151‑1 (MQ=255)
                                                                                                                          aGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTAt                      >  1:4214/1‑151 (MQ=255)
                                                                                                                          aGTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTAt                      >  1:277591/1‑151 (MQ=255)
                                                                                                                           gTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTAtt                     >  2:659908/1‑151 (MQ=255)
                                                                                                                           gTTCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTAtt                     <  1:783995/151‑1 (MQ=255)
                                                                                                                            ttCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTc                    >  1:983153/1‑151 (MQ=255)
                                                                                                                             tCATTTTCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTc                    >  1:700606/1‑150 (MQ=255)
                                                                                                                                   tCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACa              <  1:980092/150‑1 (MQ=255)
                                                                                                                                   tCAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACa              <  2:296597/150‑1 (MQ=255)
                                                                                                                                    cAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACaa             <  2:814609/150‑1 (MQ=255)
                                                                                                                                    cAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACaa             <  1:408905/150‑1 (MQ=255)
                                                                                                                                    cAACGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACaa             <  2:687441/150‑1 (MQ=255)
                                                                                                                                       cGCTGAACGCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCg         >  2:845015/1‑151 (MQ=255)
                                                                                                                                               gCCGCATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCg  <  1:983992/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCG  >  NC_000913/4407835‑4408127

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