Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,636 T→C intergenic (+282/+233) b0501 → / ← b0502 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,6360TC99.8% 1806.6 / ‑5.4 558intergenic (+282/+233)b0501/b0502hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (240/317);  total (241/317)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 577 total aligned reads displayed.

TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725
                                                                                          |                                                                                         
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         >  2:772979/1‑93 (MQ=255)
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         >  2:1631308/1‑93 (MQ=255)
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         >  2:5315872/1‑93 (MQ=255)
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         <  1:8514914/93‑1 (MQ=255)
  agagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGttt                                                                                       >  1:9824005/1‑93 (MQ=255)
  agagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGttt                                                                                       >  1:4328818/1‑93 (MQ=255)
   gagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTa                                                                                      >  1:937565/1‑93 (MQ=255)
    agagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTaa                                                                                     <  2:6116465/93‑1 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    >  2:8612435/1‑93 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    >  1:8621561/1‑93 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    <  2:1844436/93‑1 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    <  1:3422211/93‑1 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    >  2:4780159/1‑93 (MQ=255)
      agTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtg                                                                                   <  1:1313310/93‑1 (MQ=255)
       gTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTCCGCGTTTAAttg                                                                                  >  2:8732220/1‑91 (MQ=255)
        tttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtg                                                                                 >  1:2258599/1‑93 (MQ=255)
        tttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtg                                                                                 <  2:3010291/93‑1 (MQ=255)
        tttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtg                                                                                 >  2:8848698/1‑93 (MQ=255)
         ttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGa                                                                                <  1:5737455/93‑1 (MQ=255)
         ttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGa                                                                                <  2:4187925/93‑1 (MQ=255)
          tACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtgcg                                                                               >  2:2800858/1‑91 (MQ=255)
          tACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAg                                                                               <  2:7569421/93‑1 (MQ=255)
            cATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTc                                                                             >  1:9964462/1‑93 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         <  2:7233754/93‑1 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         <  1:5371940/93‑1 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         <  1:6911511/93‑1 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         >  1:7360465/1‑93 (MQ=255)
                  aTTGCACTGGTCTTTGACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATc                                                                       <  2:9314535/93‑1 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      >  2:952506/1‑93 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      <  1:2081641/93‑1 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      <  2:3568278/93‑1 (MQ=255)
                       aCTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaa                                                                  >  2:2146675/1‑93 (MQ=255)
                       aCTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaa                                                                  >  1:7854573/1‑93 (MQ=255)
                       aCTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaa                                                                  <  2:9790307/93‑1 (MQ=255)
                        cTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaaa                                                                 <  2:4596622/93‑1 (MQ=255)
                        cTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaaa                                                                 >  1:6220350/1‑93 (MQ=255)
                         tGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaat                                                                >  1:9560544/1‑93 (MQ=255)
                          ggTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaata                                                               >  2:1542611/1‑93 (MQ=255)
                            tCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGCTTAATGTGAGTCCCTATCTGAaaataaa                                                             >  2:8149070/1‑93 (MQ=255)
                              ttttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAaataaata                                                           >  1:5035738/1‑93 (MQ=255)
                              ttttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAaataaata                                                           <  2:6328339/93‑1 (MQ=255)
                                 tACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATc                                                        >  2:5105909/1‑93 (MQ=255)
                                  aCGATATCTTACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATcc                                                       <  1:4017528/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  2:9011500/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  2:747248/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      >  1:3589324/1‑93 (MQ=255)
                                    gATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCtt                                                     >  2:6908889/1‑93 (MQ=255)
                                    gATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCtt                                                     <  2:447261/93‑1 (MQ=255)
                                     atatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTc                                                    <  1:6727320/93‑1 (MQ=255)
                                      tatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTcc                                                   >  1:2646515/1‑93 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 <  1:6451115/93‑1 (MQ=255)
                                         cTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGa                                                <  1:4149159/93‑1 (MQ=255)
                                           gACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAtt                                              <  2:1411664/93‑1 (MQ=255)
                                            aCATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTa                                             >  1:7048035/1‑93 (MQ=255)
                                              aTTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaa                                           >  1:3466336/1‑93 (MQ=255)
                                              aTTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaa                                           >  1:2994370/1‑93 (MQ=255)
                                              aTTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaa                                           <  2:4060049/93‑1 (MQ=255)
                                              aTTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaa                                           >  1:1575262/1‑93 (MQ=255)
                                                    tAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaat                                     <  2:919404/93‑1 (MQ=255)
                                                    tAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaat                                     <  1:5266543/93‑1 (MQ=255)
                                                     aaTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAAtt                                    <  2:4998268/93‑1 (MQ=255)
                                                      aTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTc                                   >  2:9435946/1‑93 (MQ=255)
                                                      aTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCAGGATTAAAATAAATTc                                   >  1:1635972/1‑93 (MQ=255)
                                                       tACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCt                                  <  2:5657791/93‑1 (MQ=255)
                                                        aCATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCtt                                 <  1:5835/93‑1 (MQ=255)
                                                        aCATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCtt                                 >  2:314961/1‑93 (MQ=255)
                                                         cATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTg                                >  1:1257472/1‑93 (MQ=255)
                                                         cATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTg                                <  1:324204/93‑1 (MQ=255)
                                                           tatTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGcc                              <  1:2807117/93‑1 (MQ=255)
                                                            atTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGcag                             <  1:1083428/93‑3 (MQ=255)
                                                            atTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCg                             <  2:9079880/93‑1 (MQ=255)
                                                             ttCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCgg                            >  1:8060847/1‑93 (MQ=255)
                                                                aGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaa                         <  1:893943/93‑1 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        >  1:8271941/1‑93 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        >  2:6617637/1‑93 (MQ=255)
                                                                   cATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaaga                      <  1:3111440/93‑1 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     <  1:9083641/93‑1 (MQ=255)
                                                                     tGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaaa                    >  2:7170837/1‑93 (MQ=255)
                                                                      gCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGaaagaaag                   >  2:2037268/1‑93 (MQ=255)
                                                                       cATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAaagaaaga                  <  2:1069324/93‑1 (MQ=255)
                                                                       aattaataaAATGTTACTGCTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAagag                 <  1:9237507/92‑1 (MQ=255)
                                                                        attaattaAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAagag                 <  1:2465413/93‑1 (MQ=255)
                                                                         ttaattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAgg                <  2:6002436/93‑1 (MQ=255)
                                                                           aattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaa              >  2:399202/1‑93 (MQ=255)
                                                                           aattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaa              >  2:8354495/1‑93 (MQ=255)
                                                                           aattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaa              >  1:5597415/1‑93 (MQ=255)
                                                                            attaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaaa             <  2:237335/93‑1 (MQ=255)
                                                                             ttaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAAt            <  1:7170591/93‑1 (MQ=255)
                                                                              taaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATa           <  2:380813/93‑1 (MQ=255)
                                                                              taaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATa           <  2:7081890/93‑1 (MQ=255)
                                                                               aaaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaa          <  2:5585742/93‑1 (MQ=255)
                                                                               aaaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaa          >  1:4588637/1‑93 (MQ=255)
                                                                                aaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaaa         <  2:1155085/93‑1 (MQ=255)
                                                                                aaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaaa         >  2:6564441/1‑93 (MQ=255)
                                                                                  aTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAcc       <  2:330730/93‑1 (MQ=255)
                                                                                   tGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCa      <  1:7560234/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     >  2:1292981/1‑93 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     <  2:2600831/93‑1 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    <  1:3733589/93‑1 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    <  1:1016406/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725

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