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,8160TG97.9% 470.5 / ‑2.1 142D77A (GAT→GCT) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (62/77);  total (64/78)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

ATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATA  >  NC_000913/1571719‑1571911
                                                                                                 |                                                                                               
aTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCCTTGTCGATCCAGTTTTTGTTAATGGATAAAGcc                                                                                               <  2:365431/100‑1 (MQ=255)
 tGCCACAGCTCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTCATGGATAAAGCCa                                                                                              <  2:2570168/100‑1 (MQ=255)
  gCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCAt                                                                                             >  2:2435198/1‑100 (MQ=255)
    cacaGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCa                                                                                           <  1:818806/100‑1 (MQ=255)
      caGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTACGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAAt                                                                                         >  1:4036109/1‑100 (MQ=255)
          tCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgt                                                                                     >  1:166513/1‑100 (MQ=255)
          tCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgt                                                                                     <  1:1637712/100‑1 (MQ=255)
           cGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgtg                                                                                    <  1:37678/100‑1 (MQ=255)
           cGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTtgtg                                                                                    >  2:2054863/1‑100 (MQ=255)
            ggCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTgg                                                                                   <  1:2473239/100‑1 (MQ=255)
            ggCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTgg                                                                                   <  1:1861184/100‑1 (MQ=255)
            ggCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTgg                                                                                   >  1:1428464/1‑100 (MQ=255)
             gCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGa                                                                                  <  1:3963729/100‑1 (MQ=255)
              cgaCCATATTTACGCAACGCAGGTCGGTGGCTGCGGATTGCGGATATTCTTCTTGGTCGGTCCAGTTGTTGTTAAGGGGGAAAGCCATCAATTTGTGTAc                                                                                 <  2:3975314/98‑1 (MQ=14)
               aaCCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACa                                                                                <  2:2756052/100‑1 (MQ=255)
                   atatTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACAtttt                                                                            <  1:802854/100‑1 (MQ=255)
                    tatTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtc                                                                           <  1:363765/100‑1 (MQ=255)
                     atTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcg                                                                          >  2:1703365/1‑100 (MQ=255)
                     atTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcg                                                                          <  2:3931948/100‑1 (MQ=255)
                     atTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcg                                                                          >  1:569457/1‑100 (MQ=255)
                     atTTACGCAACGCAGGTCGATGGCTGCGAATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcg                                                                          <  2:1532794/100‑1 (MQ=255)
                      tttACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTtcgt                                                                         >  1:3200302/1‑100 (MQ=255)
                        tACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGAGCCAGTTTTTGTTAGTGGATCAAGCCATCAATTTGTGGACATTTtcgtcg                                                                       >  2:3941112/1‑100 (MQ=17)
                         acgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTcgtcgt                                                                      >  1:3316924/1‑100 (MQ=255)
                          cgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCgtcgtc                                                                     <  1:138974/100‑1 (MQ=255)
                           gctacgcaGTTCTATGGCTTCTGATTGCGAATATTCTTATTTGTCGATCCACTTTTTTTTCCTGCATAAAGCCATCAATGTGTGGACATTTTCGTCGTcc                                                                    <  2:3983955/97‑1 (MQ=255)
                           gcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCgtcgtct                                                                    <  2:3124206/100‑2 (MQ=255)
                           gcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTcc                                                                    <  1:2325338/100‑1 (MQ=255)
                            caacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTccc                                                                   >  2:3261709/1‑100 (MQ=255)
                              acgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAg                                                                 <  2:3158213/100‑1 (MQ=255)
                              acgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAg                                                                 >  1:3207419/1‑100 (MQ=255)
                               cgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAgg                                                                >  1:1562967/1‑100 (MQ=255)
                               cgcaGGACCATGGCGGCGGATTGCCGAGATTCTTCCTGGGCGATCCAGTTTTTGTTAATGAATAAAGCCATCCATTTGTGGACATTTTCGTCGTCCCAgg                                                                <  1:3891832/100‑1 (MQ=17)
                                   ggTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCAATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg                                                            >  2:3909346/1‑100 (MQ=255)
                                     tCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCCGTTGTTGTTAATGGATAATGCCATCAATTTGTGGACATTTTCGTCGTCCCGGGTCTGGc                                                          <  2:766481/100‑1 (MQ=17)
                                     tCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATCGATAAAGCCATCAATTTGTGGACATTTTCGTCGTTCCAGGTCTgac                                                          <  2:3708391/100‑3 (MQ=14)
                                     tCGATGGCTGCGGATTGCAGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGc                                                          <  2:1217294/100‑1 (MQ=255)
                                      cGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCa                                                         >  1:3230399/1‑100 (MQ=255)
                                      cGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCa                                                         <  2:2189396/100‑1 (MQ=255)
                                       gATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAg                                                        >  2:2910351/1‑100 (MQ=255)
                                        aTGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGa                                                       >  1:2873470/1‑100 (MQ=255)
                                           gCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAg                                                    >  1:4046243/1‑100 (MQ=255)
                                            cTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTTATGGATAAAGCCATCAATTTGTGGACTTTTTCGTCGTCCCAGTTCTTGCAGAAAGt                                                   <  2:3301177/100‑1 (MQ=14)
                                            cTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGt                                                   >  2:3431486/1‑100 (MQ=255)
                                               cggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGc                                                <  1:3119632/100‑1 (MQ=255)
                                                ggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGcc                                               <  1:4005888/100‑1 (MQ=255)
                                                 gattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCa                                              >  1:730687/1‑100 (MQ=255)
                                                    tgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGt                                           <  1:1110132/100‑1 (MQ=255)
                                                    tgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGt                                           <  2:2218599/100‑1 (MQ=255)
                                                    tgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGt                                           >  2:2005526/1‑100 (MQ=255)
                                                     gcggatATTCTTCTTTGTCGATCCAGTTTTTTTTTTTTTTTAAAGCCATCAATTTGTGGACATTTTCGTCTTTTCAGGTTTGGCAGAAAGTGGTCAGGtt                                          >  2:1633423/1‑100 (MQ=14)
                                                     gcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAACCTGGCCacttt                                          >  2:3352752/1‑96 (MQ=18)
                                                     gcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGtt                                          >  1:1908079/1‑100 (MQ=255)
                                                     gcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGtt                                          >  2:548274/1‑100 (MQ=255)
                                                      cggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTc                                         >  1:1278233/1‑100 (MQ=255)
                                                       ggctatTCTTTTTTATCGATCCAGTTTGTGTTCATGGATAAAGCCATGATCTTAGTGACATTGTCGTCGTCTCAGGTTTGGGAGAAAGTGGCCAGGTTCt                                        <  2:134532/97‑1 (MQ=9)
                                                       ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCt                                        >  2:157795/1‑100 (MQ=255)
                                                       ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCt                                        <  2:2328395/100‑1 (MQ=255)
                                                         atatTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGa                                      <  2:1332601/100‑1 (MQ=255)
                                                         atatTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGa                                      <  2:3082515/100‑1 (MQ=255)
                                                          tatTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGAc                                     <  1:3591750/100‑1 (MQ=255)
                                                          tatTCTTCTGTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGAc                                     <  1:187201/100‑1 (MQ=255)
                                                           atTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACg                                    >  2:1453667/1‑100 (MQ=255)
                                                           atTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACg                                    >  1:1769275/1‑100 (MQ=255)
                                                           atTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACg                                    <  1:3829/100‑1 (MQ=255)
                                                             tcttctTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGac                                  >  2:2504548/1‑99 (MQ=255)
                                                             tcttctTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAg                                  <  1:1652090/100‑1 (MQ=255)
                                                             tcttctTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAg                                  <  1:1479371/100‑1 (MQ=255)
                                                             tcttctTTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAg                                  <  1:2037418/100‑1 (MQ=255)
                                                                 cttTGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGtt                              <  2:2900826/100‑1 (MQ=255)
                                                                  tttGTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTg                             <  2:3760611/100‑1 (MQ=255)
                                                                     gTCGATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCa                          >  1:521695/1‑100 (MQ=255)
                                                                        gATCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTAACGAGCGTTGCCATCa                       >  1:3353083/1‑100 (MQ=255)
                                                                         aTCCAGTTTTTGTTTTTTTTTAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCaa                      >  2:1734967/1‑100 (MQ=14)
                                                                         aTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCaa                      <  2:1622233/100‑1 (MQ=255)
                                                                         aTCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCaa                      <  2:2589336/100‑1 (MQ=255)
                                                                          tCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAg                     >  1:4013142/1‑100 (MQ=255)
                                                                          tCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAg                     >  1:3099381/1‑100 (MQ=255)
                                                                          tCCAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAg                     <  1:1237797/100‑1 (MQ=255)
                                                                           ccAGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGa                    <  1:2817087/100‑1 (MQ=255)
                                                                             aGTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGata                  <  1:322538/100‑1 (MQ=255)
                                                                              gTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGatat                 >  2:1006184/1‑100 (MQ=255)
                                                                              gTTTTTGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGatat                 >  2:359350/1‑100 (MQ=255)
                                                                               tttttGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGAtata                <  2:1882607/100‑1 (MQ=255)
                                                                                ttttGTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATaa               >  1:2084323/1‑100 (MQ=255)
                                                                                    gTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCg           <  2:1068506/100‑1 (MQ=255)
                                                                                    gTTAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCg           >  1:2179391/1‑100 (MQ=255)
                                                                                     ttAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGt          <  2:1769470/100‑1 (MQ=255)
                                                                                      tAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTc         >  2:2340518/1‑100 (MQ=255)
                                                                                      tAATGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTc         <  1:2089600/100‑1 (MQ=255)
                                                                                       aaTGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCa        >  2:149991/1‑100 (MQ=255)
                                                                                        aTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAt       <  1:2818274/100‑1 (MQ=255)
                                                                                        aTGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAt       <  1:1172399/100‑1 (MQ=255)
                                                                                         tGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAtt      <  1:1594075/100‑1 (MQ=255)
                                                                                         tGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAtt      >  1:374674/1‑100 (MQ=255)
                                                                                         tGGATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCAtt      <  1:1103610/100‑1 (MQ=255)
                                                                                          ggATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATCATTCGTCATTg     >  2:299744/1‑100 (MQ=255)
                                                                                           gATAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGa    >  1:2040093/1‑100 (MQ=255)
                                                                                            aTAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGAt   <  1:3534684/100‑1 (MQ=255)
                                                                                             tAAAGCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATa  >  1:358619/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATA  >  NC_000913/1571719‑1571911

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: