Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,728,608 G→A 5.2% G88S (GGC→AGC)  rnt → RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,728,6080GA 5.2% 252.6 / 2.2 76G88S (GGC→AGC) rntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base G (43/29);  new base A (2/2);  total (45/31)
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.77e-01

TTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCT  >  NC_000913/1728462‑1728755
                                                                                                                                                  |                                                                                                                                                   
ttGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATcgcggc                                                                                                                                                   <  1:83080/149‑1 (MQ=255)
      ccATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTc                                                                                                                                             >  1:734622/1‑149 (MQ=255)
       cATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCa                                                                                                                                            >  2:87459/1‑149 (MQ=255)
        aTCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCCCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAg                                                                                                                                           >  4:576846/1‑149 (MQ=255)
                 aaaaTGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTCGCAACCAGAGGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGACTACGa                                                                                                                                  <  4:508387/149‑1 (MQ=255)
                    aTGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGgc                                                                                                                               >  2:327688/1‑149 (MQ=255)
                           aaCAAGGCTGGCTGATGCCGGACACCACATTACATTTCCAGGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCAc                                                                                                                        >  4:421775/1‑149 (MQ=255)
                           aaCAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCAc                                                                                                                        >  2:102635/1‑149 (MQ=255)
                           aaCAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCAc                                                                                                                        >  3:66641/1‑149 (MQ=255)
                           aaCAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCAc                                                                                                                        >  4:607621/1‑149 (MQ=255)
                             cAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGa                                                                                                                      >  3:535559/1‑149 (MQ=255)
                                    ggctgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAAttttt                                                                                                               >  3:353787/1‑149 (MQ=255)
                                    ggctgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAAttttt                                                                                                               >  2:781241/1‑149 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaaa                                                                                                            <  1:528752/149‑1 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaaa                                                                                                            <  1:500067/149‑1 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaaa                                                                                                            >  2:617044/1‑149 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaaa                                                                                                            >  4:600974/1‑149 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCAGCGCGGTCAGCGAATACGAgg                                                                                                                                <  2:723640/129‑1 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCAGCGCGGTCAGCGAATACGAgg                                                                                                                                >  3:67819/1‑129 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCAGCGCGGTCAGCGAATACGAgg                                                                                                                                <  4:67819/129‑1 (MQ=255)
                                       tgATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCAGCGCGGTCAGCGAATACGAgg                                                                                                                                >  1:723640/1‑129 (MQ=255)
                                        gATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAg                                                                                                           >  4:367134/1‑149 (MQ=255)
                                              ggACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTAc                                                                                                     <  4:412352/149‑1 (MQ=255)
                                                acacCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGt                                                                                                   >  3:214293/1‑149 (MQ=255)
                                                    cacaTTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaa                                                                                                             <  4:296885/135‑1 (MQ=255)
                                                    cacaTTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTaa                                                                                                             >  3:296885/1‑135 (MQ=255)
                                                    cacaTTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAg                                                                                               <  2:734622/149‑1 (MQ=255)
                                                    cacaTTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAg                                                                                               >  1:181681/1‑149 (MQ=255)
                                                        ttacattTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTAt                                                                                           <  3:600974/149‑1 (MQ=255)
                                                        ttacattTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTAt                                                                                           <  1:617044/149‑1 (MQ=255)
                                                          acattTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTa                                                                                         >  1:255072/1‑149 (MQ=255)
                                                           cattTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTaa                                                                                        <  1:87459/149‑1 (MQ=255)
                                                              ttCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAagc                                                                                     >  4:107783/1‑149 (MQ=255)
                                                                 cACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTaa                                                                                                 <  2:591307/134‑1 (MQ=255)
                                                                 cACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTaa                                                                                                 >  1:591307/1‑134 (MQ=255)
                                                                 cACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTaaa                                                                                       >  3:802077/1‑144 (MQ=255)
                                                                 cACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTaaa                                                                                       <  4:802077/144‑1 (MQ=255)
                                                                       gAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTg                                                                            >  4:624322/1‑149 (MQ=255)
                                                                       gAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTg                                                                            >  3:577111/1‑149 (MQ=255)
                                                                        aaCCATTTGTCGGCGCAAATTTGCAACCAGAAGCCTTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGt                                                                           >  1:715455/1‑149 (MQ=255)
                                                                         aCCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTa                                                                          >  1:755614/1‑149 (MQ=255)
                                                                           cATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAAc                                                                        <  1:16523/149‑1 (MQ=255)
                                                                             tttGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACcg                                                                      >  1:487042/1‑149 (MQ=255)
                                                                                     cgcAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTa                                                              <  3:672590/149‑1 (MQ=255)
                                                                                         aaTTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCt                                                                             >  4:614481/1‑130 (MQ=255)
                                                                                         aaTTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCt                                                                             <  3:614481/130‑1 (MQ=255)
                                                                                           tttGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTa                                                                                                              >  3:118916/1‑95 (MQ=255)
                                                                                           tttGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTa                                                                                                              <  4:118916/95‑1 (MQ=255)
                                                                                                aaCCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGcac                                                   <  4:535559/149‑1 (MQ=255)
                                                                                                  ccAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACaa                                                 >  1:405202/1‑149 (MQ=255)
                                                                                                  ccAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACaa                                                 >  1:354169/1‑149 (MQ=255)
                                                                                                     gAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAagc                                                                                     <  1:4167/110‑1 (MQ=255)
                                                                                                     gAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAagc                                                                                     >  2:4167/1‑110 (MQ=255)
                                                                                                        gCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCaa                                           <  3:576846/149‑1 (MQ=255)
                                                                                                            tCGCCTTCANCGGCATNGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAAtttt                                       <  2:282661/149‑1 (MQ=255)
                                                                                                            tCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAAtttt                                       >  3:16426/1‑149 (MQ=255)
                                                                                                             cGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTg                                      >  2:377988/1‑149 (MQ=255)
                                                                                                             cGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTg                                      >  4:570993/1‑149 (MQ=255)
                                                                                                              gCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGa                                     <  3:367134/149‑1 (MQ=255)
                                                                                                               ccTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGAt                                    <  4:214293/149‑1 (MQ=255)
                                                                                                                  tCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATcac                                 >  2:532029/1‑149 (MQ=255)
                                                                                                                  tCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATcac                                 >  2:786230/1‑149 (MQ=255)
                                                                                                                     aCGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGGAATTTTTa                                                                                                              <  3:624745/69‑1 (MQ=255)
                                                                                                                     aCGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGGAATTTTTa                                                                                                              >  4:624745/1‑69 (MQ=255)
                                                                                                                     aCGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGc                              <  1:142586/149‑1 (MQ=255)
                                                                                                                      cGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCt                             <  3:397009/149‑1 (MQ=255)
                                                                                                                        gCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCttt                           >  3:448311/1‑149 (MQ=255)
                                                                                                                        gCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCttt                           <  4:577111/149‑1 (MQ=255)
                                                                                                                         cATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTa                          <  2:755614/149‑1 (MQ=255)
                                                                                                                               cccGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATgg                    >  1:605175/1‑149 (MQ=255)
                                                                                                                                ccGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGgc                   >  2:537193/1‑149 (MQ=255)
                                                                                                                                      gATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCaa                                           <  3:225689/119‑1 (MQ=255)
                                                                                                                                      gATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCaa                                           >  4:225689/1‑119 (MQ=255)
                                                                                                                                       aTCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTat                         <  2:735864/136‑1 (MQ=255)
                                                                                                                                       aTCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTat                         >  1:735864/1‑136 (MQ=255)
                                                                                                                                        tCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTCAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGc                              >  4:129666/1‑130 (MQ=255)
                                                                                                                                        tCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGc                              <  3:129666/130‑1 (MQ=255)
                                                                                                                                                 cggcgcggTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCt  <  3:607621/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCT  >  NC_000913/1728462‑1728755

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