Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,361,104 T→C 41.6% intergenic (+488/‑72) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,361,1040TC41.6% 68.2 / 148.7 135intergenic (+488/‑72)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base C (28/28);  total (68/67)
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 137 total aligned reads displayed.

CACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACA  >  NC_000913/3361005‑3361202
                                                                                                   |                                                                                                  
cACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACAtt                                                                                                    >  1:1799749/1‑100 (MQ=255)
  ccAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATGca                                                                                                  >  2:649041/1‑100 (MQ=255)
   cAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCcac                                                                                                 <  1:1167341/100‑1 (MQ=255)
     gCACCTTATCATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCcaca                                                                                                <  1:2893512/99‑1 (MQ=255)
       ccgatcttattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAAc                                                                                             <  1:1196807/94‑1 (MQ=255)
       agattattaTTCACATGAACGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAAc                                                                                             >  1:452565/4‑100 (MQ=255)
        ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACa                                                                                            >  2:128869/1‑100 (MQ=255)
        ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACa                                                                                            <  2:1003304/100‑1 (MQ=255)
         cTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTCTTTCTGAACTAGTTATTTATTTACATTCACAAACAg                                                                                           >  2:4036171/1‑100 (MQ=255)
           tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGat                                                                                         >  2:1720316/1‑100 (MQ=255)
           tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGat                                                                                         >  2:3535722/1‑100 (MQ=255)
            attattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGata                                                                                        >  2:257349/1‑100 (MQ=255)
             ttattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatat                                                                                       <  2:1816310/100‑1 (MQ=255)
             ttattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatat                                                                                       <  1:2719926/100‑1 (MQ=255)
             cgatctacaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatat                                                                                       <  2:3050828/94‑1 (MQ=255)
                ttCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatatata                                                                                    >  1:2308917/1‑100 (MQ=255)
                  cacaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGAtatatata                                                                                  >  2:145656/1‑100 (MQ=255)
                   acaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCTGCAAGTTGGCTATTTCTGAACGAGTTATTTAGTTACATCCACAAACAGATATATATAc                                                                                 <  1:1498354/100‑1 (MQ=255)
                   acaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATAc                                                                                 >  2:446061/1‑100 (MQ=255)
                    caTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACt                                                                                >  1:2423269/1‑100 (MQ=255)
                    caTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACt                                                                                >  1:1262616/1‑100 (MQ=255)
                    caTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACt                                                                                >  2:1793809/1‑100 (MQ=255)
                      tGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAGGTTGGCTATTTTTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGc                                                                              >  1:4332165/1‑100 (MQ=255)
                      tGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGc                                                                              >  2:2801387/1‑100 (MQ=255)
                           gTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTccg                                                                         >  1:55988/1‑98 (MQ=255)
                           gTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAg                                                                         >  1:3690901/1‑100 (MQ=255)
                           gTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAg                                                                         <  1:3563501/100‑1 (MQ=255)
                            tAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGc                                                                        <  1:2317700/100‑1 (MQ=255)
                            tAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGCTATATATACTGCTTCAGc                                                                        >  1:3548868/1‑100 (MQ=255)
                                tataAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCtttt                                                                    <  2:2425072/100‑1 (MQ=255)
                                tataAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCtttt                                                                    <  1:1403057/100‑1 (MQ=255)
                                 ataAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCttttt                                                                   >  2:1941061/1‑100 (MQ=255)
                                  taAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTtttta                                                                  <  1:47046/100‑1 (MQ=255)
                                    agaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatt                                                                >  2:2986463/3‑100 (MQ=255)
                                    aaaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttatt                                                                <  2:496280/100‑1 (MQ=255)
                                    aaaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATCTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttatt                                                                <  2:1634341/100‑1 (MQ=255)
                                    aaaTTCATTAAATAAATACATCCCAGCAAGTTGACTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatt                                                                <  1:830603/100‑1 (MQ=255)
                                     aaTTCATTAAATAAATACATCCCAGCAAGTTGGCTCTTTCTGAACGAGTTATTTATTTACATCAACCAACAGATATATATACTGCTTCAGCTttttattt                                                               >  2:364119/1‑100 (MQ=255)
                                     aaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttattt                                                               <  2:4360509/100‑1 (MQ=255)
                                      cgatcttAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttatttt                                                              <  2:3251383/95‑1 (MQ=255)
                                         cATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaa                                                           <  2:2513137/100‑1 (MQ=255)
                                         cATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaa                                                           >  1:1742959/1‑100 (MQ=255)
                                         aaTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaa                                                           >  1:1941118/2‑100 (MQ=255)
                                          aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                          <  2:3719399/100‑1 (MQ=255)
                                          aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGTTATATATACTGCTTCAGCTTTTTATTTTaaaa                                                          <  2:2939308/100‑1 (MQ=255)
                                          aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                          >  2:2668571/1‑100 (MQ=255)
                                           ttAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaa                                                         <  2:4224654/100‑1 (MQ=255)
                                            taaataaaTACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttattttagatcg                                                        >  1:3251365/1‑95 (MQ=255)
                                            taaataaaTACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaaa                                                        >  2:2808354/1‑100 (MQ=255)
                                             aaataaatACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaat                                                       >  1:1097359/1‑100 (MQ=255)
                                             aaataaatACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaat                                                       >  2:3369702/1‑100 (MQ=255)
                                              aataaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCTGCTTTTTATTTTAAAaaata                                                      >  1:4043033/1‑100 (MQ=255)
                                                taaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTCCAGCTTTTTATTTTAAAaaataaa                                                    >  2:2058015/1‑100 (MQ=255)
                                                taaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATATACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataaa                                                    >  2:2544992/1‑100 (MQ=255)
                                                    taCATCCCAGCAAGTTGGCTATTTCTTAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTg                                                <  1:2766889/100‑1 (MQ=255)
                                                    taCATCCCAGCAAGTTAGCTATTGCTGAACGAGTTATTTATTTACATTCACCAACAAAAATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAAGTTg                                                <  1:3301675/100‑1 (MQ=255)
                                                       aTCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGtta                                             <  2:345850/100‑1 (MQ=255)
                                                        tCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTCACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGt                                               >  1:2713057/1‑97 (MQ=255)
                                                         cccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTTttatt                                           >  1:2154903/1‑100 (MQ=255)
                                                         cccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatt                                           >  2:475170/1‑100 (MQ=255)
                                                          ccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatta                                          >  2:1876548/1‑100 (MQ=255)
                                                          ccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAATATAAATTTGttatta                                          >  1:4161738/1‑100 (MQ=255)
                                                          ccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatta                                          >  1:127133/1‑100 (MQ=255)
                                                           cAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTaa                                         <  2:1543141/100‑1 (MQ=255)
                                                           cAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTaa                                         <  2:3782210/100‑1 (MQ=255)
                                                             gCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAgg                                       <  1:2591604/100‑1 (MQ=255)
                                                              cAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGa                                      <  2:311758/100‑1 (MQ=255)
                                                               aaGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGat                                     <  1:2544979/100‑1 (MQ=255)
                                                                 gTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGatat                                   >  2:887821/1‑100 (MQ=255)
                                                                 gTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAAAAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGatat                                   >  1:34145/1‑100 (MQ=255)
                                                                      cTATTTCTGAACTAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGtt                                <  1:3222130/98‑1 (MQ=255)
                                                                      cTATTTCTGAACTAGTTATTTATTTACATTCACAAACAAATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGtt                                <  1:3221966/98‑1 (MQ=255)
                                                                      cTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTcn                              <  1:1541953/100‑2 (MQ=255)
                                                                      cTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                              <  2:3522655/100‑1 (MQ=255)
                                                                      cTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                              >  1:3649343/1‑100 (MQ=255)
                                                                      cTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                              <  1:3749511/100‑1 (MQ=255)
                                                                       tATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                              <  1:1828941/99‑1 (MQ=255)
                                                                        aTTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                            <  1:2735415/100‑1 (MQ=255)
                                                                        aTTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                            <  1:1365931/100‑1 (MQ=255)
                                                                         tttCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCatat                           >  1:81591/1‑100 (MQ=255)
                                                                           tCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGt                         >  2:4033933/1‑100 (MQ=255)
                                                                           tCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGt                         >  2:90360/1‑100 (MQ=255)
                                                                             tGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttt                       >  2:703428/1‑100 (MQ=255)
                                                                               aaCGAGTTATTTATTTACATCCACAAACAAAAATATATACTGCTTCAGCTTTTTATTTTGAAAAATAAATTTGTTATTAAGGATATGTTCATATGttt                       <  1:2714211/98‑1 (MQ=255)
                                                                                aCGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTACTAAGGATATGTTCATATGTTTTTc                    <  1:1365955/100‑1 (MQ=255)
                                                                                 cGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCa                   >  2:1407259/1‑100 (MQ=255)
                                                                                 cGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCa                   >  1:415825/1‑100 (MQ=255)
                                                                                  gAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATCTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaa                  >  2:398503/1‑100 (MQ=255)
                                                                                  gAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaa                  >  1:2933805/1‑100 (MQ=255)
                                                                                   aGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaa                 >  2:1667245/1‑100 (MQ=255)
                                                                                   aGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaa                 >  1:4354611/1‑100 (MQ=255)
                                                                                    gTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaa                <  2:514012/100‑1 (MQ=255)
                                                                                     ttatttatTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaaa               <  2:3507947/100‑1 (MQ=255)
                                                                                     ttatttatTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaaa               <  2:2375631/100‑1 (MQ=255)
                                                                                     ttatttatTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaa                <  2:2194298/99‑1 (MQ=255)
                                                                                      tatttattTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAg              <  1:2762220/100‑1 (MQ=255)
                                                                                              taCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGACCCTcac      <  1:2111890/100‑1 (MQ=255)
                                                                                               tcATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcaca     <  2:3709210/99‑1 (MQ=255)
                                                                                                cATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCacaa    >  1:3709184/1‑100 (MQ=255)
                                                                                                  ttCACAAACAGATATATATACTGCTTAAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCacaaca  >  1:1063996/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACA  >  NC_000913/3361005‑3361202

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: