Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,390,710 T→C V312V (GTA→GTG yjeS ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,7100TC99.8% 1250.0 / ‑5.8 395V312V (GTA→GTGyjeShypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (183/211);  total (183/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 403 total aligned reads displayed.

CGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGG  >  NC_000913/4390620‑4390798
                                                                                          |                                                                                        
cgatATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCAc                                                                                        <  2:7634423/90‑1 (MQ=255)
cGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCAc                                                                                        <  1:2953722/93‑1 (MQ=255)
   tATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGc                                                                                     >  1:4357679/1‑93 (MQ=255)
      gtgCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGat                                                                                  >  1:107727/1‑93 (MQ=255)
      gtgCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGata                                                                                 >  1:4790273/1‑93 (MQ=255)
       tgCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGata                                                                                 <  2:7574351/93‑1 (MQ=255)
       tgCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGata                                                                                 <  1:3069935/93‑1 (MQ=255)
        gCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGatat                                                                                >  2:881181/1‑93 (MQ=255)
         cTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATAtt                                                                               >  1:6861885/1‑93 (MQ=255)
         cTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATAGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATAtt                                                                               >  1:6566938/1‑93 (MQ=255)
          tCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTa                                                                              >  1:4248086/1‑93 (MQ=255)
           cATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTAc                                                                             <  1:4632595/93‑1 (MQ=255)
           cATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTAc                                                                             <  1:6488278/93‑1 (MQ=255)
            aTCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACg                                                                            >  1:5274042/1‑93 (MQ=255)
            aTCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACg                                                                            <  1:6371018/93‑1 (MQ=255)
             tCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGc                                                                           <  2:7235722/93‑1 (MQ=255)
             tCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGc                                                                           >  1:4159064/1‑93 (MQ=255)
              cGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCa                                                                          <  2:6510392/93‑1 (MQ=255)
              cGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCa                                                                          <  1:569671/93‑1 (MQ=255)
               gagaAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAg                                                                         >  2:4067647/1‑93 (MQ=255)
                agaAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGACATTGCCTAATGCCACGGCGATATTACGCAGc                                                                        <  2:2169777/93‑1 (MQ=255)
                 gaAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGcc                                                                       <  1:3819040/93‑1 (MQ=255)
                  aaGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCa                                                                      >  1:3748247/1‑93 (MQ=255)
                   aGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCaa                                                                     >  2:2115113/1‑93 (MQ=255)
                   aGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCaa                                                                     <  2:189615/93‑1 (MQ=255)
                     ttggTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTGTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACg                                                                   <  1:2501446/91‑1 (MQ=255)
                     tGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACg                                                                   <  1:468816/93‑1 (MQ=255)
                      gggTGCTCTCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGc                                                                  <  1:2085390/93‑1 (MQ=255)
                       ggTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCa                                                                 <  1:7765360/93‑1 (MQ=255)
                         tGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAgg                                                               >  1:1135977/1‑93 (MQ=255)
                         tGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAgg                                                               <  1:1284336/93‑1 (MQ=255)
                         ggcacACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAgg                                                               <  2:4178087/89‑1 (MQ=255)
                          gCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAggg                                                              >  2:1660824/1‑92 (MQ=255)
                              aCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGAcc                                                          <  1:7083213/93‑1 (MQ=255)
                              aCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGAcc                                                          >  2:6227603/1‑93 (MQ=255)
                              aCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGAcc                                                          >  2:459907/1‑93 (MQ=255)
                               ccTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACTCAGCCAACGCAGGTGACCa                                                         >  1:938116/1‑93 (MQ=255)
                                cTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCaa                                                        <  2:2229512/93‑1 (MQ=255)
                                 tttACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAAt                                                       >  1:7159287/1‑93 (MQ=255)
                                  ttACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATa                                                      <  1:5475351/93‑1 (MQ=255)
                                   tACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATac                                                     <  1:156499/93‑1 (MQ=255)
                                   tACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATaa                                                     <  1:7925271/93‑2 (MQ=255)
                                      gACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATacgac                                                  >  1:3683750/1‑93 (MQ=255)
                                       aCTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATacgacg                                                 >  1:1196517/1‑93 (MQ=255)
                                        cTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATAcgacga                                                <  2:6822387/93‑1 (MQ=255)
                                        cTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATAcgacga                                                <  2:2409702/93‑1 (MQ=255)
                                         tttCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGaa                                               <  2:6039474/93‑1 (MQ=255)
                                         tttCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGaa                                               <  2:4953810/93‑1 (MQ=255)
                                         tttCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAAAGCAGGTGACCAATACGACGaa                                               <  1:1376135/93‑1 (MQ=255)
                                         cttCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGaa                                               <  2:2488606/92‑1 (MQ=255)
                                          ttCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAAt                                              <  1:2652510/93‑1 (MQ=255)
                                              aGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGcc                                          <  2:3097307/93‑1 (MQ=255)
                                                 gcTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGAt                                       <  1:7444704/93‑1 (MQ=255)
                                                  cTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATc                                      >  2:6882190/1‑93 (MQ=255)
                                                   tGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATcc                                     >  1:7613711/1‑93 (MQ=255)
                                                    gTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCt                                    >  1:1003689/1‑93 (MQ=255)
                                                    gTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCt                                    <  1:7410470/93‑1 (MQ=255)
                                                     tCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTAACCAATACGACGAATCGCCGATCCtt                                   >  1:5808003/1‑93 (MQ=255)
                                                      cAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTc                                  >  2:5810314/1‑93 (MQ=255)
                                                      cAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTc                                  <  2:271658/93‑1 (MQ=255)
                                                       aaaaTCGTTTCATCCCAAGGGGCATTTCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTcc                                 <  1:5960372/93‑1 (MQ=255)
                                                       aaaaTCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTcc                                 <  2:1424858/93‑1 (MQ=255)
                                                       aaaaTCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTcc                                 <  2:354461/93‑1 (MQ=255)
                                                       aaaaTCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTcc                                 >  2:7381282/1‑93 (MQ=255)
                                                           tCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGa                             >  1:254428/1‑93 (MQ=255)
                                                           tCGTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGa                             >  2:5980405/1‑93 (MQ=255)
                                                             gTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACt                           >  1:1718746/1‑93 (MQ=255)
                                                             gTTTCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACt                           >  2:6444684/1‑93 (MQ=255)
                                                               ttCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACTAATCGCCGATCCTTCCGTGACttt                         <  1:7520713/93‑1 (MQ=255)
                                                               ttCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACttt                         >  2:510376/1‑93 (MQ=255)
                                                               ttCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACttt                         <  2:5191727/93‑1 (MQ=255)
                                                               ttCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACttt                         <  2:4315697/93‑1 (MQ=255)
                                                               ttCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACttt                         >  1:4919457/1‑93 (MQ=255)
                                                                tCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACtttt                        >  1:2598883/1‑93 (MQ=255)
                                                                tCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACtttt                        >  2:6348607/1‑93 (MQ=255)
                                                                tCATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACtttt                        >  2:3280384/1‑93 (MQ=255)
                                                                 cATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTa                       <  2:7446697/93‑1 (MQ=255)
                                                                 cATCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTa                       >  2:1110648/1‑93 (MQ=255)
                                                                  aTCCCAAGGGGCATTGCCTAATGCCACGGCGATATTATGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTaa                      >  1:1892304/1‑93 (MQ=255)
                                                                  aTCCCAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTaa                      >  2:5721988/1‑93 (MQ=255)
                                                                     ccAAGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTaaaaa                   <  1:235371/93‑1 (MQ=255)
                                                                       aaGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACt                 >  2:7932733/1‑93 (MQ=255)
                                                                       aaGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACt                 <  2:7783699/93‑1 (MQ=255)
                                                                        aGGGGCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACtt                >  2:544763/1‑93 (MQ=255)
                                                                           ggCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtctc             >  1:3327414/1‑93 (MQ=255)
                                                                           ggCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtctc             <  2:1968864/93‑1 (MQ=255)
                                                                            gCATTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTctct            <  1:5063824/93‑1 (MQ=255)
                                                                              cttGCCTAATGCCACGGCAATATTACGCAGCCAACGCATGTGAGCAATCCAAAGAATCGCCGATCCTTCCGTGTCTTTTAAAAACTTCTCTTc          <  2:384406/92‑1 (MQ=255)
                                                                              aTTGTCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTAACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTc          <  2:2888007/93‑1 (MQ=255)
                                                                              aTTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTc          <  1:390949/93‑1 (MQ=255)
                                                                              aTTGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTc          >  2:153407/1‑93 (MQ=255)
                                                                               ttGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCg         >  2:6619965/1‑93 (MQ=255)
                                                                                tGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGc        >  2:5019599/1‑93 (MQ=255)
                                                                                tGCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGc        <  1:5401806/93‑1 (MQ=255)
                                                                                 gCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCt       >  1:2848469/1‑93 (MQ=255)
                                                                                 gCCTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCt       <  1:4577431/93‑1 (MQ=255)
                                                                                  ccTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTc      >  2:1192906/1‑93 (MQ=255)
                                                                                  ccTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTc      <  1:3526094/93‑1 (MQ=255)
                                                                                  ccTAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTc      >  2:3759928/1‑93 (MQ=255)
                                                                                    tAATGCCACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCa    >  1:2251792/1‑93 (MQ=255)
                                                                                          |                                                                                        
CGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGG  >  NC_000913/4390620‑4390798

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