Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,372,616 T→C I168T (ATC→ACC)  ycjP → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,372,6160TC99.2% 367.3 / ‑4.5 120I168T (ATC→ACC) ycjPputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (50/69);  total (50/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

Only 100 of 129 total aligned reads displayed.

TTTACGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTGGCGATGTCC  >  NC_000913/1372526‑1372704
                                                                                          |                                                                                        
tttACGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCAccc                                                                                        >  2:2873917/1‑93 (MQ=255)
tttACGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCAccc                                                                                        >  2:4654297/1‑93 (MQ=255)
    cGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGa                                                                                    <  2:1382512/93‑1 (MQ=255)
    cGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGa                                                                                    <  1:3141909/93‑1 (MQ=255)
     gACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGAt                                                                                   <  2:3117247/93‑1 (MQ=255)
     gACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGAt                                                                                   <  1:824932/93‑1 (MQ=255)
     gACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGAt                                                                                   >  2:1584609/1‑93 (MQ=255)
     gACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGAt                                                                                   >  2:2049782/1‑93 (MQ=255)
       cacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATga                                                                                 <  2:2453429/93‑1 (MQ=255)
       cacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATga                                                                                 <  2:2928247/93‑1 (MQ=255)
        acTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATgag                                                                                >  1:1575694/1‑93 (MQ=255)
        acTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATgag                                                                                <  1:692603/93‑1 (MQ=255)
        acTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATgag                                                                                <  2:208576/93‑1 (MQ=255)
         cTGAGATGGCGCCGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATgaga                                                                               <  2:3734030/93‑1 (MQ=255)
           gagaTGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATc                                                                             >  1:2623512/1‑93 (MQ=255)
                  cgcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGaagaag                                                                      <  2:394537/93‑1 (MQ=255)
                  cgcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGaagaag                                                                      <  1:3486029/93‑1 (MQ=255)
                  cgcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGaagaag                                                                      <  2:2686884/93‑1 (MQ=255)
                   gcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgc                                                                     >  2:3445396/1‑93 (MQ=255)
                   gcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgc                                                                     <  2:44725/93‑1 (MQ=255)
                      tatcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgcggc                                                                  <  1:1595692/92‑1 (MQ=255)
                       atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgcggcg                                                                 >  1:875464/1‑93 (MQ=255)
                       atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgcggcg                                                                 >  1:4225471/1‑93 (MQ=255)
                       atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAgcggcg                                                                 >  2:1942729/1‑93 (MQ=255)
                        tcatcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCga                                                                >  2:3942963/1‑93 (MQ=255)
                        tcatcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCga                                                                >  1:4079685/1‑93 (MQ=255)
                        tcatcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCga                                                                <  1:1187546/93‑1 (MQ=255)
                        tcatcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCga                                                                >  1:2490795/1‑93 (MQ=255)
                         catcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCgat                                                               <  2:4084181/93‑1 (MQ=255)
                           tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCgatga                                                             >  1:3090263/1‑93 (MQ=255)
                           tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCgatga                                                             <  2:1487287/93‑1 (MQ=255)
                           tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCgatga                                                             >  2:35461/1‑93 (MQ=255)
                             aCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGatgatg                                                           >  1:4687748/1‑93 (MQ=255)
                               cATGGTGACGCAACCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatgga                                                         <  2:3935152/93‑1 (MQ=255)
                               cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATgatgga                                                         <  1:4482988/93‑1 (MQ=255)
                               cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatgga                                                         <  1:2941037/93‑1 (MQ=255)
                               cATGGTGACGCAAACTCTGCCCACCGCCGTGTGTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatgga                                                         <  2:144283/93‑1 (MQ=255)
                                aTGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggat                                                        <  1:2806786/93‑1 (MQ=255)
                                aTGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggat                                                        <  2:2152801/93‑1 (MQ=255)
                                 tGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggatg                                                       >  2:3557758/1‑93 (MQ=255)
                                 tGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggatg                                                       <  2:2213912/93‑1 (MQ=255)
                                 tGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggatg                                                       <  2:4281014/93‑1 (MQ=255)
                                 tGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggatg                                                       <  2:2580245/93‑1 (MQ=255)
                                 tGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATgatggatg                                                       <  2:2717605/93‑1 (MQ=255)
                                  ggTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACGCCACCCCCGATGAGATCGAAGAAGCGGCGATGatggatgg                                                      <  2:3007978/93‑1 (MQ=255)
                                    tGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGct                                                    >  1:3006170/1‑92 (MQ=255)
                                    tGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGcc                                                    <  1:1778611/93‑1 (MQ=255)
                                       cGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCa                                                 >  2:4033714/1‑93 (MQ=255)
                                        gCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCaa                                                >  2:2629096/1‑93 (MQ=255)
                                         cAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAAc                                               >  2:126923/1‑93 (MQ=255)
                                         cAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAAc                                               <  1:2525067/93‑1 (MQ=255)
                                           aaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCg                                             >  2:1615300/1‑93 (MQ=255)
                                             ctctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTc                                           <  1:4167290/93‑1 (MQ=255)
                                             ctctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTc                                           >  2:3625974/1‑93 (MQ=255)
                                              tctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCt                                          >  1:294572/1‑93 (MQ=255)
                                              tctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCt                                          <  2:1388618/93‑1 (MQ=255)
                                                 gCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCa                                       <  1:4673322/93‑1 (MQ=255)
                                                  cccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaa                                      >  1:3657878/1‑93 (MQ=255)
                                                  cccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaa                                      <  1:2261709/93‑1 (MQ=255)
                                                   ccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaaa                                     >  1:2677586/1‑93 (MQ=255)
                                                   ccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaaa                                     >  2:3150236/1‑93 (MQ=255)
                                                   ccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaaa                                     <  1:4028526/93‑1 (MQ=255)
                                                   ccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCaaa                                     >  2:3467794/1‑93 (MQ=255)
                                                     aCCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAAtt                                   >  2:3424605/1‑93 (MQ=255)
                                                      ccgccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTa                                  <  1:647556/93‑1 (MQ=255)
                                                      ccgccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTa                                  <  2:4908/93‑1 (MQ=255)
                                                       cgccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTAt                                 <  1:1142948/93‑1 (MQ=255)
                                                        gccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATc                                <  1:3030023/93‑1 (MQ=255)
                                                        gccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATc                                <  2:3536937/93‑1 (MQ=255)
                                                         ccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCt                               <  1:3701579/93‑1 (MQ=255)
                                                         ccgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCt                               >  1:903529/1‑93 (MQ=255)
                                                          cgTGTTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCtt                              <  2:4527587/93‑1 (MQ=255)
                                                           gtgtTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTc                             <  1:1609027/93‑1 (MQ=255)
                                                             gtTTATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCg                           <  2:3906776/93‑1 (MQ=255)
                                                              tttATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGc                          <  2:900104/93‑1 (MQ=255)
                                                              tttATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGc                          >  1:3191856/1‑93 (MQ=255)
                                                               ttATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCa                         >  2:3238012/1‑93 (MQ=255)
                                                                tATGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCAt                        <  2:4140321/93‑1 (MQ=255)
                                                                 aTGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATc                       <  2:4465704/93‑1 (MQ=255)
                                                                 aTGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATc                       <  1:4386826/93‑1 (MQ=255)
                                                                 aTGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATc                       >  2:3160385/1‑93 (MQ=255)
                                                                  tGCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCa                      >  1:102099/1‑93 (MQ=255)
                                                                   gCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCAc                     >  1:2455601/1‑93 (MQ=255)
                                                                   gCTGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCAc                     <  2:3334699/93‑1 (MQ=255)
                                                                     tGAAAAGCTACTTCGACTCCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCg                   <  2:143735/93‑1 (MQ=255)
                                                                     tGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCg                   >  2:759651/1‑93 (MQ=255)
                                                                     tGAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCg                   >  2:869494/1‑93 (MQ=255)
                                                                      gAAAAGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGt                  >  1:3832750/1‑93 (MQ=255)
                                                                       aaaaGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTg                 <  1:4244909/93‑1 (MQ=255)
                                                                       aaaaGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTg                 >  2:1938626/1‑93 (MQ=255)
                                                                        aaaGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGc                <  1:125524/93‑1 (MQ=255)
                                                                         aaGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGcc               >  2:126677/1‑93 (MQ=255)
                                                                         aaGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGcc               <  1:781306/93‑1 (MQ=255)
                                                                          aGCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCg              >  2:3693364/1‑93 (MQ=255)
                                                                           gCTACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGt             <  1:4548814/93‑1 (MQ=255)
                                                                             tACTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTg           <  1:590371/93‑1 (MQ=255)
                                                                              aCTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTgg          <  1:904052/93‑1 (MQ=255)
                                                                              aCTTCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTgg          >  1:4164486/1‑93 (MQ=255)
                                                                                 tCGACACCACCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTGGCGa       <  2:2731054/93‑1 (MQ=255)
                                                                                  cGACACCGCCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTGGCGAt      <  2:514507/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTTACGACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTCTGCAAATTATCTTCCGCATCACCGTGCCGTTGGCGATGTCC  >  NC_000913/1372526‑1372704

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