Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,214,964 A→G I229V (ATT→GTT)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,9640AG99.2% 425.1 / NA 121I229V (ATT→GTT) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (61/59);  total (62/59)
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.00e+00

Only 100 of 121 total aligned reads displayed.

TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCG  >  NC_000913/4214868‑4215059
                                                                                                |                                                                                               
tttCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTg                                                                                              <  2:943707/100‑1 (MQ=255)
   ccGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCt                                                                                           >  2:33478/1‑100 (MQ=255)
    cGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCtt                                                                                          >  2:2826505/1‑100 (MQ=255)
     ggCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCttt                                                                                         >  2:446736/1‑100 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                        <  1:141009/100‑1 (MQ=255)
        aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtgtg                                                                                      >  2:2343564/1‑100 (MQ=255)
        aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                        <  1:2504187/98‑1 (MQ=255)
          cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGAc                                                                                    >  2:344820/1‑100 (MQ=255)
          cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGAc                                                                                    >  1:2072124/1‑100 (MQ=255)
           gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACg                                                                                   >  2:2174248/1‑100 (MQ=255)
             tGATTCGTGATTACACCGATCTGGAAATTATGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                 <  1:1210892/100‑1 (MQ=255)
              gATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                <  1:27751/100‑1 (MQ=255)
               aTTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGcc                                                                               >  2:120909/1‑100 (MQ=255)
               aTTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGcc                                                                               >  1:94665/1‑100 (MQ=255)
                    tGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATccc                                                                          <  1:1058285/100‑1 (MQ=255)
                    tGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATccc                                                                          >  1:595947/1‑100 (MQ=255)
                    tGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAAATAAGCGCGTTGCCTTTGTGACGGGCCATccc                                                                          >  2:503071/1‑100 (MQ=255)
                     gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCg                                                                         >  2:1458968/1‑100 (MQ=255)
                      aTTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGa                                                                        >  1:901244/1‑100 (MQ=255)
                        tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAat                                                                      <  1:1903052/100‑1 (MQ=255)
                        tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAat                                                                      <  2:2660007/100‑1 (MQ=255)
                         acacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                     >  2:2645077/1‑100 (MQ=255)
                         acacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                     >  1:1399418/1‑100 (MQ=255)
                         acacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                     <  1:1609927/100‑1 (MQ=255)
                         acacCGATCTGGAAATTCTGGCAGAGACAGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGagat                                                                     >  1:211633/1‑97 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                    >  1:1413930/1‑100 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                    >  2:32477/1‑100 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                    >  1:236960/1‑100 (MQ=255)
                           acCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatg                                                                   <  2:697728/100‑1 (MQ=255)
                           acCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatg                                                                   >  2:190229/1‑100 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                  <  2:2187067/100‑1 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                  <  2:466565/100‑1 (MQ=255)
                                  tGGAAATTCCGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCa                                                            >  2:1976294/1‑100 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCCGTAGAGATAAGCGCGTTGCCTTTGTGACTGGCCATCCCGAATATGATGCGCa                                                            >  2:1501129/1‑99 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                           >  1:2032161/1‑100 (MQ=255)
                                    gAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaaa                                                          >  2:1988813/1‑100 (MQ=255)
                                     aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                         <  1:2110010/100‑1 (MQ=255)
                                     aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                         <  1:644543/100‑1 (MQ=255)
                                     aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                         >  2:2782904/1‑100 (MQ=255)
                                       aTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGc                                                       >  1:2941100/1‑100 (MQ=255)
                                        ttCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCt                                                      >  2:1931781/1‑100 (MQ=255)
                                          cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                    <  1:1805933/100‑1 (MQ=255)
                                            ggCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGgcg                                                  <  1:33479/100‑1 (MQ=255)
                                             gCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGgcgc                                                 >  2:1429912/1‑100 (MQ=255)
                                              cAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCa                                                >  1:1733311/1‑100 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                               <  2:17908/100‑1 (MQ=255)
                                                gagaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCa                                                <  2:2749978/98‑1 (MQ=255)
                                                gagaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCACGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAgg                                              <  1:697880/100‑1 (MQ=255)
                                                  gaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGaa                                            >  1:2004728/1‑100 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                          >  1:2354161/1‑100 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                          <  2:798284/100‑1 (MQ=255)
                                                      gaagaaGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtttt                                        >  1:2548525/1‑100 (MQ=255)
                                                        agaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                      >  2:1902502/1‑100 (MQ=255)
                                                         gaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTcc                                     <  1:2907074/100‑1 (MQ=255)
                                                         gaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGAAGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                      >  1:1459257/1‑99 (MQ=255)
                                                         gaagGGGATGCATAACTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTcc                                     <  1:540406/100‑1 (MQ=255)
                                                           agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgc                                   >  2:2068174/1‑100 (MQ=255)
                                                           agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgc                                   >  2:1139990/1‑100 (MQ=255)
                                                            ggggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgcg                                  <  2:2032143/100‑1 (MQ=255)
                                                             gggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGAGCCATCCCGCATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGa                                 >  2:2326726/1‑100 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                <  2:385086/100‑1 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                >  1:21817/1‑100 (MQ=255)
                                                               gATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGTGCAAACGCTGGCGCAGGAATTTTTCCGCGAtg                               <  1:2874845/100‑1 (MQ=255)
                                                                 tGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGAAAACGCTGGCGCAGGAATTTTTCCGAGAtgtg                             >  1:248532/1‑100 (MQ=255)
                                                                  tcttatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTgg                            <  2:2905764/97‑1 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTTCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                           >  1:2875383/1‑100 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGaga                         >  1:2905790/1‑98 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         >  1:587959/1‑100 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         <  1:1577532/100‑1 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         <  1:120910/100‑1 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGGGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         <  2:1093240/100‑1 (MQ=255)
                                                                       tCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                       <  2:94664/100‑1 (MQ=255)
                                                                         tGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                     <  1:2645100/100‑1 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGGCGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                     <  2:1399406/99‑1 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                    >  2:1065452/1‑100 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGan                   <  1:758006/100‑2 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                   <  1:2242564/100‑1 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCACGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                   <  1:445165/100‑1 (MQ=255)
                                                                            ttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACt                  <  1:1526835/100‑1 (MQ=255)
                                                                            ttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACt                  <  1:1429924/100‑1 (MQ=255)
                                                                            ttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACt                  <  2:2371101/100‑1 (MQ=255)
                                                                             tGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                 >  1:768868/1‑100 (MQ=255)
                                                                             tGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                 <  2:236957/100‑1 (MQ=255)
                                                                                 aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGGCCGACTAGAcc             >  1:843616/1‑100 (MQ=255)
                                                                                 aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc             >  2:139214/1‑100 (MQ=255)
                                                                                 aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc             >  1:2257830/1‑100 (MQ=255)
                                                                                 aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc             >  1:1369561/1‑100 (MQ=255)
                                                                                    aaaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCgg          >  2:215810/1‑100 (MQ=255)
                                                                                    aaaGATAAGCGCATTGGCTTTGTGACGGGCTATACCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCTGGACTAGAc              >  2:2187655/1‑96 (MQ=255)
                                                                                     gatctAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa         <  1:619385/96‑1 (MQ=255)
                                                                                     aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa         <  1:344825/100‑1 (MQ=255)
                                                                                     aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa         >  1:307802/1‑100 (MQ=255)
                                                                                     aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTGTTCCGCGATGTGGAAGCCGGACTAGACCCGGa         <  2:1883104/100‑1 (MQ=255)
                                                                                      aGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt        <  1:1922201/100‑1 (MQ=255)
                                                                                       gATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg       <  2:2178650/100‑1 (MQ=255)
                                                                                         tAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTa     >  2:820875/1‑100 (MQ=255)
                                                                                         tAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCAGGATGTa     >  1:1007564/1‑100 (MQ=255)
                                                                                          aaGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGaagatc    >  2:619378/1‑95 (MQ=255)
                                                                                            gcgcGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCg  <  1:2253633/100‑1 (MQ=255)
                                                                                            gcgcGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCg  >  1:117001/1‑100 (MQ=255)
                                                                                                |                                                                                               
TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCG  >  NC_000913/4214868‑4215059

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