Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,588 C→G R86T (AGG→ACG)  gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,5880CG99.7% 1024.9 / ‑4.6 295R86T (AGG→ACG) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (162/132);  minor base T (0/1);  total (162/133)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

Only 100 of 304 total aligned reads displayed.

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677
                                                                                          |                                                                                         
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         <  1:1578485/93‑1 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  2:3685648/1‑93 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  2:3993458/1‑93 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGAGTAACGTc                                                                                         >  2:2787110/1‑93 (MQ=255)
 ccTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                                        >  2:6213474/1‑93 (MQ=255)
 ccTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                                        <  2:5092845/93‑1 (MQ=255)
  ctctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGc                                                                                       <  2:5450061/93‑1 (MQ=255)
   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      >  1:4520660/1‑93 (MQ=255)
    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     >  1:1291513/1‑93 (MQ=255)
        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 <  1:5798911/93‑1 (MQ=255)
        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 >  2:5258162/1‑93 (MQ=255)
         ccGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                                                <  2:4888855/93‑1 (MQ=255)
         ccGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                                                >  1:4175942/1‑93 (MQ=255)
          cGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGa                                                                               <  1:535448/93‑1 (MQ=255)
             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  2:2406242/1‑93 (MQ=255)
             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  2:341253/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:2002833/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  1:3496404/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:2123138/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           <  2:4130417/93‑1 (MQ=255)
               aTCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcg                                                                          >  1:721935/1‑93 (MQ=255)
                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         >  2:3604942/1‑93 (MQ=255)
                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         >  1:4052653/1‑93 (MQ=255)
                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         <  2:5679112/93‑1 (MQ=255)
                 cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                                        <  2:2286393/93‑1 (MQ=255)
                  aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                                       >  2:1289149/1‑93 (MQ=255)
                  aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                                       >  2:1012427/1‑93 (MQ=255)
                  aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                                       >  1:2249490/1‑93 (MQ=255)
                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      <  1:3770828/93‑1 (MQ=255)
                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    <  2:2972086/93‑1 (MQ=255)
                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    <  2:6449211/93‑1 (MQ=255)
                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    <  1:975522/93‑1 (MQ=255)
                      gTAACGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                                   >  1:3991011/1‑93 (MQ=255)
                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                                                 <  1:4409479/93‑1 (MQ=255)
                          cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                                               >  2:3236753/1‑93 (MQ=255)
                          cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                                               <  1:6364879/93‑1 (MQ=255)
                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgag                                                              <  1:1749690/93‑1 (MQ=255)
                              cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgaggga                                                           >  1:1859823/1‑93 (MQ=255)
                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGATGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                                          <  1:5065496/93‑1 (MQ=255)
                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                                          <  1:856610/93‑1 (MQ=255)
                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                                          >  1:3348055/1‑93 (MQ=255)
                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  1:4904785/1‑93 (MQ=255)
                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                                      >  2:939282/1‑93 (MQ=255)
                                    agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAata                                                     >  1:5561082/1‑93 (MQ=255)
                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                                   >  2:3892121/1‑93 (MQ=255)
                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                                  >  2:5904493/1‑93 (MQ=255)
                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                                  <  1:5938622/93‑1 (MQ=255)
                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                                 <  2:6181219/93‑1 (MQ=255)
                                          tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                                               <  1:5004786/93‑1 (MQ=255)
                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                                              <  1:4395139/93‑1 (MQ=255)
                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  1:1183395/93‑1 (MQ=255)
                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGGGGGAGGAATATCTTGCAGACGt                                         >  2:4727549/1‑93 (MQ=255)
                                                cgcttGGCAACAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  1:1480818/90‑1 (MQ=255)
                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                        >  1:2194083/1‑93 (MQ=255)
                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       <  2:5206521/93‑1 (MQ=255)
                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       >  2:4578069/1‑93 (MQ=255)
                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAGCGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:4217641/93‑1 (MQ=255)
                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      >  1:5212316/1‑93 (MQ=255)
                                                   ttGGCAATAAAATCACCCGCATTGACCGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:829614/93‑1 (MQ=255)
                                                    tGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAg                                     >  2:2168743/1‑93 (MQ=255)
                                                    tGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAg                                     >  1:4607058/1‑93 (MQ=255)
                                                        aagaaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                 <  1:2468555/90‑1 (MQ=255)
                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               <  2:1499209/93‑1 (MQ=255)
                                                           aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                              <  2:4325674/93‑1 (MQ=255)
                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             >  2:1698787/1‑93 (MQ=255)
                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             >  2:1122429/1‑93 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  2:5613051/1‑93 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  1:3448273/1‑93 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  2:187239/1‑93 (MQ=255)
                                                               tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc                          <  1:685156/93‑1 (MQ=255)
                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       <  2:5776167/93‑1 (MQ=255)
                                                                  cccGCACTGCCGGCACGACGTCACGTCGCCGTAGCGCCACGGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACCGGCCCaaa                       <  1:2293538/93‑1 (MQ=255)
                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                      <  2:3159450/93‑1 (MQ=255)
                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                      <  1:3319659/93‑1 (MQ=255)
                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     >  1:904237/1‑93 (MQ=255)
                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     <  2:68941/93‑1 (MQ=255)
                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     <  2:501366/93‑1 (MQ=255)
                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     >  1:6288818/1‑93 (MQ=255)
                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    >  2:1462985/1‑93 (MQ=255)
                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  1:3755533/93‑1 (MQ=255)
                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   <  2:4264246/93‑1 (MQ=255)
                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTATCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  1:2742777/1‑93 (MQ=255)
                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  2:1609632/1‑93 (MQ=255)
                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  1:3008277/1‑93 (MQ=255)
                                                                         tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                <  1:2401349/93‑1 (MQ=255)
                                                                          gACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCa               >  1:478676/1‑93 (MQ=255)
                                                                             ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc            >  2:5403612/1‑93 (MQ=255)
                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           <  1:6237769/93‑1 (MQ=255)
                                                                                aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  2:1961026/1‑93 (MQ=255)
                                                                                aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  1:3075870/1‑93 (MQ=255)
                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:3189200/1‑93 (MQ=255)
                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:513399/1‑93 (MQ=255)
                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:6475420/1‑93 (MQ=255)
                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        <  2:3268255/93‑1 (MQ=255)
                                                                                     gTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    >  1:3288399/1‑93 (MQ=255)
                                                                                      tAACGTCGCCGTCGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   >  1:1053192/1‑93 (MQ=255)
                                                                                      tAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   >  1:2008821/1‑93 (MQ=255)
                                                                                      tAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   <  2:2549324/93‑1 (MQ=255)
                                                                                      gaaCGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   <  1:4970845/92‑1 (MQ=255)
                                                                                       aaCGTCGCCGTAGCGACAGCGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAtt  <  1:28306/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

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