Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,610,720 T→A 8.5% L436Q (CTG→CAG)  prfC → peptide chain release factor RF3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,610,7200TA 8.5% 171.5 / 6.2 59L436Q (CTG→CAG) prfCpeptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GTTCACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAG  >  NC_000913/4610571‑4610865
                                                                                                                                                     |                                                                                                                                                 
gTTCACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCt                                                                                                                                                   <  2:419254/150‑1 (MQ=255)
gTTCACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCt                                                                                                                                                   >  1:1215244/1‑150 (MQ=255)
  tCACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGa                                                                                                                                                 >  2:291527/1‑150 (MQ=255)
   cACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGAt                                                                                                                                                >  1:857515/1‑150 (MQ=255)
       ggTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                            <  2:851671/150‑1 (MQ=255)
         tATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTgg                                                                                                                                          <  1:291527/150‑1 (MQ=255)
            tCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGc                                                                                                                                       <  1:851881/150‑1 (MQ=255)
            tCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGc                                                                                                                                       >  2:851881/1‑150 (MQ=255)
            tCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGc                                                                                                                                       >  2:898475/1‑150 (MQ=255)
             ccGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCACCAACGATCAGATCGttgttgga                                                                                                                                      <  2:956306/150‑3 (MQ=255)
             ccGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCa                                                                                                                                      >  1:403529/1‑150 (MQ=255)
              cGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAg                                                                                                                                     >  1:1190759/1‑150 (MQ=255)
              cGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGTTGTTGGAg                                                                                                                                     <  2:1287906/150‑1 (MQ=255)
              cGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGGGTTCCGTCCAATCTCCACCAACGATCAGATCGTTGGTGGAg                                                                                                                                     <  2:1044302/150‑1 (MQ=255)
                  cTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTgg                                                                                                                                 <  1:1353238/150‑1 (MQ=255)
                     cGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGCCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGgtgt                                                                                                                              >  2:394650/1‑150 (MQ=255)
                        aCCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGtgct                                                                                                                           >  1:1335303/1‑150 (MQ=255)
                         ccAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGtgctg                                                                                                                          >  2:68148/1‑150 (MQ=255)
                               cTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                    <  1:979466/150‑1 (MQ=255)
                                 gTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGa                                                                                                                  >  1:1283127/1‑150 (MQ=255)
                                   tCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGAtg                                                                                                                <  2:856646/150‑1 (MQ=255)
                                      gtcgtATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATgtgg                                                                                                             <  2:645237/150‑1 (MQ=255)
                                          tATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTggtggt                                                                                                         <  2:1000445/150‑1 (MQ=255)
                                                  tGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGcc                                                                                                 <  2:959780/150‑1 (MQ=255)
                                                   gAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCt                                                                                                >  1:898086/1‑150 (MQ=255)
                                                   gAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCt                                                                                                >  1:902919/1‑150 (MQ=255)
                                                        aTCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAaga                                                                                           >  1:1093107/1‑150 (MQ=255)
                                                          ccGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAg                                                                                                       <  1:990677/136‑1 (MQ=255)
                                                          ccGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAg                                                                                                       >  2:990677/1‑136 (MQ=255)
                                                           cGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCg                                                                                        >  1:1132657/1‑150 (MQ=255)
                                                           cGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCg                                                                                        <  2:403529/150‑1 (MQ=255)
                                                              tGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAAt                                                                                     <  1:141553/150‑1 (MQ=255)
                                                                              gctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgt                                                                     >  1:1080110/1‑150 (MQ=255)
                                                                                 cAAAGGGCTGGTACAGCTTTCCCAAGAGGGCGCGGTGCAGGTGGTCAGTCCAAACTCCAACAACGATCTGATCCATGGTGCAGTTGGTGTGCTTCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTa                                                                  >  2:607279/1‑150 (MQ=255)
                                                                                   aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATg                                                                <  2:1132657/150‑1 (MQ=255)
                                                                                     gggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGag                                                              >  2:397732/1‑150 (MQ=255)
                                                                                     gggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAACAGCGAATACCACGTTGAAGCAGTGTATGag                                                              >  2:634309/1‑150 (MQ=255)
                                                                                             aCAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCaa                                                      <  2:368075/150‑1 (MQ=255)
                                                                                              cAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAAc                                                     >  2:345433/1‑150 (MQ=255)
                                                                                               aGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACg                                                    >  2:942254/1‑150 (MQ=255)
                                                                                                  tttCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTg                                                 >  1:1077204/1‑150 (MQ=255)
                                                                                                     ccGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCa                                              >  2:1188979/1‑150 (MQ=255)
                                                                                                        aaGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTg                                           <  2:1111266/150‑1 (MQ=255)
                                                                                                                 cgGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTggg                                  <  1:397732/150‑1 (MQ=255)
                                                                                                                   gTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTa                                <  2:1215244/150‑1 (MQ=255)
                                                                                                                     gCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                              >  1:534701/1‑150 (MQ=255)
                                                                                                                       aGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAAt                            <  2:1093107/150‑1 (MQ=255)
                                                                                                                            ttCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCTAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                       <  2:902919/150‑1 (MQ=255)
                                                                                                                            ttCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                       <  2:898086/150‑1 (MQ=255)
                                                                                                                            gcgactCCAATGTCGATCAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                       <  1:1208664/145‑1 (MQ=255)
                                                                                                                              ccGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGa                     >  2:367912/1‑150 (MQ=255)
                                                                                                                               cGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCACTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGAc                    <  2:1190759/150‑1 (MQ=255)
                                                                                                                                    aaTCTCCAACAACGATCTGAGCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCgaa               <  2:705792/150‑1 (MQ=255)
                                                                                                                                      tctcCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCgaaga             <  2:335463/150‑1 (MQ=255)
                                                                                                                                         gcaccaacGATCAGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAAt          >  2:161310/5‑150 (MQ=255)
                                                                                                                                            acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCg       <  2:1283127/150‑1 (MQ=255)
                                                                                                                                              aacGATCAGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGaa     >  1:1417556/1‑150 (MQ=255)
                                                                                                                                               acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAag    <  1:640275/150‑1 (MQ=255)
                                                                                                                                                cGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAaga   >  1:99822/1‑150 (MQ=255)
                                                                                                                                                 gATCTGATCGTTGGTGCAGTTGGTGTGTTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAagag  >  2:1003002/1‑150 (MQ=255)
                                                                                                                                                 gATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAagag  <  2:1335303/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTCACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAG  >  NC_000913/4610571‑4610865

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