Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,624,020 G→A intergenic (‑170/+85) dgcZ ← / ← ydeI diguanylate cyclase, zinc‑sensing/hydrogen peroxide resistance OB fold protein; putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,0200GA98.8% 559.2 / NA 160intergenic (‑170/+85)dgcZ/ydeIdiguanylate cyclase, zinc‑sensing/hydrogen peroxide resistance OB fold protein; putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (74/84);  minor base C (1/1);  total (75/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

ATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAA  >  NC_000913/1623924‑1624116
                                                                                                |                                                                                                
aTAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAACAAAGTGTTTAATTCTGTAAATTACCTCTCCAt                                                                                               >  2:4707074/1‑100 (MQ=255)
aTAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACAt                                                                                               <  1:2518148/100‑1 (MQ=255)
aTAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACAt                                                                                               <  2:1514446/100‑1 (MQ=255)
 tAGAAATCCGCAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAGATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACAtt                                                                                              >  2:4678776/1‑100 (MQ=255)
     aaTCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCg                                                                                          <  1:2399391/100‑1 (MQ=255)
     aaTCCAAAGTTATGATACATCCCTGAGAGAACATTGTTATTCAACGCAAAATCAGGGGACAAAAAAGTGTTTAATTATCTAAATTACCTCTCCCTTATCg                                                                                          <  2:4838717/100‑1 (MQ=255)
         cGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGtaaa                                                                                      >  2:293035/1‑100 (MQ=255)
          gAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGtaaat                                                                                     <  1:2697112/100‑1 (MQ=255)
          gAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGtaaat                                                                                     <  1:3461557/100‑1 (MQ=255)
          gAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGtaaat                                                                                     <  2:440637/100‑1 (MQ=255)
               tATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAg                                                                                >  1:1803729/1‑100 (MQ=255)
               tATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAg                                                                                >  1:2214688/1‑100 (MQ=255)
               tATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAg                                                                                <  2:4208671/100‑1 (MQ=255)
               tATGAGTCATCTCGGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAg                                                                                <  2:134881/100‑1 (MQ=255)
                 tGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGa                                                                              >  2:5026174/1‑100 (MQ=255)
                  taGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGAt                                                                             <  1:5066627/99‑1 (MQ=255)
                  gAGTCATCTATGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGAt                                                                             <  1:4194353/100‑1 (MQ=255)
                   aGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGAta                                                                            >  2:4987866/1‑99 (MQ=255)
                   aGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATg                                                                            <  2:1852284/100‑1 (MQ=255)
                    gTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGa                                                                           >  2:576976/1‑100 (MQ=255)
                    gTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGa                                                                           <  1:1078751/100‑1 (MQ=255)
                      cATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACa                                                                         >  1:810027/1‑100 (MQ=255)
                      cATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACa                                                                         >  1:1464906/1‑100 (MQ=255)
                         ctctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAAt                                                                      <  1:2634562/100‑1 (MQ=255)
                         ctctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAAt                                                                      >  1:3022810/1‑100 (MQ=255)
                         ctctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAAt                                                                      >  1:4420229/1‑100 (MQ=255)
                         atctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAGTAAAAGGATGACAAAt                                                                      <  2:1381045/99‑1 (MQ=255)
                          tctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATa                                                                     >  2:48013/1‑100 (MQ=255)
                          tctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATa                                                                     >  1:3123252/1‑100 (MQ=255)
                          tctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATa                                                                     >  2:204441/1‑100 (MQ=255)
                          tctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATa                                                                     >  2:2550326/1‑100 (MQ=255)
                           ctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAg                                                                    >  2:1707584/1‑100 (MQ=255)
                           ctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAg                                                                    <  1:4754762/100‑1 (MQ=255)
                            tGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGc                                                                   <  2:3390496/100‑1 (MQ=255)
                             gagaTAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCa                                                                  <  1:1637895/100‑1 (MQ=255)
                             gagaTAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATGCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCa                                                                  <  2:3515609/100‑1 (MQ=255)
                                 tAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAAc                                                              >  2:4740322/1‑100 (MQ=255)
                                  aaCATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAAcc                                                             <  2:417660/100‑1 (MQ=255)
                                   aCATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAAccc                                                            >  2:13829/1‑100 (MQ=255)
                                   aCATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAAccc                                                            <  1:1490501/100‑1 (MQ=255)
                                      ttGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAAt                                                         <  1:3718318/100‑1 (MQ=255)
                                       tgtgATTTACAACAAAATAAGCGGATAAACAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATa                                                        <  1:2826527/100‑1 (MQ=255)
                                        gtgATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATAc                                                       <  1:2835887/100‑1 (MQ=255)
                                           aTTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCt                                                    >  1:3911115/1‑100 (MQ=255)
                                           aTTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCt                                                    >  2:1221680/1‑100 (MQ=255)
                                            tttAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTa                                                   >  1:1435209/1‑100 (MQ=255)
                                             ttAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTaa                                                  >  1:1898543/1‑100 (MQ=255)
                                               aaaaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATg                                                <  2:2184947/100‑1 (MQ=255)
                                               aaaaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATg                                                <  1:1372353/100‑1 (MQ=255)
                                               aaaaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATg                                                <  2:595674/100‑1 (MQ=255)
                                                aaaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATgg                                               <  1:4819929/100‑1 (MQ=255)
                                                 aaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGc                                              <  2:348656/100‑1 (MQ=255)
                                                 aaCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGc                                              >  1:1981144/1‑100 (MQ=255)
                                                  aCAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTAACTCTACATTATCGTAAATAAAATAATAACAAATAGCATAACCCAATACCCTAATGGcc                                             >  2:5056331/1‑100 (MQ=255)
                                                   cAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGccc                                            >  2:3248188/1‑100 (MQ=255)
                                                   cAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGccc                                            >  2:1833914/1‑100 (MQ=255)
                                                   cAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGccc                                            <  1:4178187/100‑1 (MQ=255)
                                                    aaaaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCa                                           <  1:2010849/100‑1 (MQ=255)
                                                     aaaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                          <  1:5105262/100‑1 (MQ=255)
                                                     aaaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                          >  1:5058397/1‑100 (MQ=255)
                                                     aaaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                          <  1:48689/100‑1 (MQ=255)
                                                      aaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCagt                                         <  1:205089/100‑1 (MQ=255)
                                                      aaTCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCagt                                         >  1:3379647/1‑100 (MQ=255)
                                                        tCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCagtag                                       <  2:23804/100‑1 (MQ=255)
                                                            cGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCa                                   >  1:3281324/1‑100 (MQ=255)
                                                              gATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAgg                                 <  2:2974662/100‑1 (MQ=255)
                                                              gATAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAgg                                 >  1:5132086/1‑100 (MQ=255)
                                                               aTAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGc                                <  1:2617694/100‑1 (MQ=255)
                                                               aTAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGc                                <  2:1780197/100‑1 (MQ=255)
                                                                tAAAAAAGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGcc                               <  1:1234149/100‑1 (MQ=255)
                                                                  aaaaaGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGTATAACCCAATACCCTAATGGCCCATTAGTTCAGGCCAt                             >  2:369537/1‑100 (MQ=255)
                                                                   aaaaGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATc                            <  1:612319/100‑1 (MQ=255)
                                                                    aaaGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCa                           >  2:2763715/1‑100 (MQ=255)
                                                                     aaGTGTTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAg                          <  2:1321572/100‑1 (MQ=255)
                                                                     aaGTGCTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAg                          >  2:313374/1‑100 (MQ=255)
                                                                       gtgtTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATCGCATACCCCAATACCCTAGTGGCCGAGTAGTTCAGGCCATCAGGc                        >  1:4826938/1‑100 (MQ=255)
                                                                       gtgtTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGc                        <  1:2406194/100‑1 (MQ=255)
                                                                       gtgtTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGc                        <  1:679063/100‑1 (MQ=255)
                                                                         gtTTAATTCTGTAAATTACCTCTACGTTATCGTAAATAAAAGGATGAAAAATAGCATAACCCAATACCCTAATGACCCAGTAGTTCAGGCCATCATGCTa                      >  1:4615794/1‑100 (MQ=255)
                                                                         gtTTAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTa                      <  1:366337/100‑1 (MQ=255)
                                                                            tAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatt                   >  2:1736772/1‑100 (MQ=255)
                                                                            tAATTCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatt                   >  1:2105300/1‑100 (MQ=255)
                                                                                tCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatttatt               <  1:676361/100‑1 (MQ=255)
                                                                                tCTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatttatt               <  1:3870140/100‑1 (MQ=255)
                                                                                 cTGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatttattt              >  2:664048/1‑100 (MQ=255)
                                                                                  tGTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttt             <  1:1036470/100‑1 (MQ=255)
                                                                                   gTAAATTACCTCTACATTATCGTAAATAACAGGATGTCAAATATCATAACCCAATACCCTCGAGGCCCAGTAGTTCAGGCCATCAGGCTATtttattttt            <  2:385662/100‑1 (MQ=255)
                                                                                   gTAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttattttt            >  2:4770681/1‑100 (MQ=255)
                                                                                    tAAATTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttattttta           <  1:3144376/100‑1 (MQ=255)
                                                                                     aaaTTACCTCTACATTACCGTAAATAAAAGGATGACAAATAGCAGAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttttat          <  2:816691/100‑1 (MQ=255)
                                                                                      aaTTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttttatt         <  2:2711431/100‑1 (MQ=255)
                                                                                      aaTTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttttatt         >  1:788275/1‑100 (MQ=255)
                                                                                       aTTACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATttatttttattt        <  2:2064665/100‑1 (MQ=255)
                                                                                        ttACCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTc       <  2:732075/100‑1 (MQ=255)
                                                                                          aCCTCTACATTATCGTAAATAAAAGGATGACAAATAGCGTAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTg     <  2:3342277/100‑1 (MQ=255)
                                                                                          aCCTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTg     <  1:2526308/100‑1 (MQ=255)
                                                                                           ccTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGc    <  1:3284657/100‑1 (MQ=255)
                                                                                           ccTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGc    <  2:4859998/100‑1 (MQ=255)
                                                                                           ccTCTACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGACCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGc    >  1:4505990/1‑100 (MQ=255)
                                                                                             tctACATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCaa  >  1:3524242/1‑100 (MQ=255)
                                                                                                |                                                                                                
ATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAA  >  NC_000913/1623924‑1624116

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