Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 442,389 G→A 5.2% V281I (GTA→ATA)  thiI → tRNA uridine 4‑sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913442,3890GA 5.2% 245.4 / 2.1 76V281I (GTA→ATA) thiItRNA uridine 4‑sulfurtransferase
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.31e-01

TTGGCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGCACATCATCAACC  >  NC_000913/442240‑442527
                                                                                                                                                     |                                                                                                                                          
ttGGCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCg                                                                                                                                            <  2:1155675/150‑1 (MQ=255)
ttGGCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACAGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTAAGATGGGCGTTATCCTCAAACGAATGATGGTGCGTGCAGCATATAAAGTGGCTGAACGTTACTGCg                                                                                                                                            >  2:953497/1‑150 (MQ=255)
  ggCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTa                                                                                                                                          >  2:796269/1‑150 (MQ=255)
  ggCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTa                                                                                                                                          >  2:791258/1‑150 (MQ=255)
   gCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTAc                                                                                                                                         >  1:572429/1‑150 (MQ=255)
       cTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGgc                                                                                                                                     >  1:1118122/1‑150 (MQ=255)
        tCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCGGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGgcg                                                                                                                                    >  1:790457/1‑150 (MQ=255)
          ccACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCATACAGGCGCt                                                                                                                                  >  1:469316/1‑150 (MQ=255)
          ccACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCATACAGGCGCt                                                                                                                                  >  1:1286367/1‑150 (MQ=255)
            aCCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTgg                                                                                                                                <  2:218059/150‑1 (MQ=255)
               gcgtgcgtTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCa                                                                                                                             >  1:883110/1‑150 (MQ=255)
                        ttGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAg                                                                                                                    <  1:1308260/150‑1 (MQ=255)
                         tGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAgc                                                                                                                   >  1:889842/1‑150 (MQ=255)
                          gTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAgcg                                                                                                                  <  1:287491/150‑1 (MQ=255)
                            cGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCt                                                                                                                >  1:989576/1‑150 (MQ=255)
                               tATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCgg                                                                                                             <  1:1081079/150‑1 (MQ=255)
                                 ttAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGcc                                                                                                           >  1:97359/1‑150 (MQ=255)
                                 ttAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCATACAGGCGCTGGTCACCGGCGAAGCGCTCGGcc                                                                                                           <  2:469316/150‑1 (MQ=255)
                                 ttAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCATACAGGCGCTGGTCACCGGCGAAGCGCTCGGcc                                                                                                           <  2:1286367/150‑1 (MQ=255)
                                    aTTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAgg                                                                                                        >  1:406432/1‑150 (MQ=255)
                                     tttCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGgt                                                                                                       <  1:802286/150‑1 (MQ=255)
                                      ttCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGgtg                                                                                                      >  1:390058/1‑150 (MQ=255)
                                        cGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTc                                                                                                    <  1:179209/150‑1 (MQ=255)
                                         gAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTcc                                                                                                   <  2:1118122/150‑1 (MQ=255)
                                          aaCCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTcca                                                                                                  <  1:65065/150‑1 (MQ=255)
                                              ggTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTccagcca                                                                                              >  1:194217/1‑150 (MQ=255)
                                                 cgtcgGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGAc                                                                                           >  2:1150234/1‑150 (MQ=255)
                                                  gtcgGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAgg                                                                                                        >  1:1175604/1‑136 (MQ=255)
                                                  gtcgGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAgg                                                                                                        <  2:1175604/136‑1 (MQ=255)
                                                   tcgGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGc                                                                                         <  2:63995/150‑1 (MQ=255)
                                                     ggggAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTgg                                                                                                                                <  1:186356/109‑1 (MQ=255)
                                                     ggggAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTgg                                                                                                                                >  2:186356/1‑109 (MQ=255)
                                                      gggAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTgc                                                                                      >  1:1377418/1‑149 (MQ=255)
                                                      gggAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCCCTGGTCACCCGCGAAACGCTCGGCCACGTCTCCACCAAGACCCTGa                                                                                      >  1:609251/1‑150 (MQ=255)
                                                          aaTTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCaa                                                                                  <  2:1125300/150‑1 (MQ=255)
                                                            ttCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAAcc                                                                                <  1:644661/150‑1 (MQ=255)
                                                                 gagaAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTgcgc                                                                           <  1:325826/150‑1 (MQ=255)
                                                                  agaAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGcc                                                                          <  2:390058/150‑1 (MQ=255)
                                                                     aaaTCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtga                                                                       >  2:665532/1‑150 (MQ=255)
                                                                        tCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtgattg                                                                    >  2:53938/1‑150 (MQ=255)
                                                                        tCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtgattg                                                                    <  1:1176670/150‑1 (MQ=255)
                                                                          gacgacGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtgattgat                                                                  >  2:205043/1‑150 (MQ=255)
                                                                             gacgGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAAc                                                               <  2:883110/150‑1 (MQ=255)
                                                                              acgGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACg                                                              <  1:194806/150‑1 (MQ=255)
                                                                                   cAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTcc                                                         <  1:1409326/150‑1 (MQ=255)
                                                                                    aGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCg                                                        >  2:945858/1‑150 (MQ=255)
                                                                                       tgcgCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGaca                                                     <  1:1271766/147‑1 (MQ=255)
                                                                                       tGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGaca                                                     >  2:21215/1‑150 (MQ=255)
                                                                                        gggCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGacac                                                    >  1:1324004/1‑150 (MQ=255)
                                                                                         ggCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACg                                                   >  1:504757/1‑150 (MQ=255)
                                                                                         ggCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACg                                                   >  2:86895/1‑150 (MQ=255)
                                                                                                 ccTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGAc                                                                                     >  2:142397/1‑108 (MQ=255)
                                                                                                 ccTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGAc                                                                                     <  1:142397/108‑1 (MQ=255)
                                                                                                 ccTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtgattgat                                                                  >  1:626808/1‑127 (MQ=255)
                                                                                                 ccTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCGCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCtgattgat                                                                  <  2:626808/127‑1 (MQ=255)
                                                                                                  cTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTg                                          >  2:295665/1‑150 (MQ=255)
                                                                                                    cAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGCCCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCg                                        >  2:625147/1‑150 (MQ=255)
                                                                                                    cAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCg                                        >  2:1251688/1‑150 (MQ=255)
                                                                                                      aaCGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTc                                      >  1:472799/1‑150 (MQ=255)
                                                                                                        cGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCg                                    >  2:453640/1‑150 (MQ=255)
                                                                                                             gatgGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGAt                               >  1:614836/1‑150 (MQ=255)
                                                                                                                   gcgtgcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTa                         >  1:1211305/1‑150 (MQ=255)
                                                                                                                    cgtgcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTAc                        <  1:205043/150‑1 (MQ=255)
                                                                                                                     gtgcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACg                       >  2:399799/1‑150 (MQ=255)
                                                                                                                     gtgcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACg                       <  2:572429/150‑1 (MQ=255)
                                                                                                                      tgcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGa                      >  2:597877/1‑150 (MQ=255)
                                                                                                                       gcCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGAc                     >  2:1074273/1‑150 (MQ=255)
                                                                                                                        ccGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACa                    <  2:406432/150‑1 (MQ=255)
                                                                                                                        ccGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACa                    >  2:432581/1‑150 (MQ=255)
                                                                                                                         cGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACaa                   >  2:994680/1‑150 (MQ=255)
                                                                                                                         cGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACaa                   <  2:1211305/150‑1 (MQ=255)
                                                                                                                             tCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAagag               <  1:53938/150‑1 (MQ=255)
                                                                                                                              cTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGc              >  1:1053676/1‑150 (MQ=255)
                                                                                                                                   gTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGCAcatc         <  1:86895/150‑1 (MQ=255)
                                                                                                                                      gCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGCACatcatc      >  2:764355/1‑150 (MQ=255)
                                                                                                                                          aaCGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGCACATCATCAAcc  >  1:363298/1‑150 (MQ=255)
                                                                                                                                                aCGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACg                                                              <  2:515111/84‑1 (MQ=255)
                                                                                                                                                aCGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACg                                                              >  1:515111/1‑84 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TTGGCAGCTCCCACCGCGTGCGTTTTGTCGCTATTAATTTCGAACCGGTCGTCGGGGAAATTCTCGAGAAAATCGACGACGGTCAGATGGGCGTTATCCTCAAACGTATGATGGTGCGTGCCGCATCTAAAGTGGCTGAACGTTACGGCGTACAGGCGCTGGTCACCGGCGAAGCGCTCGGCCAGGTGTCCAGCCAGACGCTGACCAACCTGCGCCTGATTGATAACGTCTCCGACACGCTGATCCTGCGTCCGCTGATCTCTTACGACAAAGAGCACATCATCAACC  >  NC_000913/442240‑442527

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: