Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,571,816 T→G D77A (GAT→GCT)  gadB ← glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,571,8160TG98.2% 769.4 / 2.2 221D77A (GAT→GCT) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base G (118/98);  total (119/102)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 221 total aligned reads displayed.

ATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATA  >  NC_000913/1571719‑1571911
                                                                                                 |                                                                                               
  gCCACAGATAGGCAAACATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGCTATTCTTCTTTGTAGATCCAGTTTTTGTTAATGGATAAAGccct                                                                                             >  2:528922/1‑98 (MQ=14)
   ccACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATc                                                                                            <  1:5049241/100‑1 (MQ=255)
   ccACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATc                                                                                            <  1:1684904/100‑1 (MQ=255)
   ccACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAACCCATc                                                                                            <  2:4843564/100‑1 (MQ=255)
    nacaGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCa                                                                                           >  1:1547880/2‑100 (MQ=255)
    cacaGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCa                                                                                           >  1:2395721/1‑100 (MQ=255)
    cacaGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCa                                                                                           >  2:264291/1‑100 (MQ=255)
    cacaGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCa                                                                                           >  1:2811778/1‑100 (MQ=255)
        gATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAAttt                                                                                       >  1:3155911/1‑100 (MQ=255)
        gATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAAttt                                                                                       >  1:4117252/1‑100 (MQ=255)
          tCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgt                                                                                     <  1:290955/100‑1 (MQ=255)
           cGGCAACCATATTTAGGCAACGCAGGTCGACGGCTGCGGCTTGCGCATATTCTTCTTGGTCGATCCAGTTTTTGTTGATGGATAAAGCCAGCAATTtgtg                                                                                    >  1:4848318/1‑100 (MQ=14)
           cGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgtg                                                                                    <  1:182837/100‑1 (MQ=255)
           cGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgtg                                                                                    <  1:1425971/100‑1 (MQ=255)
            ggCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGTATAAAGCCATCAATTTGTgg                                                                                   <  1:3484761/100‑1 (MQ=255)
             gCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGa                                                                                  <  1:4562288/100‑1 (MQ=255)
              cAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGAc                                                                                 >  2:2955822/1‑100 (MQ=255)
              cAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGAc                                                                                 >  1:3566579/1‑100 (MQ=255)
               aaCCATATTTACGCAACGCAGGTCGATGGCTGCGGCTTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACa                                                                                <  2:4637469/100‑1 (MQ=255)
                   atatTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACAtttt                                                                            <  2:5187373/100‑1 (MQ=255)
                     atTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcg                                                                          >  1:202028/1‑100 (MQ=255)
                      tttACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcgt                                                                         <  2:5410653/100‑1 (MQ=255)
                      tttACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcgt                                                                         <  2:344876/100‑1 (MQ=255)
                         acgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTcgtcgt                                                                      >  1:4363391/1‑100 (MQ=255)
                          cgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCgtcgtc                                                                     <  2:4522736/100‑1 (MQ=255)
                          cgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCgtcgtc                                                                     <  1:3713059/100‑1 (MQ=255)
                          cgcaacgcaGGTCGATGACTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCgtcgtc                                                                     >  1:1883055/1‑100 (MQ=255)
                            caacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTccc                                                                   <  1:2669969/100‑1 (MQ=255)
                            caacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTccc                                                                   >  1:4724006/1‑100 (MQ=255)
                            caacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTccc                                                                   <  2:3721356/100‑1 (MQ=255)
                             aacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCa                                                                  >  1:1590711/1‑100 (MQ=255)
                             aacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCa                                                                  <  1:1178951/100‑1 (MQ=255)
                               cgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGCTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAgg                                                                <  1:5297656/100‑1 (MQ=255)
                               cgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAgg                                                                <  2:2200202/100‑1 (MQ=255)
                                   ggTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTGGATGCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg                                                            >  2:5158671/1‑100 (MQ=17)
                                   ggTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGTTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg                                                            >  2:4484037/1‑100 (MQ=255)
                                   ggTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg                                                            >  1:1079626/1‑100 (MQ=255)
                                     ttgATGGCTGCGGATTGCGGATCTTCTTCTTTGTCGATCCACTTTTTGTTACTGCATAAAGCCATCAATTTGTGGAAATTTTCGTCGTCCCAGGTCTGGc                                                          <  2:1200666/98‑1 (MQ=17)
                                     tCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGc                                                          <  2:392734/100‑1 (MQ=255)
                                     tCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGc                                                          <  1:1769280/100‑1 (MQ=255)
                                      tgATGGCTGCGGATTGCGGATATTCTTCTTGGTCGATCCAGTTTTTGTTAATGGATAAAGCCATGAATTTGTGGACATTTTCGTCGGCCCAGGTCTGGCa                                                         <  2:786489/99‑1 (MQ=17)
                                      cGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCa                                                         <  1:2993779/100‑1 (MQ=255)
                                      cGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCa                                                         >  2:80826/1‑100 (MQ=255)
                                         tGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGaa                                                      <  2:1794572/100‑1 (MQ=255)
                                           tcTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAg                                                    <  2:4398804/99‑1 (MQ=255)
                                           gCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAg                                                    <  1:890080/100‑1 (MQ=255)
                                           gCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAg                                                    <  1:5074634/100‑1 (MQ=255)
                                           gCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAg                                                    >  1:1863133/1‑100 (MQ=255)
                                             tgcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTg                                                  >  2:3868511/1‑100 (MQ=255)
                                                ngattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGcc                                               >  1:1146201/2‑100 (MQ=255)
                                                ggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGcn                                               <  2:2395429/100‑2 (MQ=255)
                                                 gattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCa                                              <  1:5508849/100‑1 (MQ=255)
                                                  attgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAg                                             <  1:245885/100‑1 (MQ=255)
                                                  attgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAg                                             <  1:3963245/100‑1 (MQ=255)
                                                  attgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAg                                             <  2:3834073/100‑1 (MQ=255)
                                                   tcgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATGAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAgg                                            <  2:5473393/98‑1 (MQ=255)
                                                    tgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGt                                           <  2:4103570/100‑1 (MQ=255)
                                                     gcggatATTCTTCTTTGTCGATCCAGTTTTTTTTTTTTTTTTTAGCCATCAATTTGTGGACATTTTCGTCGTCCGAGGTCTGGCAGAAAGTGGTCAGCtt                                          >  2:1710850/1‑100 (MQ=14)
                                                      cggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTc                                         >  1:980883/1‑100 (MQ=255)
                                                       ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCt                                        <  2:1873617/100‑1 (MQ=255)
                                                       ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCt                                        >  1:3794708/1‑100 (MQ=255)
                                                        gatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTg                                       >  1:3029330/1‑100 (MQ=255)
                                                          tatTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGAc                                     <  1:1488832/100‑1 (MQ=255)
                                                          tatTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGAc                                     <  2:3572407/100‑1 (MQ=255)
                                                            ttcttcTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGa                                   >  2:688819/1‑100 (MQ=255)
                                                              cttcttTGTCGATCCAGTTTTTGTTAATGGATAAAGCCGTCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGc                                 >  1:4937222/1‑100 (MQ=255)
                                                              cttcttTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGc                                 >  1:5400727/1‑100 (MQ=255)
                                                                tcttTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGt                               <  1:4133258/100‑1 (MQ=255)
                                                                  tttGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTg                             <  2:3215761/100‑1 (MQ=255)
                                                                   ttGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGc                            >  1:145644/1‑100 (MQ=255)
                                                                   ttGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGc                            <  1:697349/100‑1 (MQ=255)
                                                                     gTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCa                          <  1:3449413/100‑1 (MQ=255)
                                                                       cGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATc                        >  2:903447/1‑100 (MQ=255)
                                                                       cGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATc                        >  2:3118689/1‑100 (MQ=255)
                                                                       cGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGGTGCCATc                        >  1:256911/1‑100 (MQ=255)
                                                                        gATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCa                       >  2:2173047/1‑100 (MQ=255)
                                                                        gATCCAGTTATTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCa                       >  2:5341431/1‑100 (MQ=255)
                                                                         aTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCaa                      >  1:1350989/1‑100 (MQ=255)
                                                                         aTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCaa                      <  1:3617561/100‑1 (MQ=255)
                                                                          tCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAg                     <  2:2760064/100‑1 (MQ=255)
                                                                          tCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAg                     >  2:751067/1‑100 (MQ=255)
                                                                            cAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGat                   <  1:1298948/100‑1 (MQ=255)
                                                                             aGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGata                  <  2:1515652/100‑1 (MQ=255)
                                                                               tttttGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGAtata                <  2:2261520/100‑1 (MQ=255)
                                                                                   tGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATTATTc            <  2:4676079/100‑1 (MQ=255)
                                                                                     ttAATGGATACAGCCATCTATTTGTGGAACTTTTCGTAGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGCGCGTTGCCATCAAGATATAATTCGt          <  2:1626979/100‑1 (MQ=14)
                                                                                     ttAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGt          >  2:337431/1‑100 (MQ=255)
                                                                                      tAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTc         >  2:3191182/1‑100 (MQ=255)
                                                                                       aaTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCa        >  2:1817516/1‑100 (MQ=255)
                                                                                       aaTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCa        >  2:1421564/1‑100 (MQ=255)
                                                                                       aaTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCa        >  1:4630517/1‑100 (MQ=255)
                                                                                       aaTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCAGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCa        >  1:5373066/1‑100 (MQ=255)
                                                                                        aTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAt       <  1:2238461/100‑1 (MQ=255)
                                                                                         tGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAtt      >  1:4165207/1‑100 (MQ=255)
                                                                                          ggATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTg     <  1:2319082/100‑1 (MQ=255)
                                                                                          ggATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTg     <  1:5230110/100‑1 (MQ=255)
                                                                                          ggATAAAGCCATCAATTTGTGGACATTTTCGTCGACCCAGGTCTGGCAGAAAGAGGCCAGGTTCTGACGAGCGTAGCCATCAAGATATAATTCGTCATTg     >  1:5173388/1‑100 (MQ=17)
                                                                                           gATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGa    <  2:4230255/100‑1 (MQ=255)
                                                                                            aTAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGAt   <  2:2534553/100‑1 (MQ=255)
                                                                                             tAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATa  >  2:4349273/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATA  >  NC_000913/1571719‑1571911

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