Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,807 C→G 8.8% G515R (GGT→CGT)  tynA ← copper‑containing amine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,8070CG 8.8% 1033.0 / 71.2 332G515R (GGT→CGT) tynAcopper‑containing amine oxidase
Reads supporting (aligned to +/- strand):  ref base C (152/149);  new base G (15/14);  total (167/165)
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.41e-03

Only 100 of 339 total aligned reads displayed.

CAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTCCGGC  >  NC_000913/1449658‑1449955
                                                                                                                                                     |                                                                                                                                                    
cAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCCTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAcc                                                                                                                                                      >  1:212/1‑150 (MQ=255)
 aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                                     <  2:2273372/150‑1 (MQ=255)
     cGAGGCGGAAATTATAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATc                                                                                                                                                >  2:23291/1‑150 (MQ=255)
       aGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCg                                                                                                                                               <  1:253513/150‑1 (MQ=255)
        ggCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGa                                                                                                                                              <  2:883766/150‑1 (MQ=255)
        ggCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGa                                                                                                                                              >  1:742024/1‑150 (MQ=255)
          cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGaag                                                                                                                                            <  1:478244/150‑3 (MQ=255)
           ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                           <  1:857472/150‑1 (MQ=255)
            gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGcc                                                                                                                                          <  2:1473745/150‑1 (MQ=255)
            gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCATCGATGCCGGTGCTACGGGCATCGAAGcg                                                                                                                                          <  1:2442821/150‑2 (MQ=255)
            gAAATTATAAATATGTTGGGGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTGGCCGTCTCATCGTGCATGGTTTTCGGTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGcc                                                                                                                                          <  2:1583056/150‑1 (MQ=255)
             aaaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcgg                                                                                                                                         <  1:34922/150‑3 (MQ=255)
             aaaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCa                                                                                                                                         <  2:718196/150‑1 (MQ=255)
             aaaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCATCGATGCCGGTGCTACGGGCATCGAAGcgg                                                                                                                                         <  1:1839565/150‑3 (MQ=255)
              aaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcggt                                                                                                                                        <  2:1163443/150‑4 (MQ=255)
              aaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcggt                                                                                                                                        <  2:230368/150‑4 (MQ=255)
               aTTATAAATATGTTGGGGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcggtg                                                                                                                                       <  1:1663536/150‑5 (MQ=255)
                    aaaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGc                                                                                                                                  <  1:917465/150‑1 (MQ=255)
                     aaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGcc                                                                                                                                 >  2:1942208/1‑150 (MQ=255)
                           ttGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTc                                                                                                                           >  1:2255108/1‑150 (MQ=255)
                            tGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCa                                                                                                                          <  2:1848934/150‑1 (MQ=255)
                            tGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCa                                                                                                                          >  2:949864/1‑150 (MQ=255)
                             ggTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCAt                                                                                                                         <  1:103217/150‑1 (MQ=255)
                                    agtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGAt                                                                                                                  >  1:1091562/1‑150 (MQ=255)
                                     gtaCCCACGATATTGTGATCGATAAGCGTGCAGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTACCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATc                                                                                                                 >  2:2174320/1‑150 (MQ=255)
                                      taCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATcc                                                                                                                >  2:753373/1‑150 (MQ=255)
                                      taCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATcc                                                                                                                <  2:67135/150‑1 (MQ=255)
                                       aCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                               >  2:2158848/1‑150 (MQ=255)
                                         ccACGATCTTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCTGCATCGATGCCAATAGTGCCGTTTTCATGGAATATCCAGt                                                                                                             <  2:377394/150‑1 (MQ=255)
                                         ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGt                                                                                                             >  1:629217/1‑150 (MQ=255)
                                          cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTc                                                                                                            <  1:1373429/150‑1 (MQ=255)
                                           aCGATATTGTGATCGATAAGCGGGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCa                                                                                                           <  2:1402262/150‑1 (MQ=255)
                                            cGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCaa                                                                                                          >  1:1134400/1‑150 (MQ=255)
                                                      aTCGATAAGCGTGCCGTAGCGCGTGTAATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                                <  2:333205/150‑1 (MQ=255)
                                                       tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCa                                                                                               >  1:1640061/1‑150 (MQ=255)
                                                         gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                             <  2:2356861/150‑1 (MQ=255)
                                                          aTAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAg                                                                                            <  2:1259034/150‑1 (MQ=255)
                                                            aaGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATTCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGtt                                                                                          >  2:1888386/1‑150 (MQ=255)
                                                             aGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTa                                                                                         >  1:2299459/1‑150 (MQ=255)
                                                               cGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAcc                                                                                       <  1:1605747/150‑1 (MQ=255)
                                                               cGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAcc                                                                                       >  2:509238/1‑150 (MQ=255)
                                                                gTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAccc                                                                                      <  1:2392439/150‑1 (MQ=255)
                                                                  gCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                                    <  1:2043286/150‑1 (MQ=255)
                                                                       aGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                               <  2:2476873/150‑1 (MQ=255)
                                                                        gcgcGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                              <  2:283847/150‑1 (MQ=255)
                                                                          gcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGa                                                                            <  2:1309595/150‑1 (MQ=255)
                                                                           cgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGAt                                                                           >  2:2370069/1‑150 (MQ=255)
                                                                            gtgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATc                                                                          <  2:1062934/150‑1 (MQ=255)
                                                                                cATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgc                                                                      <  1:355696/150‑1 (MQ=255)
                                                                                 aTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgcg                                                                     >  2:61430/1‑150 (MQ=255)
                                                                                    tttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcac                                                                  <  1:146646/150‑1 (MQ=255)
                                                                                     ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCATTCACCGCTTCGATGCCCGTAGCACCGGCATCGAAGCCAATAGTGCCGTTTTCATGGAAGATCCAGTAAAAAATGTAGTCATAGATACCCACTGTACTGATCCAGCGcacc                                                                 >  2:2449637/1‑150 (MQ=255)
                                                                                        gCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACt                                                              <  1:723925/150‑1 (MQ=255)
                                                                                         ccGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTa                                                             <  2:1638666/150‑1 (MQ=255)
                                                                                           gTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAAc                                                           >  1:528583/1‑150 (MQ=255)
                                                                                               cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCCCTCTACTGATCCAGCGCACCACTAACTccc                                                       >  2:2098916/1‑150 (MQ=255)
                                                                                               cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATACAGCGCACCACTAACTccc                                                       >  1:1697756/1‑150 (MQ=255)
                                                                                                aTCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCg                                                      >  2:451819/1‑150 (MQ=255)
                                                                                                 tCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCgg                                                     <  1:413152/150‑1 (MQ=255)
                                                                                                   ttGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCg                                                   <  1:1783123/149‑1 (MQ=255)
                                                                                                   gTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCg                                                   >  1:2406480/1‑150 (MQ=255)
                                                                                                      cATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgc                                                                      >  1:1820234/1‑128 (MQ=255)
                                                                                                      cATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTc                                                <  1:593863/150‑1 (MQ=255)
                                                                                                      cATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTc                                                >  1:2336877/1‑150 (MQ=255)
                                                                                                           ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTAc                                           >  2:2409599/1‑150 (MQ=255)
                                                                                                            tttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATcc                                                                                                                <  1:197190/80‑1 (MQ=255)
                                                                                                              tCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGa                                        >  2:54507/1‑150 (MQ=255)
                                                                                                                gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACg                                      >  1:544866/1‑150 (MQ=255)
                                                                                                                 cTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGt                                     >  2:170849/1‑150 (MQ=255)
                                                                                                                 cTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGt                                     >  1:698784/1‑150 (MQ=255)
                                                                                                                 cTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATACCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGt                                     >  1:1218805/1‑150 (MQ=255)
                                                                                                                   ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCAACACTAACTCCCGGCGTTCGGTACTGACGTTg                                   >  2:1987880/1‑150 (MQ=255)
                                                                                                                    tAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTgg                                  <  2:2257954/150‑1 (MQ=255)
                                                                                                                      acacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGc                                <  1:79558/150‑1 (MQ=255)
                                                                                                                      acacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGc                                >  1:1533434/1‑150 (MQ=255)
                                                                                                                       cacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCt                               >  2:31180/1‑150 (MQ=255)
                                                                                                                         acTTTCACGGCTTCGATGCCCGTAGAACCGGCATCGATGCGAATAGTGCAGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTgg                             <  2:1079206/149‑1 (MQ=255)
                                                                                                                          cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGc                            <  2:1335040/150‑1 (MQ=255)
                                                                                                                            ttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGccc                          <  2:1671023/150‑1 (MQ=255)
                                                                                                                              cACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAGAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt                        >  2:2024837/1‑150 (MQ=255)
                                                                                                                              cACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCCCTGTACTGATCCAGCGCACCCCTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt                        >  1:2089563/1‑150 (MQ=255)
                                                                                                                               aCCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAtt                       >  2:1530422/1‑150 (MQ=255)
                                                                                                                                   catCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt                   >  2:1833339/3‑150 (MQ=255)
                                                                                                                                   cTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt                   <  1:1287727/150‑1 (MQ=255)
                                                                                                                                    ttCGATGCCCGTAGCACCGGCATCGATGCCAATAGTCCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTg                  <  2:1764868/150‑1 (MQ=255)
                                                                                                                                    atCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTg                  >  1:1003195/2‑150 (MQ=255)
                                                                                                                                    atCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTg                  >  1:971826/2‑150 (MQ=255)
                                                                                                                                      cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt                >  1:1890390/1‑150 (MQ=255)
                                                                                                                                          gCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCtt            >  2:926537/1‑150 (MQ=255)
                                                                                                                                           cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta           >  1:1726645/1‑150 (MQ=255)
                                                                                                                                            ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGCAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtat          <  1:2356475/150‑1 (MQ=255)
                                                                                                                                             cGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATACAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtata         <  2:1697756/150‑1 (MQ=255)
                                                                                                                                              gTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATAc        <  1:991698/150‑1 (MQ=255)
                                                                                                                                                aGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTc      <  2:1818857/150‑1 (MQ=255)
                                                                                                                                                  ctccGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTccg    <  1:215428/148‑1 (MQ=255)
                                                                                                                                                  cACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTccg    >  1:757493/1‑150 (MQ=255)
                                                                                                                                                   aCCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTccgg   <  1:926537/150‑1 (MQ=255)
                                                                                                                                                   aCCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCAACACTAACTCCCGGCGTTCGGTACTGAACTTGTGCTGGCAAATTAACTAATGCTTAAACtgaga   >  1:1057775/1‑146 (MQ=255)
                                                                                                                                                   aCCGGCATCGATGCCAATAGTGCCGTTTTCATAGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTccgg   >  1:2227893/1‑150 (MQ=255)
                                                                                                                                                    ccGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTccggc  >  2:1748215/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTCCGGC  >  NC_000913/1449658‑1449955

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