Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,618,274 C→T G212D (GGC→GAC)  marC ← UPF0056 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,618,2740CT98.6% 261.3 / ‑4.6 70G212D (GGC→GAC) marCUPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (35/34);  total (35/35)
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.98e-01

ATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAA  >  NC_000913/1618127‑1618417
                                                                                                                                                   |                                                                                                                                               
aTGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGt                                                                                                                                              <  2:72333/151‑1 (MQ=255)
  gtgGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGtt                                                                                                                                             <  2:409208/150‑1 (MQ=255)
  gtgGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGtt                                                                                                                                             <  2:121560/150‑1 (MQ=255)
    gggCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaat                                                                                                                                          <  1:558422/151‑1 (MQ=255)
     ggCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaat                                                                                                                                          >  1:151387/1‑150 (MQ=255)
        gcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaataata                                                                                                                                      <  1:745344/151‑1 (MQ=255)
           gTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAAc                                                                                                                                   >  1:755264/1‑151 (MQ=255)
               cTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGc                                                                                                                                >  1:15434/1‑150 (MQ=255)
                tGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGcac                                                                                                                              >  1:141541/1‑151 (MQ=255)
                 gCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGcac                                                                                                                              >  2:769698/1‑150 (MQ=255)
                  cTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCacac                                                                                                                            >  1:501809/1‑151 (MQ=255)
                  attttGCCGCCTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGTTGTCGTTAATAATAAACTGCacac                                                                                                                            >  1:581203/2‑151 (MQ=255)
                       gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCcat                                                                                                                        <  1:799798/150‑1 (MQ=255)
                         cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCcataca                                                                                                                     <  1:434409/151‑1 (MQ=255)
                            gTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCatacatac                                                                                                                  <  1:163692/151‑1 (MQ=255)
                             ttAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATAcc                                                                                                                 >  2:279765/1‑151 (MQ=255)
                              tAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATAcc                                                                                                                 <  1:165527/150‑1 (MQ=255)
                                aTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATATTAAACTGCACACCCATACATACcagc                                                                                                              >  1:492582/1‑151 (MQ=255)
                                 tttGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACcagc                                                                                                              >  2:190774/1‑150 (MQ=255)
                                  ttGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACcagc                                                                                                              >  2:90008/1‑149 (MQ=255)
                                       cAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATc                                                                                                       >  2:655504/1‑151 (MQ=255)
                                           cATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                                                                                                   <  1:148170/151‑1 (MQ=255)
                                           cATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                                                                                                   <  1:240032/151‑1 (MQ=255)
                                            aTTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                                  >  2:695014/1‑151 (MQ=255)
                                               atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                                  >  1:697961/1‑148 (MQ=255)
                                               atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAg                                                                                                <  1:58947/150‑1 (MQ=255)
                                               atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAg                                                                                                <  2:512359/150‑1 (MQ=255)
                                                 atTCCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGAc                                                                                              >  1:637330/1‑150 (MQ=255)
                                                       ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATTCCAGCAGGAATCCCATCAGACGGGAGAt                                                                                       <  1:384800/151‑1 (MQ=255)
                                                                 gTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATg                                                                             >  2:607332/1‑151 (MQ=255)
                                                                      acTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAcc                                                                        >  2:576886/1‑151 (MQ=255)
                                                                      acTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAcc                                                                        >  1:798395/1‑151 (MQ=255)
                                                                         cGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAccc                                                                       <  2:832556/149‑1 (MQ=255)
                                                                         cGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAccc                                                                       <  1:89985/149‑1 (MQ=255)
                                                                         cGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCtt                                                                     <  2:121243/151‑1 (MQ=255)
                                                                          gaagaaCAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTg                                                                    >  2:291018/1‑151 (MQ=255)
                                                                                cAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacc                                                              <  2:521565/151‑1 (MQ=255)
                                                                                 aGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacca                                                             <  2:249888/151‑1 (MQ=255)
                                                                                    cGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGcc                                                          >  2:129297/1‑151 (MQ=255)
                                                                                      cAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGc                                                        <  1:300141/151‑1 (MQ=255)
                                                                                      cAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGc                                                        >  2:344949/1‑151 (MQ=255)
                                                                                         tAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                      >  1:566870/1‑150 (MQ=255)
                                                                                            ttGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGCAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATa                                                     <  1:503622/148‑1 (MQ=255)
                                                                                             tGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCtt                                                                     >  2:431904/1‑131 (MQ=255)
                                                                                             tGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCtt                                                                     <  1:431871/131‑1 (MQ=255)
                                                                                             tGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAAtt                                                  <  1:531013/150‑1 (MQ=255)
                                                                                              gAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTgc                                                <  2:19203/151‑1 (MQ=255)
                                                                                               aaaGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTgcg                                               <  2:420918/151‑1 (MQ=255)
                                                                                                  ggCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCg                                            >  2:444322/1‑151 (MQ=255)
                                                                                                    cccATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                           >  2:89693/1‑150 (MQ=255)
                                                                                                     ccATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                           >  1:438678/1‑149 (MQ=255)
                                                                                                         tCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTg                                      <  1:693300/150‑1 (MQ=255)
                                                                                                          cGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCg                                    <  1:612660/151‑1 (MQ=255)
                                                                                                            ggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGt                                   >  2:24384/1‑150 (MQ=255)
                                                                                                            ggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTa                                  <  1:293339/151‑1 (MQ=255)
                                                                                                            ggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTa                                  <  2:20416/151‑1 (MQ=255)
                                                                                                              ccTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAgg                                >  2:715939/1‑151 (MQ=255)
                                                                                                               cTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGc                               >  2:533982/1‑151 (MQ=255)
                                                                                                                       aTGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                       <  1:731722/151‑1 (MQ=255)
                                                                                                                         ggTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaa                     <  2:30153/151‑1 (MQ=255)
                                                                                                                            aCGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAAt                  >  2:44534/1‑151 (MQ=255)
                                                                                                                              gTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATaa                <  1:87494/151‑1 (MQ=255)
                                                                                                                              gTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATa                 >  1:688890/1‑150 (MQ=255)
                                                                                                                                   aaTGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcc           >  1:175272/1‑151 (MQ=255)
                                                                                                                                    aTGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcca          <  1:590640/151‑1 (MQ=255)
                                                                                                                                          tCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAgg     <  2:759701/150‑1 (MQ=255)
                                                                                                                                           ccAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaa   <  1:484614/151‑1 (MQ=255)
                                                                                                                                            cAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaa   >  1:795127/1‑150 (MQ=255)
                                                                                                                                             aGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaaa  >  2:518293/1‑150 (MQ=255)
                                                                                                                                             aGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaaa  >  2:356048/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAA  >  NC_000913/1618127‑1618417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑