Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,871,503 A→C D76E (GAT→GAG)  dgoT ← D‑galactonate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,871,5030AC92.9% 518.6 / 14.7 168D76E (GAT→GAG) dgoTD‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (76/80);  total (83/85)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGATCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCGAACTC  >  NC_000913/3871405‑3871598
                                                                                                  |                                                                                               
ccAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGAt                                                                                                >  2:750543/1‑100 (MQ=255)
 cAATTCATGACATTAAACCGGTGGCCAAGCCCTCTATTAAAGTCTCCACTGCCGAGCCCAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTc                                                                                               >  2:750087/1‑100 (MQ=255)
   aTTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTa                                                                                             >  1:249614/1‑100 (MQ=255)
     tAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTaaa                                                                                           >  2:1223684/1‑100 (MQ=255)
      aaTGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTaaaa                                                                                          <  1:674285/100‑1 (MQ=255)
      aaTGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTaaaa                                                                                          <  2:3881092/100‑1 (MQ=255)
        tGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaac                                                                                        >  1:2479035/1‑100 (MQ=255)
          aCATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaacca                                                                                      >  1:2332909/1‑100 (MQ=255)
          aCATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaacca                                                                                      <  2:674621/100‑1 (MQ=255)
           cATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaaccaa                                                                                     >  1:1652026/1‑100 (MQ=255)
            aTTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaaccaac                                                                                    <  1:1376199/100‑1 (MQ=255)
            aTTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAaaccaac                                                                                    <  1:1836895/100‑1 (MQ=255)
                 gCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAAccgccg                                                                               >  2:2778208/1‑100 (MQ=255)
                  cccGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCgg                                                                              >  2:172062/1‑100 (MQ=255)
                   ccGTGGCAAAGCCCTGGAATAAAGTCGCCACTATCCAGCCAAATACCGCAATAAAATAAGTCGCGCGAGAACCTACGCGCTCTAAAAACCCACCGCCggg                                                                             <  2:1758880/100‑1 (MQ=255)
                    cGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCgggg                                                                            <  1:153997/100‑1 (MQ=255)
                     gTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGGACCTACGCGATCTAAAAACCAACCGCCGGGGa                                                                           >  1:1215715/1‑100 (MQ=255)
                     gTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGa                                                                           >  1:1900589/1‑100 (MQ=255)
                        gCAAAGCCCTGGAATAAAGTCGCCACTGACCATCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAATCCAACCGCCGGGGATCt                                                                        <  1:3734961/100‑1 (MQ=255)
                        gCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCt                                                                        <  1:945080/100‑1 (MQ=255)
                         cAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTg                                                                       >  1:736839/1‑100 (MQ=255)
                          aaaGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGa                                                                      <  2:1707071/100‑1 (MQ=255)
                           aaGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGAc                                                                     <  2:3489064/100‑1 (MQ=255)
                           aaGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGAc                                                                     >  1:635802/1‑100 (MQ=255)
                           aaGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGAc                                                                     <  2:3276758/100‑1 (MQ=255)
                             gCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACAt                                                                   <  1:2820809/100‑1 (MQ=255)
                                cTGGAATAAAGTCGCCACTGACCAGCCAAATATAGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGc                                                                >  1:3851133/1‑100 (MQ=255)
                                 tGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCg                                                               >  1:577081/1‑100 (MQ=255)
                                 tGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCg                                                               <  1:601324/100‑1 (MQ=255)
                                  ggAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGt                                                              <  1:3667164/100‑1 (MQ=255)
                                  ggAATAAAGACGCCAATGCCCAGCTAAATAAAGCAATAAAATAAATCACGCGAGAACATACGAGACCTAAAAATAAACCGCCGGGGATCTGACAtcgaga                                                              >  2:3882570/1‑95 (MQ=255)
                                    aaTAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGtat                                                            <  1:3499985/100‑1 (MQ=255)
                                     aTAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGtata                                                           <  1:1262102/100‑1 (MQ=255)
                                     aTAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGtata                                                           <  2:3841158/100‑1 (MQ=255)
                                         aGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGc                                                       <  1:2989886/100‑1 (MQ=255)
                                          gTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGcc                                                      >  2:1099958/1‑100 (MQ=255)
                                          gTCGCCACTGACCAGCCAAATATCGCAACAAAATAAGTCGCGCGAGAACCTACGCGATCTAAAAACCAACCGCCGGGGATCTGAAATAGCGTATAAAGcc                                                      <  2:3045416/100‑1 (MQ=255)
                                           tCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCa                                                     <  1:2528417/100‑1 (MQ=255)
                                             gCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAACCCAACCGCAGTGAATCTGACATAGGGGATACAGCCaaa                                                   >  2:3280622/1‑98 (MQ=255)
                                             gCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGCCATAGCGTATAAAGCCAgg                                                   >  2:3579398/1‑100 (MQ=255)
                                             gCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAgg                                                   <  1:3504939/100‑1 (MQ=255)
                                             gCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAgg                                                   >  2:3550397/1‑100 (MQ=255)
                                              ccACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGc                                                  >  2:2301884/1‑100 (MQ=255)
                                               cACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCg                                                 >  2:3822639/1‑100 (MQ=255)
                                                aCTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGa                                                <  1:2084980/100‑1 (MQ=255)
                                                 cTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGaa                                               >  2:1436889/1‑100 (MQ=255)
                                                  tGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAg                                              <  1:4033937/100‑1 (MQ=255)
                                                  tGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAg                                              <  2:3244925/100‑1 (MQ=255)
                                                   gACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAgg                                             >  1:2848538/1‑100 (MQ=255)
                                                    aCCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCATCCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGc                                            <  1:1819050/100‑1 (MQ=255)
                                                     ccAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAATCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGcc                                           >  1:2148103/1‑100 (MQ=255)
                                                     ccAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGcc                                           >  1:3824644/1‑100 (MQ=255)
                                                     ccAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGcc                                           >  2:1212261/1‑100 (MQ=255)
                                                        gCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGaa                                        >  2:1415155/1‑100 (MQ=255)
                                                          cAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGaaaa                                      >  1:3804325/1‑100 (MQ=255)
                                                          cAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGaaaa                                      >  1:3587296/1‑100 (MQ=255)
                                                           aaaTATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAAt                                     <  1:1523259/100‑1 (MQ=255)
                                                            aaTATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATa                                    <  1:43427/100‑1 (MQ=255)
                                                            aaTATCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATa                                    <  1:3512289/100‑1 (MQ=255)
                                                             atatCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATAc                                   <  2:2440928/100‑1 (MQ=255)
                                                              tatCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACa                                  <  1:2218059/100‑1 (MQ=255)
                                                              tatCGCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACa                                  <  1:1868200/100‑1 (MQ=255)
                                                                  gcgttaacataaGGCACGCGAGAACCTACGCGCGTTAAAAGCCAGCCGCCGGGCATCTGACAGAGCGTTTAAAGCCATGCGAAGGCCGAAAATACATAGc                              <  2:3946253/96‑1 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAACATACATCGc                              >  1:3986654/1‑100 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              <  1:1030252/100‑1 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              <  1:229732/100‑1 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              >  2:832514/1‑100 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              >  1:2988043/1‑100 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              >  1:3124887/1‑100 (MQ=255)
                                                                  gCAATAAAATAAGTCACGCGAGAACCTACGCGATCTAAAAACAAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGc                              >  1:522820/1‑100 (MQ=255)
                                                                   cAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGcc                             >  2:2042207/1‑100 (MQ=255)
                                                                   cAATAAAATAAGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGcc                             <  1:3613518/100‑1 (MQ=255)
                                                                        aaataaGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCAttt                        >  2:1635215/1‑100 (MQ=255)
                                                                        aaataaGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCAttt                        <  2:3724424/100‑1 (MQ=255)
                                                                         aataaGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTc                       >  2:55418/1‑100 (MQ=255)
                                                                         aagaaGTCACGCGAGAACCTACGCGCGTTAAAAACCAACCGCCGGGGATCGGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTc                       <  2:3503398/97‑1 (MQ=255)
                                                                          ataaGTCACGCAACAACCTAAGCGCTCTAAATACCATCCTCCGGGGATCTGACATAGCGTATAAAGCCAGACGAAGGCCGAAAATACATAGCCCATTTcc                      <  2:733665/100‑1 (MQ=255)
                                                                             aGTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCt                   >  2:2536847/1‑100 (MQ=255)
                                                                              gTCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCtt                  >  2:1563742/1‑100 (MQ=255)
                                                                               tCACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCttt                 <  2:655171/100‑1 (MQ=255)
                                                                                cACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTg                >  1:3376049/1‑100 (MQ=255)
                                                                                cACGCGAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTg                >  2:2903321/1‑100 (MQ=255)
                                                                                 aCGCGAGAACCTACGCGATCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTgg               <  2:2153740/100‑1 (MQ=255)
                                                                                  cgcgAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGt              <  2:2128000/100‑1 (MQ=255)
                                                                                  cgcgAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGt              <  2:1030152/100‑1 (MQ=255)
                                                                                  cgcgAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGt              >  1:4006067/1‑100 (MQ=255)
                                                                                  cgcgAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGt              >  2:3107138/1‑100 (MQ=255)
                                                                                  cgcgAGAACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATACCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGt              <  2:1133773/100‑1 (MQ=255)
                                                                                   gcgAGAACCTACGCGCTCTAAAAACCACCCGCCGGTGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTa             <  1:767793/100‑1 (MQ=255)
                                                                                     gagaACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAAt           <  1:2184272/100‑1 (MQ=255)
                                                                                     gagaACCTACGCGATCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAAt           <  1:2649628/100‑1 (MQ=255)
                                                                                       gaACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGc         >  2:377573/1‑100 (MQ=255)
                                                                                       gaACCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGc         <  1:3904214/100‑1 (MQ=255)
                                                                                        aaCCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGcc        <  1:2913236/100‑1 (MQ=255)
                                                                                        aaCCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGcc        <  1:1392373/100‑1 (MQ=255)
                                                                                         aCCTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCg       >  1:1046836/1‑100 (MQ=255)
                                                                                          ccTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCGa      >  1:4074505/1‑100 (MQ=255)
                                                                                          ccTACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCGa      >  2:878276/1‑100 (MQ=255)
                                                                                            tACGCGCTCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCGAAc    <  1:2362195/100‑1 (MQ=255)
                                                                                             atgcgcTCTAAAAACCAACCGCCGGGGATATGACATAGCGTATAAAGACAGGCGAATGCCGAAAATACATAGCCCATCTCCGCTTTGGTAATGCCGAAct   >  2:250697/3‑100 (MQ=255)
                                                                                                  |                                                                                               
CCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAATAAAATAAGTCACGCGAGAACCTACGCGATCTAAAAACCAACCGCCGGGGATCTGACATAGCGTATAAAGCCAGGCGAAGGCCGAAAATACATAGCCCATTTCCGCTTTGGTAATGCCGAACTC  >  NC_000913/3871405‑3871598

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