Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,974,266 A→G 5.2% N583N (AAT→AAC)  aas ← acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,974,2660AG 5.2% 297.6 / 4.9 95N583N (AAT→AAC) aasacyltransferase
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base G (3/2);  total (48/47)
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.44e-01

GCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAAT  >  NC_000913/2974117‑2974409
                                                                                                                                                     |                                                                                                                                               
gCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACa                                                                                                                                                 >  1:1355414/1‑150 (MQ=255)
   aGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATtc                                                                                                                                              >  1:985069/1‑150 (MQ=255)
     ttgttCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATtctc                                                                                                                                            <  1:69930/150‑1 (MQ=255)
        ttCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTcggc                                                                                                                                         >  2:636352/1‑150 (MQ=255)
          ccaccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTcggcgg                                                                                                                                       >  2:995683/1‑150 (MQ=255)
          ccaccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTcggcgg                                                                                                                                       >  1:1469565/1‑150 (MQ=255)
          ccaccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCACCGCGAACATTCTcggcgg                                                                                                                                       >  2:1336508/1‑150 (MQ=255)
           caccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGt                                                                                                                                      >  1:1916810/1‑150 (MQ=255)
             ccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTgg                                                                                                                                    <  1:636867/150‑1 (MQ=255)
             ccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTgg                                                                                                                                    <  2:1947023/150‑1 (MQ=255)
             ccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGAGAATGTTCgcggcgaaat                                                                                                                                    <  1:39773/150‑5 (MQ=255)
             ccaTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGAGAATGTTCgcggcgaaat                                                                                                                                    <  1:23311/150‑5 (MQ=255)
               aTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGc                                                                                                                                  <  2:1103064/150‑1 (MQ=255)
               aTTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGc                                                                                                                                  <  1:1638319/150‑1 (MQ=255)
                  tCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACt                                                                                                                               <  1:1988790/150‑1 (MQ=255)
                   ccAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACtt                                                                                                                              >  2:1334099/1‑150 (MQ=255)
                    cAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTc                                                                                                                             >  2:922687/1‑150 (MQ=255)
                       cGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAg                                                                                                                          >  1:1837083/1‑150 (MQ=255)
                        gACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGt                                                                                                                         >  1:1719788/1‑150 (MQ=255)
                        gACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGt                                                                                                                         <  2:441615/150‑1 (MQ=255)
                            ccATTTCGCCTGCAATTTTGGCAAAGCGATTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCAAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTacac                                                                                                                     >  2:137430/1‑150 (MQ=255)
                                 tccccTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGc                                                                                                                <  2:841357/147‑1 (MQ=255)
                                      tGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTtctc                                                                                                           >  2:1623166/1‑150 (MQ=255)
                                        cAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCa                                                                                                         <  2:373733/150‑1 (MQ=255)
                                          aTTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCAcc                                                                                                       <  2:1631655/150‑1 (MQ=255)
                                            tttGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCg                                                                                                     <  2:799346/150‑1 (MQ=255)
                                             ttGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGc                                                                                                    <  1:1725765/150‑1 (MQ=255)
                                               ggCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAg                                                                                                  >  1:304354/1‑150 (MQ=255)
                                                gCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGa                                                                                                 >  2:115690/1‑150 (MQ=255)
                                                gCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGa                                                                                                 >  1:1471971/1‑150 (MQ=255)
                                                 cAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGAt                                                                                                >  2:1156765/1‑150 (MQ=255)
                                                         tttGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTc                                                                                        <  1:1064344/150‑1 (MQ=255)
                                                           tGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCAt                                                                                      >  1:560957/1‑150 (MQ=255)
                                                           tGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCAt                                                                                      <  2:1431533/150‑1 (MQ=255)
                                                            gcgcGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCAtt                                                                                     <  2:1719788/150‑1 (MQ=255)
                                                            gcgcGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCAtt                                                                                     >  2:1866925/1‑150 (MQ=255)
                                                             cgaatcccTGAATCTGCAAAAAGCCATGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCTCAGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTa                                                                                    >  2:1659471/6‑150 (MQ=255)
                                                                ggCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGt                                                                                 <  2:1940166/150‑1 (MQ=255)
                                                                   ccTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCg                                                                              >  2:604978/1‑150 (MQ=255)
                                                                     tGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGc                                                                            <  2:1585491/150‑1 (MQ=255)
                                                                     tGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGc                                                                            <  1:636352/150‑1 (MQ=255)
                                                                       aaTCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGccc                                                                          <  2:1112283/150‑1 (MQ=255)
                                                                          cTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCttt                                                                       >  1:1421840/1‑150 (MQ=255)
                                                                          cTGCACAAAGCCCTGCTCGTCAAAACCCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCttt                                                                       >  2:1798650/1‑150 (MQ=255)
                                                                           tGCTCAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTc                                                                      >  1:1765996/1‑150 (MQ=255)
                                                                           tGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCGAGTGTCATACCAGCCTCGCTCCATTTCACCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTc                                                                      <  1:1336508/150‑1 (MQ=255)
                                                                           tGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTc                                                                      >  2:952673/1‑150 (MQ=255)
                                                                            gCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCa                                                                     >  1:1962507/1‑150 (MQ=255)
                                                                             cacaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAg                                                                    >  1:990692/1‑150 (MQ=255)
                                                                             cacaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAg                                                                    >  2:643157/1‑150 (MQ=255)
                                                                              acaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGt                                                                   <  2:1469565/150‑1 (MQ=255)
                                                                                aaaGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTg                                                                 <  1:1866925/150‑1 (MQ=255)
                                                                                 aaGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGc                                                                >  2:1517585/1‑150 (MQ=255)
                                                                                    cccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAgg                                                             >  2:1032905/1‑150 (MQ=255)
                                                                                    cccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAgg                                                             <  1:995683/150‑1 (MQ=255)
                                                                                    cccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAgg                                                             >  1:824262/1‑150 (MQ=255)
                                                                                     ccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGc                                                            <  1:824084/150‑1 (MQ=255)
                                                                                     ccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCCAACATTCTCGGCGGTGGGCACTTCCACTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGc                                                            >  2:384617/1‑150 (MQ=255)
                                                                                           cGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGc                                                      <  2:1734809/150‑1 (MQ=255)
                                                                                            gTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGcc                                                     >  1:928877/1‑150 (MQ=255)
                                                                                                 aaCGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCtt                                                <  1:1032905/150‑1 (MQ=255)
                                                                                                   cGCACAATATCGCCAGTGTCATACCAGCCGCGGTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCg                                              <  1:1329835/150‑1 (MQ=255)
                                                                                                   cGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCg                                              >  2:483622/1‑150 (MQ=255)
                                                                                                   cGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCg                                              >  2:1950706/1‑150 (MQ=255)
                                                                                                   cGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCg                                              <  1:1398415/150‑1 (MQ=255)
                                                                                                    gCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGa                                             <  2:17939/150‑1 (MQ=255)
                                                                                                     cacaATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGAt                                            <  2:1962507/150‑1 (MQ=255)
                                                                                                     cacaATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGAt                                            <  2:1765996/150‑1 (MQ=255)
                                                                                                       caATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATAc                                          <  1:606248/150‑1 (MQ=255)
                                                                                                             cGCCAGTGTCATACCAGCCGCGCTCCATATCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGAAGGCGTCCGCCCTCTTCGATACCAGGGa                                    >  2:1519298/1‑150 (MQ=255)
                                                                                                                cAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCg                                 <  1:1818260/150‑1 (MQ=255)
                                                                                                                  gtgtCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGac                               <  2:1837083/150‑1 (MQ=255)
                                                                                                                    gtCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCCCCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGacca                             >  2:1060331/1‑148 (MQ=255)
                                                                                                                    gtCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGacaa                             <  2:928877/150‑1 (MQ=255)
                                                                                                                     gcATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGacaac                            <  1:1334099/149‑1 (MQ=255)
                                                                                                                      cATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGacaaca                           >  1:89968/1‑150 (MQ=255)
                                                                                                                        tACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAgg                         >  2:1927926/1‑150 (MQ=255)
                                                                                                                        tACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAgg                         >  1:1078897/1‑150 (MQ=255)
                                                                                                                          ccAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcg                       >  1:576127/1‑150 (MQ=255)
                                                                                                                          ccAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcg                       >  2:1536854/1‑150 (MQ=255)
                                                                                                                          ccAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcg                       >  1:159988/1‑150 (MQ=255)
                                                                                                                            aGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcgcg                     <  1:1517585/150‑1 (MQ=255)
                                                                                                                            aGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcgcg                     <  2:824262/150‑1 (MQ=255)
                                                                                                                              ccGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCa                   <  2:1421840/150‑1 (MQ=255)
                                                                                                                              ccGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCa                   <  1:1156765/150‑1 (MQ=255)
                                                                                                                                 cgcTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATcc                <  1:1623166/150‑1 (MQ=255)
                                                                                                                                   cTCCATTTCGCCGCGCACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCCCCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAt              >  2:1653973/1‑150 (MQ=255)
                                                                                                                                   cTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAt              >  1:37164/1‑150 (MQ=255)
                                                                                                                                    tgcccaccgccgaGAATGTTCGCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAtt             >  2:648178/8‑150 (MQ=255)
                                                                                                                                       atatCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCt          <  1:1519298/147‑1 (MQ=255)
                                                                                                                                        tttCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTg         <  2:1916810/150‑1 (MQ=255)
                                                                                                                                          ccgccgaGAATGTTCGCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGt       >  2:198455/2‑150 (MQ=255)
                                                                                                                                          ccgccgaGAATGTTCGCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGt       >  2:10596/2‑150 (MQ=255)
                                                                                                                                            gccgcGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAg     >  2:1032520/1‑150 (MQ=255)
                                                                                                                                            gccgcGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAg     <  1:952673/150‑1 (MQ=255)
                                                                                                                                               gcgAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCGCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAAt  <  1:137430/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAAT  >  NC_000913/2974117‑2974409

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: ‑