Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,352,410 C→T I48I (ATC→ATT gltB → glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,352,4100CT99.0% 726.1 / ‑3.9 208I48I (ATC→ATTgltBglutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (118/88);  total (120/88)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 219 total aligned reads displayed.

GAAGTGGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATTCTCGCCG  >  NC_000913/3352321‑3352501
                                                                                         |                                                                                           
gAAGTGGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGcc                                                                                          <  2:4702336/93‑1 (MQ=255)
   gTGGGGTTCCCGCAGAGCCTGGGGGAGGCTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCcac                                                                                       >  1:3296045/1‑93 (MQ=255)
    tGGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCcaca                                                                                      >  2:782617/1‑93 (MQ=255)
     ggggTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACAt                                                                                     <  1:4245173/93‑1 (MQ=255)
     ggggTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACAt                                                                                     <  1:4263419/93‑1 (MQ=255)
      ggggTCCCGGCGAGCCTGGGGGAGGGTCACGATATGTTGGACGATAAATCCCTTGAGAGGGATAACTGGGGGTTTGGGCTTATTGGCCCCATa                                                                                    <  2:4609440/93‑1 (MQ=255)
       ggTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAg                                                                                   <  1:3640969/93‑1 (MQ=255)
          tCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAg                                                                                <  1:4463146/93‑1 (MQ=255)
           cccGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAgg                                                                               >  2:3816873/1‑93 (MQ=255)
             cGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGGTAACTGTGGTTTCGGCCGGATTGCCCACATAGAAGGCg                                                                             >  1:604734/1‑93 (MQ=255)
                agagCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAAc                                                                          >  2:2888919/1‑93 (MQ=255)
                agagCCTGGGGGAGGTTCACCTTATGTTGTACGATAAATCCCTTGAGAGGCGTAACTGTGCTTTCGGCCTGATTGCCCACATAGAAGGCGAAc                                                                          >  2:4113703/1‑93 (MQ=255)
                 gagCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAAcc                                                                         >  2:423894/1‑93 (MQ=255)
                  agCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCt                                                                        <  1:462576/93‑1 (MQ=255)
                  agCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGCCCACATAGAAGGCGAACCt                                                                        >  2:4026827/1‑93 (MQ=255)
                   gCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTa                                                                       <  2:1577356/93‑1 (MQ=255)
                   gCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTa                                                                       >  1:3557787/1‑93 (MQ=255)
                   gCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAAACTa                                                                       >  2:1918624/1‑93 (MQ=255)
                   gCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTGCGGCCTGATTGCCCACATAGAAGGCGAACCTa                                                                       <  2:2142639/93‑1 (MQ=255)
                    ccTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAg                                                                      >  1:636710/1‑93 (MQ=255)
                     cTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGc                                                                     >  1:3467692/1‑93 (MQ=255)
                        ggggAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCcac                                                                  <  2:1997620/93‑1 (MQ=255)
                        ggggAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCcac                                                                  <  1:2489022/93‑1 (MQ=255)
                        ggggAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACAAAGAAGGCGAACCTAGCcac                                                                  >  2:372502/1‑93 (MQ=255)
                         gggAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCcaca                                                                 >  2:708183/1‑93 (MQ=255)
                          ggAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACaa                                                                <  1:800485/93‑1 (MQ=255)
                            aGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAgg                                                              >  1:3069594/1‑93 (MQ=255)
                             ggTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGt                                                             >  2:534857/1‑93 (MQ=255)
                              gTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTa                                                            >  1:503095/1‑93 (MQ=255)
                               ttCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAg                                                           <  1:726773/93‑1 (MQ=255)
                                tCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGt                                                          >  2:834512/1‑93 (MQ=255)
                                tCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGt                                                          >  2:1442118/1‑93 (MQ=255)
                                tCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGt                                                          >  1:1944026/1‑93 (MQ=255)
                                 cACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTg                                                         >  1:881294/1‑93 (MQ=255)
                                 cACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTg                                                         >  2:615761/1‑93 (MQ=255)
                                  aCGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGc                                                        <  2:4655482/93‑1 (MQ=255)
                                  aCGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGc                                                        <  2:1125489/93‑1 (MQ=255)
                                  aCGATATGTTGTACGATAAATACCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGATAGTGc                                                        >  1:3140970/1‑93 (MQ=255)
                                   cGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGATGGCGAACCTAGCCACAAGGTAGTGCg                                                       >  1:2045305/1‑93 (MQ=255)
                                    gATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGt                                                      >  2:3785509/1‑93 (MQ=255)
                                    gATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGt                                                      >  1:3167217/1‑93 (MQ=255)
                                      tatGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTAc                                                    <  1:4717215/93‑1 (MQ=255)
                                       atGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACt                                                   <  2:1188977/93‑1 (MQ=255)
                                       atGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACt                                                   >  2:2432902/1‑93 (MQ=255)
                                       atGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACt                                                   <  2:3137356/93‑1 (MQ=255)
                                        tgttgtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTg                                                  >  2:183059/1‑93 (MQ=255)
                                         gttgtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGc                                                 <  1:405810/93‑1 (MQ=255)
                                         gttgtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGc                                                 >  1:821760/1‑93 (MQ=255)
                                         gttgtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGc                                                 >  1:1520294/1‑93 (MQ=255)
                                           tgtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCaa                                               >  2:4167724/1‑93 (MQ=255)
                                            gtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAAt                                              <  2:3396698/93‑1 (MQ=255)
                                            gtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAAt                                              >  2:1006236/1‑93 (MQ=255)
                                            gtACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAAt                                              <  2:1177951/93‑1 (MQ=255)
                                             tACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTGGCCACAAGGTAGTGCGTACTGCAATa                                             <  1:936925/93‑1 (MQ=255)
                                              aCGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATac                                            <  1:99464/93‑1 (MQ=255)
                                               cGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATaca                                           >  1:3242733/1‑93 (MQ=255)
                                               cGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATaca                                           <  2:3017100/93‑1 (MQ=255)
                                               cGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATaca                                           >  1:1730272/1‑93 (MQ=255)
                                                gATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATacac                                          >  1:2885117/1‑93 (MQ=255)
                                                gATAAATCCCTTGAGAGGGATAACTGTGGTTTAGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATacac                                          <  1:4672758/93‑1 (MQ=255)
                                                 aTAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATacacc                                         >  2:2320326/1‑92 (MQ=255)
                                                  tAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGc                                        >  2:3605016/1‑93 (MQ=255)
                                                  tAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGc                                        >  1:1783686/1‑93 (MQ=255)
                                                  tAAACCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACATAGCCACAAGGTAGCGCGTACTGCAATACACGc                                        >  2:2411236/1‑93 (MQ=255)
                                                    aaTCCCTTGAGATGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCAc                                      >  1:2291632/1‑93 (MQ=255)
                                                    aaTCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCAc                                      >  1:695022/1‑93 (MQ=255)
                                                     aTCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACt                                     <  2:3072435/93‑1 (MQ=255)
                                                     aTCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACt                                     >  2:2313701/1‑93 (MQ=255)
                                                      tCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTg                                    >  2:2207561/1‑93 (MQ=255)
                                                         cTTGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGcc                                 >  1:2437248/1‑93 (MQ=255)
                                                          ttGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGccc                                >  1:567682/1‑93 (MQ=255)
                                                           tGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCg                               <  1:1435829/93‑1 (MQ=255)
                                                           tGAGAGGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCg                               >  1:206755/1‑93 (MQ=255)
                                                             agagGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCa                             <  2:1509488/93‑1 (MQ=255)
                                                              gagGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCAt                            >  1:4144908/1‑93 (MQ=255)
                                                              gagGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCAt                            >  2:4403822/1‑93 (MQ=255)
                                                               agGGATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATg                           >  2:3223371/1‑93 (MQ=255)
                                                                gggATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgc                          >  1:1494671/1‑93 (MQ=255)
                                                                 ggATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgca                         <  2:3335062/93‑1 (MQ=255)
                                                                 ggATAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgca                         >  2:639426/1‑93 (MQ=255)
                                                                   aTAACTGTGTTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgcagc                       <  2:1749724/93‑1 (MQ=255)
                                                                   aTAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgcagc                       <  1:2569694/93‑1 (MQ=255)
                                                                   aTAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgcagc                       <  1:3497087/93‑1 (MQ=255)
                                                                    tAACTGTGGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATgcagcc                      >  1:3825512/1‑92 (MQ=255)
                                                                        tgtgGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGt                  <  2:2896633/93‑1 (MQ=255)
                                                                         gtgGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTg                 >  2:1611491/1‑93 (MQ=255)
                                                                         gtgGTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGGGCGGACTTGAATACACGCACTGGCCCGCATGCAGCATCGTg                 >  2:2622136/1‑93 (MQ=255)
                                                                            gTTTCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGgcg              >  1:3943107/1‑93 (MQ=255)
                                                                             tttCGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGgcgc             >  1:1221582/1‑93 (MQ=255)
                                                                                cGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGAt          >  2:3435781/1‑93 (MQ=255)
                                                                                cGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGAt          >  2:1281745/1‑93 (MQ=255)
                                                                                cGGCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGAt          <  1:1426899/93‑1 (MQ=255)
                                                                                 ggCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGAtt         <  2:4280601/93‑1 (MQ=255)
                                                                                  gCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATtc        >  1:4297710/1‑93 (MQ=255)
                                                                                  gCCTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGAGGCGCGATtc        >  1:1602888/1‑93 (MQ=255)
                                                                                   ccTGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATtct       <  1:4109210/93‑1 (MQ=255)
                                                                                     tGATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATTCTCg     >  1:4172421/1‑93 (MQ=255)
                                                                                      gATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATTCTCGc    >  1:2916145/1‑93 (MQ=255)
                                                                                      gATTGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATTCTCGc    <  2:2599443/93‑1 (MQ=255)
                                                                                        tCGCCCACATAGAAGGCGAACTTTGCCACAAGGTAGTGCGTACTGCAATACACACACTGGCCCGCATGCAGCACCGTGGCGAGAGTCTCGCCg  >  1:2757657/1‑93 (MQ=255)
                                                                                         |                                                                                           
GAAGTGGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGCCCACATAGAAGGCGAACCTAGCCACAAGGTAGTGCGTACTGCAATACACGCACTGGCCCGCATGCAGCACCGTGGCGCGATTCTCGCCG  >  NC_000913/3352321‑3352501

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