Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,103 C→T 8.4% H112Y (CAT→TAT)  hyaD → putative hydrogenase 1 maturation protease HyaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,1030CT 8.4% 213.2 / 7.3 71H112Y (CAT→TAT) hyaDputative hydrogenase 1 maturation protease HyaD
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGC  >  NC_000913/1035955‑1036250
                                                                                                                                                    |                                                                                                                                                   
gATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATTAGCCTGCATCAGAACAGTTTCTCCGATGTGTTGGCGCTGGCGGATATCCGCGGAct                                                                                                                                                    <  2:1600907/150‑2 (MQ=255)
   tctcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCt                                                                                                                                                 >  2:491863/1‑150 (MQ=255)
   tctcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATcc                                                                                                                                                 <  1:1704130/150‑2 (MQ=255)
    ctcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                                                                >  2:902562/1‑150 (MQ=255)
       gATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCa                                                                                                                                             <  1:534184/150‑1 (MQ=255)
        aTGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTTCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAg                                                                                                                                            <  1:220890/150‑1 (MQ=255)
        aTGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAg                                                                                                                                            <  1:1186469/150‑1 (MQ=255)
        aTGCCATTGACTACGGGCTGGAACCTGGAACCCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACATTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAg                                                                                                                                            <  2:1981304/150‑1 (MQ=255)
         tGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGc                                                                                                                                           <  1:997465/150‑1 (MQ=255)
          gCCATTGACTACGGGCTGGAACCTGGAACGCTGCGTACCTATGCCGGAGAACGCTTTCCGGCTTATCTCAGCGCGATGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                                                                          <  2:832404/150‑2 (MQ=255)
           ccATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcac                                                                                                                                         >  1:1642299/1‑150 (MQ=255)
                gACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGATGTGTTTGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                                                                                    <  2:1196537/150‑1 (MQ=255)
                gACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                                                                                    >  2:1226242/1‑150 (MQ=255)
                       ggCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcg                                                                                                                             <  2:1817106/150‑1 (MQ=255)
                       ggCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcg                                                                                                                             >  1:7587/1‑150 (MQ=255)
                           ggaacctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCgg                                                                                                                         >  2:1944543/1‑150 (MQ=255)
                           ggaacctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCgg                                                                                                                         >  1:20484/1‑150 (MQ=255)
                           ggaacctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCgg                                                                                                                         >  2:1042150/1‑150 (MQ=255)
                                ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                                                                    >  2:595360/1‑150 (MQ=255)
                                ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                                                                    >  2:655838/1‑150 (MQ=255)
                                   gaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAc                                                                                                                 <  2:1560465/150‑1 (MQ=255)
                                       gCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGc                                                                                                             <  1:112245/150‑1 (MQ=255)
                                        cTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCa                                                                                                            >  2:1673406/1‑150 (MQ=255)
                                           cGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATg                                                                                                         <  2:499686/150‑1 (MQ=255)
                                               ccTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                                                                                     <  1:1042150/150‑1 (MQ=255)
                                                 tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                                                                   <  1:491863/150‑1 (MQ=255)
                                                 tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                                                                   <  1:469579/150‑1 (MQ=255)
                                                  aTGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACa                                                                                                                                                    <  1:1160107/100‑1 (MQ=255)
                                                  aTGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACa                                                                                                                                                    >  2:1160107/1‑100 (MQ=255)
                                                  aTGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgacg                                                                                                  >  2:195091/1‑150 (MQ=255)
                                                     ccGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACt                                                                                               <  1:1374908/150‑1 (MQ=255)
                                                          gaACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggc                                                                                          >  2:602236/1‑150 (MQ=255)
                                                            aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                                                                        <  1:263631/150‑1 (MQ=255)
                                                             cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGt                                                                                       >  2:1071588/1‑150 (MQ=255)
                                                                   ccGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                                                                 >  2:417654/1‑150 (MQ=255)
                                                                      gCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                                                                              <  1:902562/150‑1 (MQ=255)
                                                                          aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                                                                          >  2:1012119/1‑150 (MQ=255)
                                                                          aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                                                                          <  1:1877474/150‑1 (MQ=255)
                                                                              cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGa                                                                                <  1:1212434/140‑1 (MQ=255)
                                                                              cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGa                                                                                >  2:1212434/1‑140 (MQ=255)
                                                                              cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc                                                                      <  2:1942584/150‑1 (MQ=255)
                                                                                 cgcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACg                                                                   >  2:242138/1‑150 (MQ=255)
                                                                                  gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACgg                                                                  >  1:1908599/1‑150 (MQ=255)
                                                                                    gaagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGa                                                                >  1:1591408/1‑150 (MQ=255)
                                                                                     aagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg                                                               <  1:595360/150‑1 (MQ=255)
                                                                                       gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTGGACGACTACGGCGGTAGCCTGAGCGAACTGGCACTGGAGCa                                                             <  2:7587/150‑1 (MQ=255)
                                                                                        aaaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCaa                                                            >  2:1552488/1‑150 (MQ=255)
                                                                                                    cATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGTCCGCTGCg                                                >  2:399376/1‑150 (MQ=255)
                                                                                                       cAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGaa                                             <  2:188822/150‑1 (MQ=255)
                                                                                                         gAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCCACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACa                                           >  2:1705988/1‑150 (MQ=255)
                                                                                                                  ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCt                                  >  2:1913283/1‑150 (MQ=255)
                                                                                                                      gAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGgcg                              <  2:1642299/150‑1 (MQ=255)
                                                                                                                      gAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGgcg                              <  2:127473/150‑1 (MQ=255)
                                                                                                                         gtgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAg                           <  1:245966/150‑1 (MQ=255)
                                                                                                                          tgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGc                          >  2:1031468/1‑150 (MQ=255)
                                                                                                                          tggtGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGc                          <  1:863512/147‑1 (MQ=255)
                                                                                                                           gtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCATTGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCt                         <  2:1148716/150‑1 (MQ=255)
                                                                                                                              ggCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgc                      <  1:1334042/150‑1 (MQ=255)
                                                                                                                                tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctg                    >  1:1208016/2‑150 (MQ=255)
                                                                                                                                tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctg                    >  2:1771303/2‑150 (MQ=255)
                                                                                                                                 gcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGGCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctgc                   <  2:1340668/150‑1 (MQ=255)
                                                                                                                                  cTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCg                  >  1:89042/1‑150 (MQ=255)
                                                                                                                                     gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTgg               <  1:172471/150‑1 (MQ=255)
                                                                                                                                     gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgcca               <  2:89042/150‑4 (MQ=255)
                                                                                                                                      cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGcaag                       <  1:702018/141‑4 (MQ=255)
                                                                                                                                      cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGcaag                       >  2:702018/1‑138 (MQ=255)
                                                                                                                                      cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag              <  2:20484/150‑5 (MQ=255)
                                                                                                                                      cGGATAGCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag              <  1:911810/150‑5 (MQ=255)
                                                                                                                                           atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAtt         >  2:402357/1‑150 (MQ=255)
                                                                                                                                               gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGc     >  2:145806/1‑150 (MQ=255)
                                                                                                                                               gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGc     <  1:1226242/150‑1 (MQ=255)
                                                                                                                                                  gACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc  >  2:1088697/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGC  >  NC_000913/1035955‑1036250

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: