Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,473,614 C→T 12.4% W156* (TGG→TGA rpsG ← 30S ribosomal subunit protein S7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,473,6140CT12.4% 149.0 / 9.1 63W156* (TGG→TGArpsG30S ribosomal subunit protein S7
Reads supporting (aligned to +/- strand):  ref base C (22/32);  new base T (3/5);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGTTTAC  >  NC_000913/3473473‑3473756
                                                                                                                                             |                                                                                                                                              
aCCGACGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAAcc                                                                                                                                               <  1:536895/143‑1 (MQ=255)
 ccGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCACTTTAAGTAGCCCAAAGCGAGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGACAACCa                                                                                                                                              >  2:572314/1‑143 (MQ=255)
 ccGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAAccg                                                                                                                                              >  2:443044/1‑142 (MQ=255)
  cGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCaa                                                                                                                                             <  2:278749/143‑1 (MQ=255)
  cGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGGTTCCTCGTTTATCTTTTAGGCGTGCCGTTTAGGTAGCTGAAAGCGGGCAGCTTTGTGGAAGCGCCCGCCGGGGGAAGAAAACTCCGAAGGGATAACCaa                                                                                                                                             >  1:821695/1‑143 (MQ=255)
  cGATGTTACGGTAGCGTGCGATGGGTGGTGTACGAGGCATGTGTTGCCGCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGGTAACCaa                                                                                                                                             >  1:494802/1‑143 (MQ=255)
  cGATGTTACGGTAGCGTGCGATGGGTGGTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCaa                                                                                                                                             <  2:220521/143‑1 (MQ=255)
         aCGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAACTAGCCCATCGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAgt                                                                                                                                      <  2:794274/143‑1 (MQ=255)
           ggTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAgtgt                                                                                                                                    <  1:38266/143‑1 (MQ=255)
                cGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGaa                                                                                                                               <  1:560052/143‑1 (MQ=255)
                 gTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAAc                                                                                                                              <  1:724589/143‑1 (MQ=255)
                         ggTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTGGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCttgtt                                                                                                                      <  1:586001/143‑1 (MQ=255)
                               gacccaGCGATTTGTTGCCTCGTTTATCTTTTCGGCCTTCAATTTAAGGAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTCGTGATTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTc                                                                                                                <  1:441031/138‑1 (MQ=255)
                                 aCGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAGCTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCgg                                                                                                              <  2:193793/143‑1 (MQ=255)
                                 aCGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCgg                                                                                                              <  1:184022/143‑1 (MQ=255)
                                   gAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGcc                                                                                                            >  1:796791/1‑143 (MQ=255)
                                   gAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCa                                                                                                           >  2:656677/1‑143 (MQ=255)
                                   gAGCCATTTGTTTCCTAGTTTATCGTTTAGGCGTTCAATTTAAGTAGCTCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGcc                                                                                                            >  2:643265/1‑143 (MQ=255)
                                   gAGCCATTTGTTTCATCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCGCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTGGGTTGCCGTCGCCTTGTTTCCTTAGGcc                                                                                                            >  1:841789/1‑143 (MQ=255)
                                          ttGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGACCGCCTTGTAGGCTTCGGCCATACGGt                                                                                                     >  2:501217/1‑143 (MQ=255)
                                          ttGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGt                                                                                                     <  1:85085/143‑1 (MQ=255)
                                           tGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTg                                                                                                    <  1:568151/143‑1 (MQ=255)
                                              ttCGTCGGTTATCTTTTAGGCGTGCAGTTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAAc                                                                                                 <  2:857986/143‑1 (MQ=255)
                                              ttCCTCGTTTATCTTGTAGGCGTTCAAGTTAAGTAGCTCAAAGCGGGCTGCTCACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAAc                                                                                                 <  2:316104/143‑1 (MQ=255)
                                               tCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAAGTCCGAAGGGATAACCAACGGTAGTGTGCGAGCGCCTTGTTGGCTTCGGCCATACGGTGAACg                                                                                                >  1:532935/1‑143 (MQ=255)
                                                  tCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTTAAACCCCAAATGGTTTGTTAACACAACAATGAGAACGCCTTGTTGGCTTCGGCCAAACGATTAACGTCt                                                                                             >  1:197439/1‑143 (MQ=255)
                                                  tCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCt                                                                                             >  1:515596/1‑143 (MQ=255)
                                                     tttATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTGCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGACCAGACGGTGGACGTCGTCa                                                                                          >  2:147628/1‑143 (MQ=255)
                                                      ttATGTTTTATGCGTTCAATTCCAGCAGCCCAAAGCGGGCTGCTGACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCAc                                                                                         <  1:592585/143‑1 (MQ=255)
                                                        aTCTTTGAGGCGTTCAATGTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGt                                                                                       <  1:457240/143‑1 (MQ=255)
                                                         tCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGtt                                                                                      >  1:912308/1‑143 (MQ=255)
                                                           ttttAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTc                                                                                    <  2:246607/143‑1 (MQ=255)
                                                           ttttAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTc                                                                                    <  2:615123/143‑1 (MQ=255)
                                                           ttttAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGGCTTcgggttgc                                                                                    >  1:78494/1‑136 (MQ=255)
                                                               aGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTaa                                                                                <  1:160282/143‑1 (MQ=255)
                                                                 ccGTTCAATTTCACTTGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGACGGCATCCTCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACt                                                                              <  2:794245/142‑1 (MQ=255)
                                                                    ttCAATGTAACCAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCa                                                                           <  1:1037054/143‑1 (MQ=255)
                                                                     tCAATTTAGGTAGCCCACAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAGCTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAg                                                                          <  2:54863/143‑1 (MQ=255)
                                                                       ggtgtaggTAGCCGAAAGCGGGGTGCTCACTGGCAGCGCGCGCCGGGTGGCTAAGACGCCGAAGGGATACCCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTa                                                                        <  2:691652/136‑1 (MQ=255)
                                                                              gTAGCCCAAAGCGGGCCGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGt                                                                 <  1:715810/143‑1 (MQ=255)
                                                                                  cccAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTc                                                             >  2:1052598/1‑143 (MQ=255)
                                                                                  cacAAACCGCGCTGCGTACTGGAAGCGCCCGGCTGGTGACTAAATCTCCGAAGGGATAACCAACGGTAGTGTGCGATCGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTCCACGTTTCTTAACTGCAGTACCTTTGTTTTc                                                             <  1:140074/141‑1 (MQ=255)
                                                                                    cAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTg                                                           <  1:778982/143‑1 (MQ=255)
                                                                                       cgccggcggcTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCGACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTCTCTTAACTGCAGTACCTTTGTTTTCTgcag                                                        <  2:524251/139‑1 (MQ=255)
                                                                                             ctgctTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTGTCTGCAGCAACAg                                                  >  1:805070/1‑143 (MQ=255)
                                                                                                ctTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTCTGTGTGCTGCAGCATCAGata                                               >  2:1012573/1‑141 (MQ=255)
                                                                                                 ttACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCCCGTTTCTTACCTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAg                                              <  1:827885/143‑1 (MQ=255)
                                                                                                 ttACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAg                                              >  1:602753/1‑143 (MQ=255)
                                                                                                  tACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgt                                             <  1:538404/143‑1 (MQ=255)
                                                                                                    cTGGAAGCGCCCGCCTGGCGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGaaaaata                                           >  1:707488/1‑139 (MQ=255)
                                                                                                     tGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTTTCTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttcg                                          >  2:1018302/1‑143 (MQ=255)
                                                                                                     tGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttcg                                          >  1:567393/1‑143 (MQ=255)
                                                                                                        aaGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGGATCAGAAGgttcgttc                                       >  1:7801/1‑143 (MQ=255)
                                                                                                        aaGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttcgttc                                       >  1:25703/1‑143 (MQ=255)
                                                                                                          gcgcCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGc                                     >  1:574314/1‑143 (MQ=255)
                                                                                                           cgcccgccTAGCGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCCCGGTGTGTACGGTGGTCGTCTTCACGCACCCTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGcc                                    <  1:490059/143‑1 (MQ=255)
                                                                                                            tcccgccTGGTGACTAAAACGCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCa                                   <  1:448597/142‑1 (MQ=255)
                                                                                                               cgccTGGTGACTAAATCTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGgc                                <  1:643265/143‑1 (MQ=255)
                                                                                                                  cTGGTGACTAAAGCTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCa                             <  2:446136/143‑1 (MQ=255)
                                                                                                                  cTGGTGACTAAAACTCCCAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCa                             <  1:83891/143‑1 (MQ=255)
                                                                                                                          tAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAGCGTCTTCACGTTTCGTAACTGCAGTCCCTTTGTTTTCTGCAGCATCAGAAAGTCCGTTCGACAGGCGCAGAGcccag                     >  2:618549/1‑140 (MQ=255)
                                                                                                                          gaagccgccGAAGGGATAACCGGCGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATg                     <  2:59148/136‑1 (MQ=255)
                                                                                                                               cTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGAttt                <  1:572314/143‑1 (MQ=255)
                                                                                                                                  cGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATc             >  2:944051/1‑143 (MQ=255)
                                                                                                                                            aCCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACCGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGTTTa   <  1:557617/143‑1 (MQ=255)
                                                                                                                                             ccAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCCTGGATTTATCACCGCGtttcc  >  1:56026/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                              
ACCGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGTTTAC  >  NC_000913/3473473‑3473756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: