Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,361,104 T→C intergenic (+488/‑72) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,361,1040TC97.7% 582.1 / ‑0.6 171intergenic (+488/‑72)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (77/90);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 171 total aligned reads displayed.

CCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAA  >  NC_000913/3361007‑3361200
                                                                                                 |                                                                                                
   gCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCcaca                                                                                              <  1:3317732/99‑1 (MQ=255)
     aCCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAAc                                                                                           <  1:1210338/100‑1 (MQ=255)
      ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACa                                                                                          <  1:2174152/100‑1 (MQ=255)
      ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACa                                                                                          <  1:2452148/100‑1 (MQ=255)
      ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACa                                                                                          >  1:1948864/1‑100 (MQ=255)
      ccTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATCTCTGAACGAGTTATTTATTTACATCCACAAACa                                                                                          >  2:1558805/1‑100 (MQ=255)
       cTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAg                                                                                         >  1:2631181/1‑100 (MQ=255)
       cTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAg                                                                                         >  1:787680/1‑100 (MQ=255)
         tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACa                                                                                          <  2:2028523/97‑1 (MQ=255)
         tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGat                                                                                       >  2:140184/1‑100 (MQ=255)
         tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGat                                                                                       >  2:685257/1‑100 (MQ=255)
         tattatTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGat                                                                                       >  2:1698965/1‑100 (MQ=255)
          attattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGata                                                                                      >  2:2877862/1‑100 (MQ=255)
           ttattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatat                                                                                     >  2:2916296/1‑100 (MQ=255)
            tattCACATGATCGTAACTATAACATTCATTAAATAAATACATCACAGCACGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatata                                                                                    <  2:1984867/100‑1 (MQ=255)
            tattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatata                                                                                    <  2:1090159/100‑1 (MQ=255)
            tattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatata                                                                                    <  1:425075/100‑1 (MQ=255)
            tattCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatata                                                                                    <  1:1478323/100‑1 (MQ=255)
               ctacaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGatatatat                                                                                 <  1:3238686/98‑1 (MQ=255)
                cacaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGCtatatata                                                                                >  1:2397303/1‑100 (MQ=255)
                 acaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGAtatatatag                                                                               >  2:3238673/1‑99 (MQ=255)
                 acaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATAc                                                                               <  1:982427/100‑1 (MQ=255)
                 acaTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATAc                                                                               <  1:173113/100‑1 (MQ=255)
                  caTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACt                                                                              >  1:1051252/1‑100 (MQ=255)
                  caTGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACt                                                                              <  1:250790/100‑1 (MQ=255)
                    tGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGc                                                                            >  2:2158638/1‑100 (MQ=255)
                       tCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTc                                                                         <  2:3328117/100‑1 (MQ=255)
                        cGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCa                                                                        >  1:2177343/1‑100 (MQ=255)
                        cGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCa                                                                        <  1:1799810/100‑1 (MQ=255)
                          tAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGc                                                                      >  2:71390/1‑100 (MQ=255)
                            aCTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCtt                                                                    >  2:402224/1‑100 (MQ=255)
                             cTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCttt                                                                   >  2:2603511/1‑100 (MQ=255)
                             cTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCttt                                                                   >  1:1590294/1‑100 (MQ=255)
                              tataAAATTCATTAAATAAATACATCCCTGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCtttt                                                                  <  2:1173392/100‑1 (MQ=255)
                               ataAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCttttt                                                                 <  2:245319/100‑1 (MQ=255)
                                taAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTtttta                                                                >  1:253379/1‑100 (MQ=255)
                                 aaaaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttat                                                               >  2:2559760/1‑100 (MQ=255)
                                 aaaaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttat                                                               <  1:1493939/100‑1 (MQ=255)
                                  aaaTTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCCGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTttttatt                                                              >  2:1431460/1‑100 (MQ=255)
                                       cATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaa                                                         >  2:3028873/1‑100 (MQ=255)
                                       cATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaa                                                         <  2:1807198/100‑1 (MQ=255)
                                        atgaaataaatACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                        <  2:3248377/97‑1 (MQ=255)
                                        aTTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAAAAGATATATATACTGCTTCTGCTTTTTATTTTaaaa                                                        >  2:3228784/1‑100 (MQ=255)
                                        aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                        <  1:2598731/100‑1 (MQ=255)
                                        aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                        <  1:2599145/100‑1 (MQ=255)
                                        aTTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaa                                                        <  2:2307682/100‑1 (MQ=255)
                                          taaataaaTACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaaa                                                      <  1:550994/100‑1 (MQ=255)
                                          taaataaaTACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaaa                                                      >  1:250101/1‑100 (MQ=255)
                                           aaataaatACATCCCATCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaat                                                     <  2:1492769/100‑1 (MQ=255)
                                             ataaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataa                                                   >  1:1881263/1‑100 (MQ=255)
                                              taaataCATCCCAGCAAGTTGGCTATTTCTGAACTAGTTATTTATTTACATCCACTAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataaa                                                  <  1:1642840/100‑1 (MQ=255)
                                              taaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataaa                                                  >  1:3118149/1‑100 (MQ=255)
                                               aaataCATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataaat                                                 <  2:160307/100‑1 (MQ=255)
                                                    cATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaaagatcggaaga                                            >  1:289369/1‑90 (MQ=255)
                                                    cATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGtt                                            >  1:1146318/1‑100 (MQ=255)
                                                    cATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGtt                                            <  1:453429/100‑1 (MQ=255)
                                                     aTCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGtta                                           <  2:1948860/100‑1 (MQ=255)
                                                       cccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatt                                         >  2:2764620/1‑100 (MQ=255)
                                                       cccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatt                                         >  1:3321280/1‑100 (MQ=255)
                                                        ccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatta                                        >  1:1005798/1‑100 (MQ=255)
                                                        ccAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatta                                        <  1:1483419/100‑1 (MQ=255)
                                                          aGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAg                                      <  2:308961/100‑1 (MQ=255)
                                                           gCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTaa                                       <  2:1367298/98‑1 (MQ=255)
                                                           gCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAgg                                     <  1:3395688/100‑1 (MQ=255)
                                                            cAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGa                                    <  2:3118137/100‑1 (MQ=255)
                                                                ttGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATATTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATg                                <  2:1522137/100‑1 (MQ=255)
                                                                ttGGCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATg                                >  2:2888180/1‑100 (MQ=255)
                                                                 tcttccgatCTGAACGAGTTATTTATTTACATCCACAAACAAATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGt                               <  1:2842636/92‑1 (MQ=255)
                                                                  ggCTATTTCTGAATGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGtt                              >  2:421446/1‑100 (MQ=255)
                                                                  ggCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGtt                              <  1:2877873/100‑1 (MQ=255)
                                                                   gCTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTc                             <  1:257163/100‑1 (MQ=255)
                                                                    cTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGtt                              <  2:2722060/98‑1 (MQ=255)
                                                                    cTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                            >  2:1756720/1‑100 (MQ=255)
                                                                    cTATTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                            <  1:2764631/100‑1 (MQ=255)
                                                                      aTTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAGAAAATAAATTTGTTATTAAGGATATGTTCata                          <  2:2631173/100‑1 (MQ=255)
                                                                      aTTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                          <  1:434625/100‑1 (MQ=255)
                                                                      aTTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                          <  1:71389/100‑1 (MQ=255)
                                                                      aTTTCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                          <  1:343107/100‑1 (MQ=255)
                                                                       tttCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCatat                         <  2:659144/100‑1 (MQ=255)
                                                                        ttCTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATg                        >  1:3005689/1‑100 (MQ=255)
                                                                          cTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtt                      >  2:2730613/1‑100 (MQ=255)
                                                                          cTGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtn                      <  1:946264/100‑2 (MQ=255)
                                                                           tGAACGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttt                     >  1:569569/1‑100 (MQ=255)
                                                                             aaCGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttttt                   <  1:2005367/100‑1 (MQ=255)
                                                                              atgAGTTATTTATATACATCCACAAACAGATGTATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTc                  >  1:2289059/3‑100 (MQ=255)
                                                                              aCGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTc                  <  1:2399551/100‑1 (MQ=255)
                                                                              aCGAGTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTc                  <  1:685257/100‑1 (MQ=255)
                                                                                  gTTATTTATTTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaa              <  2:3349395/100‑1 (MQ=255)
                                                                                    tatttattTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAg            >  2:223519/1‑100 (MQ=255)
                                                                                    tatttattTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAg            <  1:140183/100‑1 (MQ=255)
                                                                                    tatttattTACATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAg            >  2:1949825/1‑100 (MQ=255)
                                                                                       ttatttaCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGaan         <  2:2212148/100‑2 (MQ=255)
                                                                                         atttaCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCt       >  2:465438/1‑100 (MQ=255)
                                                                                           ttaCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTca     <  2:2171756/100‑1 (MQ=255)
                                                                                            taCATCCACATACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcac    <  1:1739999/100‑1 (MQ=255)
                                                                                            taCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcac    <  1:771975/100‑1 (MQ=255)
                                                                                            taCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcac    >  1:2160593/1‑100 (MQ=255)
                                                                                            taCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcac    <  2:787680/100‑1 (MQ=255)
                                                                                             aCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcaca   <  2:1051251/100‑1 (MQ=255)
                                                                                             aCATCCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTcaca   <  2:2804498/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAA  >  NC_000913/3361007‑3361200

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: