Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,388 G→T 12.1% P186Q (CCG→CAG)  pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,797,3880GT12.1% 197.4 / 14.5 80P186Q (CCG→CAG) pheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base T (4/6);  total (36/44)
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.81e-01

TTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGC  >  NC_000913/1797240‑1797534
                                                                                                                                                    |                                                                                                                                                  
tttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATTGTCGCACCAACCGGAACGATATCCGGTTGAACAAGCgg                                                                                                                                                   >  2:508660/1‑150 (MQ=255)
  tcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCCTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACAAGCggaa                                                                                                                                                 >  2:130581/1‑148 (MQ=255)
   cATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAg                                                                                                                                                <  1:841961/150‑1 (MQ=255)
     tCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTACTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                                                                              <  1:846402/150‑1 (MQ=255)
      ccACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                                                                             >  2:1163329/1‑150 (MQ=255)
      ccACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                                                                             >  1:584095/1‑150 (MQ=255)
         caGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGtt                                                                                                                                          <  1:311695/150‑1 (MQ=255)
         caGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt                                                                                                                                          >  1:96828/1‑150 (MQ=255)
         caGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt                                                                                                                                          <  2:766281/150‑1 (MQ=255)
          tgCGGAGTTGGCGCTTTAAGGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTTAACCAGCGGCAGCTGGTTc                                                                                                                                         <  1:1641237/149‑1 (MQ=255)
          cgCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTc                                                                                                                                         <  2:810494/149‑1 (MQ=255)
          aGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGGCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                                                                         <  1:1578506/150‑1 (MQ=255)
           gCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                                                                        >  1:533552/1‑150 (MQ=255)
              gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                                                                     <  1:1434800/150‑4 (MQ=255)
              gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                                                                     <  1:268098/150‑4 (MQ=255)
              gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGCTGAACCAGCTGCCGCTGGTTCaacc                                                                                                                                     <  1:640378/150‑4 (MQ=255)
              gAGGTGGCGCTTTAACGTTAAGGCCTTTTACCACACGGCCAAGATACCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCa                                                                                                                                     <  1:1376140/150‑1 (MQ=255)
                gTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACg                                                                                                                                   >  2:1257778/1‑150 (MQ=255)
                  tGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGc                                                                                                                                 <  1:383503/150‑1 (MQ=255)
                        tttAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTc                                                                                                                           >  1:1688990/1‑150 (MQ=255)
                           aaCGTTAATGCCTTTTACCCCACGGCCAATATAAGGCGGGCTGGCGTATGGCGCTTCGACTGTTATCGGCAGCGTGTCGTCGATGGTGGCTCCAACCGGAACGTTTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcg                                                                                                                        <  2:1533691/150‑1 (MQ=255)
                             cGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgcg                                                                                                                      >  2:385994/1‑150 (MQ=255)
                                 aaTGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAac                                                                                                                  >  2:1081575/1‑150 (MQ=255)
                                     ccTTTTACCCCACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGGCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCa                                                                                                              <  2:140716/150‑1 (MQ=255)
                                     ccTTTTACCACCCGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCTGCAGCGTGTCGTCGATGGTCTCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCa                                                                                                              <  2:135485/150‑1 (MQ=255)
                                        tttACCACACGGCCAAGATATCGCGGGCAGGCTTCCTGCGCGTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACTCGCAACACCAatg                                                                                                           <  2:832927/150‑1 (MQ=255)
                                           aCCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                                                                        <  2:445995/150‑1 (MQ=255)
                                           aCCACACGGCCAAGATAACGCGGGCAGGCGTCCGGCGCTTCGACGGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                                                                        <  2:1234003/150‑1 (MQ=255)
                                             cacaCGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                      <  2:195420/150‑1 (MQ=255)
                                              acacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGTCAGCGTGTCGTCGATGGTCGCACCAACCTGAACGATTTCCGCTTGAACCAGAGGCAGCTGGTTAAGCACGGCAACGTCACGCGCAACACCAATGATGACt                                                                                                     >  1:181364/1‑150 (MQ=255)
                                              acacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCt                                                                                                     >  2:1205561/1‑150 (MQ=255)
                                                  gtccAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGGCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                 <  1:433975/148‑1 (MQ=255)
                                                   gCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCTCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGAGCAACACGAATGATGCCTAAGCa                                                                                                >  1:844569/1‑150 (MQ=255)
                                                   gCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACCCCAATGATGCCTAAGCa                                                                                                >  1:1476944/1‑150 (MQ=255)
                                                       aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCg                                                                                            >  2:190154/1‑150 (MQ=255)
                                                       aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCATGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCg                                                                                            >  2:946289/1‑150 (MQ=255)
                                                        gATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCgg                                                                                           <  1:1469086/150‑1 (MQ=255)
                                                           aaCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCCGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                        >  2:614705/1‑150 (MQ=255)
                                                               cgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCCCGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGtt                                                                                    >  1:438954/1‑150 (MQ=255)
                                                                gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                                                                   >  1:562496/1‑150 (MQ=255)
                                                                   cAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGc                                                                                <  2:1192827/150‑1 (MQ=255)
                                                                    aGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                                               >  1:775041/1‑150 (MQ=255)
                                                                      gCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTc                                                                             <  2:1407551/150‑1 (MQ=255)
                                                                           cGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGATCCAGCGGCAGCTGGTTCAGCACGGCCACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCCCGGTTTGGCGTCACGCt                                                                        >  1:968946/1‑150 (MQ=255)
                                                                            ggCGCTGCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTg                                                                       <  1:1299140/150‑1 (MQ=255)
                                                                               gcTTCGACTGTAATCGGCAGAGTGTCGTCGATGGTCGCACCAACGGGATCGATTTCCGGTTGAACCAGCTGCAGCTGTTTCTGCACGGATACGTCACGCGCAACACTAATGATGCCTAAGCAGTCGGCTCGGTTTGGGGTGACGCTGAtt                                                                    <  2:967913/150‑1 (MQ=255)
                                                                                 ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                                                                  <  1:1163329/150‑1 (MQ=255)
                                                                                  tCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCg                                                                 <  2:63630/150‑1 (MQ=255)
                                                                                    gACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                                                               <  2:1060171/150‑1 (MQ=255)
                                                                                     aCTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg                                                              <  2:1344789/150‑1 (MQ=255)
                                                                                        ttAATCGGCAGCGTGTCGTCGAGGGTCACACCATCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCATCACGGCAACGTCACGCGCAACACCAATTATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg                                                           <  2:1390654/149‑1 (MQ=255)
                                                                                        gTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg                                                           >  1:383749/1‑150 (MQ=255)
                                                                                          aaTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgtt                                                         >  2:414578/1‑150 (MQ=255)
                                                                                               gCAGCGTGTCGTCGATCGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGGCACGCTGATTTCGATGGTGTTGtcat                                                    >  1:1632306/1‑150 (MQ=255)
                                                                                                cAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtcatc                                                   <  1:385994/150‑1 (MQ=255)
                                                                                                  gCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCaa                                                 <  1:190154/150‑1 (MQ=255)
                                                                                                   cGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAg                                                <  1:1448003/150‑1 (MQ=255)
                                                                                                        cgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGGAGCTGGTTCAGCACGGCAACGTCACGCGCAACAACAATGATGCCTAAGCAGTCGGCACGGTATGGACTCACGCTGATATCGATGGTGTTGTCATCAAGTTTCa                                           >  1:239914/1‑150 (MQ=255)
                                                                                                        cgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCa                                           <  1:680934/150‑1 (MQ=255)
                                                                                                         gtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAg                                          >  2:1267851/1‑150 (MQ=255)
                                                                                                          tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAgg                                         <  2:790360/150‑1 (MQ=255)
                                                                                                          tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAgg                                         <  1:1081575/150‑1 (MQ=255)
                                                                                                           cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGt                                        >  2:993616/1‑150 (MQ=255)
                                                                                                             aTGGTCGCCCCACCCGCAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATTATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTAt                                      <  1:956317/150‑1 (MQ=255)
                                                                                                               ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTc                                    >  2:1711445/1‑150 (MQ=255)
                                                                                                                  cGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACg                                 <  2:228515/150‑1 (MQ=255)
                                                                                                                     ttcAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGat                              >  2:1358523/3‑150 (MQ=255)
                                                                                                                        aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCTACTTCACGCGCAACACCACTGATGCCTAAGCAGTCGGCACGGTTTGGCGGCACGCTGATTTCGTTTGTGTTGTCATCTAGTTTCAGGTATTCACGGCTATc                           <  2:1636101/150‑1 (MQ=255)
                                                                                                                           cGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCggc                        <  1:945234/150‑2 (MQ=255)
                                                                                                                             gAACGATTTACGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGc                      >  1:1719996/1‑150 (MQ=255)
                                                                                                                              aaCGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGcc                     >  2:1393000/1‑150 (MQ=255)
                                                                                                                                  aTTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATc                 <  1:614705/150‑1 (MQ=255)
                                                                                                                                   tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCg                >  1:1620340/1‑150 (MQ=255)
                                                                                                                                   tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCg                <  1:1205561/150‑1 (MQ=255)
                                                                                                                                     ccgtgctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACCCCAATGATGCCTAAGCAGTCGGCACGGGTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgc              >  2:231344/7‑150 (MQ=255)
                                                                                                                                      cgtgctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCCCGCTGATTTCGATGGTGTTGTCATCTAGTTTCAGTTATTCACGGATATCGGTGCCAATCGgcg             >  2:1050232/6‑150 (MQ=255)
                                                                                                                                      ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgcg             <  1:159842/150‑1 (MQ=255)
                                                                                                                                        tgctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCa           >  1:410048/4‑150 (MQ=255)
                                                                                                                                        tgctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGGTTGGCGTCACGCTGATGTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCa           >  1:259451/4‑150 (MQ=255)
                                                                                                                                            gAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCg       >  1:1058428/1‑150 (MQ=255)
                                                                                                                                              aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgca     >  2:1622173/1‑150 (MQ=255)
                                                                                                                                                 aGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGAGGCCTACGCAGTCGGCACGGGTTGGCGTCACGCTGATTTCGATGGTGCTGTCATCAAGTTTCAGGGATTCACGGATATCGGTGCCAATCGGCGCATCCgcaggc  <  1:858616/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGC  >  NC_000913/1797240‑1797534

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